A346475 a(n) = A342919(A276086(n)).

0, 1, 1, 5, 1, 7, 1, 7, 1, 31, 13, 41, 1, 1, 11, 37, 2, 47, 1, 11, 7, 43, 19, 53, 2, 13, 17, 49, 11, 59, 1, 3, 5, 41, 17, 55, 1, 59, 71, 247, 53, 317, 19, 73, 23, 289, 127, 359, 13, 29, 113, 331, 37, 401, 11, 101, 67, 373, 169, 443, 1, 11, 13, 47, 5, 61, 17, 23, 43, 277, 121, 347, 1, 83, 107, 319, 71, 389, 31, 97, 2, 361, 163

Offset: 0

Author

Antti Karttunen, Jul 21 2021

Keywords

• nonn
• base
• look

Comments

For n >= 1, each term a(n) is a divisor of A342002(n).

Links

• Antti Karttunen, Table of n, a(n) for n = 0..13860
• Index entries for sequences related to primorial base

Crossrefs

Cf. A001615, A003415, A276086, A327860, A342919, A346470, A346471.

Cf. also A342002, A345930, A346474 for sequences with similar scatter plots.

Programs

PARI

A001615(n) = if(1==n,n, my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[i, 2] + f[i, 1]^(f[i, 2]-1))); \\ After code in A001615
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A342919(n) = { my(u=A003415(n)); (u/gcd(u, A001615(n))); };
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A346475(n) = A342919(A276086(n));

Formula

a(n) = A342919(A276086(n)).

a(n) = A327860(n) / gcd(A327860(n), A346470(n)).