login
a(n) = A069359(n) / gcd(A003415(n), A069359(n)).
3

%I #10 Jun 30 2021 19:47:32

%S 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,

%T 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,

%U 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

%N a(n) = A069359(n) / gcd(A003415(n), A069359(n)).

%H Antti Karttunen, <a href="/A344757/b344757.txt">Table of n, a(n) for n = 2..12012</a>

%H Antti Karttunen, <a href="/A344757/a344757.txt">Data supplement: n, a(n) computed for n = 2..65537</a>

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

%o (PARI)

%o A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));

%o A069359(n) = (n*sumdiv(n, d, isprime(d)/d)); \\ From A069359

%o A344757(n) = { my(u=A069359(n)); (u/gcd(u,A003415(n))); };

%Y Cf. A003415, A069359, A340070, A344756.

%Y Cf. also A344697.

%K nonn

%O 2,11

%A _Antti Karttunen_, May 28 2021