A344757 a(n) = A069359(n) / gcd(A003415(n), A069359(n)).

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 5, 1, 7, 1, 1, 1, 5, 1, 1, 1, 9, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 7, 1, 1, 1, 13, 8, 1, 1, 5, 1, 7, 1, 15, 1, 5, 1, 9, 1, 1, 1, 31, 1, 1, 10, 1, 1, 1, 1, 19, 1, 1, 1, 5, 1, 1, 8, 21, 1, 1, 1, 7, 1, 1, 1, 41, 1, 1, 1, 13, 1, 31, 1, 25, 1, 1, 1, 5, 1, 9, 14, 1, 1, 1, 1, 15

Offset: 2

Author

Antti Karttunen, May 28 2021

Keywords

• nonn

Links

• Antti Karttunen, Table of n, a(n) for n = 2..12012
• Antti Karttunen, Data supplement: n, a(n) computed for n = 2..65537

Crossrefs

Cf. A003415, A069359, A340070, A344756.

Cf. also A344697.

Programs

PARI

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A069359(n) = (n*sumdiv(n, d, isprime(d)/d)); \\ From A069359
A344757(n) = { my(u=A069359(n)); (u/gcd(u,A003415(n))); };

Formula

a(n) = A069359(n) / A340070(n) = A069359(n) / gcd(A003415(n), A069359(n)).