login
a(n) = A008472(n) / gcd(A007947(n), A008472(n)).
3

%I #9 Feb 01 2021 19:41:08

%S 0,1,1,1,1,5,1,1,1,7,1,5,1,9,8,1,1,5,1,7,10,13,1,5,1,15,1,9,1,1,1,1,

%T 14,19,12,5,1,21,16,7,1,2,1,13,8,25,1,5,1,7,20,15,1,5,16,9,22,31,1,1,

%U 1,33,10,1,18,8,1,19,26,1,1,5,1,39,8,21,18,3,1,7,1,43,1,2,22,45,32,13,1,1,20,25,34,49,24,5,1,9

%N a(n) = A008472(n) / gcd(A007947(n), A008472(n)).

%H Antti Karttunen, <a href="/A340678/b340678.txt">Table of n, a(n) for n = 1..16384</a>

%H Antti Karttunen, <a href="/A340678/a340678.txt">Data supplement: n, a(n) computed for n = 1..65537</a>

%F a(n) = A008472(n) / A099636(n) = A008472(n) / gcd(A007947(n), A008472(n)).

%o (PARI)

%o A007947(n) = factorback(factorint(n)[, 1]);

%o A008472(n) = vecsum(factor(n)[, 1]);

%o A340678(n) = { my(v=A008472(n)); (v/gcd(A007947(n), v)); };

%Y Cf. A007947, A008472, A099636, A340677.

%Y Cf. also A082343.

%K nonn

%O 1,6

%A _Antti Karttunen_, Feb 01 2021