|
%I #7 Dec 28 2020 17:48:06
%S 1,1,1,1,1,2,1,1,1,4,1,2,1,6,4,1,1,4,1,4,3,10,1,2,2,12,4,2,1,8,1,1,5,
%T 16,12,4,1,18,12,4,1,12,1,10,4,22,1,2,3,16,16,4,1,8,20,6,9,28,1,8,1,
%U 30,12,1,3,4,1,16,11,8,1,4,1,36,16,6,15,24,1,4,1,40,1,12,16,42,28,10,1,16,4,22,15,46,36,2,1,36
%N a(n) = A003958(n) / gcd(n-1, A003958(n)).
%H Antti Karttunen, <a href="/A340082/b340082.txt">Table of n, a(n) for n = 1..8191</a>
%H Antti Karttunen, <a href="/A340082/a340082.txt">Data supplement: n, a(n) computed for n = 1..65537</a>
%F a(n) = A003958(n) / A340081(n) = A003958(n) / gcd(n-1, A003958(n)).
%o (PARI)
%o A003958(n) = if(1==n,n,my(f=factor(n)); for(i=1,#f~,f[i,1]--); factorback(f));
%o A340082(n) = { my(u=A003958(n)); u/gcd(n-1, u); };
%Y Cf. A003958, A340081, A340083, A340085 (gives the odd part).
%Y Cf. also A160595, A340072.
%K nonn
%O 1,6
%A _Antti Karttunen_, Dec 28 2020
|
