A353574 a(n) = A342002(n) / gcd(A342002(n), A342002(A276154(n))).

1, 1, 5, 1, 7, 1, 7, 2, 31, 13, 41, 1, 9, 11, 37, 2, 47, 1, 11, 7, 43, 19, 53, 1, 13, 17, 49, 11, 59, 1, 9, 5, 41, 17, 55, 2, 59, 71, 247, 53, 317, 19, 73, 46, 289, 127, 359, 13, 87, 113, 331, 74, 401, 11, 101, 67, 373, 169, 443, 1, 11, 13, 47, 5, 61, 17, 69, 43, 277, 121, 347, 2, 83, 107, 319, 71, 389, 31, 97, 16

Offset

1,3

Comments

Compare the scatter plot to that of A342002.

Links

Antti Karttunen, Table of n, a(n) for n = 1..11550

Index entries for sequences computed from indices in prime factorization

Index entries for sequences related to primorial base

Formula

a(n) = A342002(n) / A353573(n) = A342002(n) / gcd(A342002(n), A353572(n)).

Prog

(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003557(n) = (n/factorback(factorint(n)[, 1]));
A003961(n) = { my(f = factor(n)); for(i=1, #f~, f[i, 1] = nextprime(f[i, 1]+1)); factorback(f); }; \\ From A003961
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A342001(n) = (A003415(n) / A003557(n));
A353571(n) = A342001(A003961(n));
A353574(n) = { my(u=A342001(A276086(n))); u/gcd(u, A353571(A276086(n))); };

Crossrefs

Cf. A003415, A003557, A003961, A342001, A342002, A353572, A353573, A353575.

Sequence in context: A318676 A346475 A346474 * A265293 A348502 A340089

Adjacent sequences: A353571 A353572 A353573 * A353575 A353576 A353577

Keyword

nonn,base,easy,look

Author

Antti Karttunen, Apr 29 2022

Status

approved