2,2

The p-part of a number k is the highest power of p that divides k. For example, the 2-part of 24 is 8, the 3-part is 3.

One can think of A360519 as a chain of circles, each circle linked to its neighbors to the left and the right. The n-th term b(n) = A360519(n) is a product L(n)*M(n)*R(n), where L(n) is the part of b(n) sharing primes with the term to the left, R(n) the part sharing primes with the term to the right, and M(n) is the rest of b(n).

By definition of A360519, the set of primes in L(n) is disjoint from the primes in R(n).

N. J. A. Sloane, Table of n, a(n) for n = 2..20000

# Suppose bW is a list of the terms of A360519.

# Then f3(bW[n-1], bW[n], bW[n+1]); returns [L(n), M(n), R(n)] where:

with(numtheory);

f3:=proc(a, b, c)

local lefta, midb, rightc, i, p, pa, pc, ta, tb, tc, t1, t2;

ta:=a; tb:=b; tc:=c;

# left

t1:=igcd(a, b);

t2:=factorset(t1);

t2:=convert(t2, list);

lefta:=1;

for i from 1 to nops(t2) do

p:=t2[i];

while (tb mod p) = 0 do lefta:=lefta*p; tb:=tb/p; od;

od:

# right

t1:=igcd(b, c);

t2:=factorset(t1);

t2:=convert(t2, list);

rightc:=1;

for i from 1 to nops(t2) do

p:=t2[i];

while (tb mod p) = 0 do rightc:=rightc*p; tb:=tb/p; od;

od:

# middle

midb:=b/(lefta*rightc);

[lefta, midb, rightc];

end; # N. J. A. Sloane, Mar 09 2023

Cf. A360519, A361118, A361129, A361130.

Sequence in context: A286362 A361118 A361130 * A309200 A361170 A097964

Adjacent sequences: A361125 A361126 A361127 * A361129 A361130 A361131

Scott R. Shannon, Rémy Sigrist, and N. J. A. Sloane, Mar 09 2023

approved