A318829 - OEIS

%I #11 Sep 09 2018 19:14:28

%S 1,1,1,1,1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,

%T 2,1,1,1,2,1,1,1,1,1,2,1,1,1,1,1,2,1,1,1,2,1,1,1,1,1,1,1,2,1,1,1,1,1,

%U 1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,4,1,2,1,1,1,2,1,1,1,2,1,1,1,2,1,1,1,1,1,2

%N a(n) = A063994(n) / A049559(n) = (1/gcd(n-1, phi(n))) * Product_{primes p dividing n} gcd(p-1, n-1).

%C Records occur at: 1, 15, 85, 247, 671, 949, 1105, 1387, 2047, 2821, 9471, 11305, 13747, 13981, 29341, 40885, 51319, 63973, ...

%H Antti Karttunen, <a href="/A318829/b318829.txt">Table of n, a(n) for n = 1..65537</a>

%F a(n) = A063994(n) / A049559(n).

%F a(n) = A160595(n) / A247074(n).

%o (PARI)

%o A049559(n) = gcd(eulerphi(n), n-1); \\ From A049559

%o A063994(n) = { my(f=factor(n)[,1]); prod(i=1, #f, gcd(f[i]-1, n-1)); };

%o A318829(n) = (A063994(n)/A049559(n));

%Y Cf. A049559, A063994, A160595, A247074, A318828.

%K nonn

%O 1,15

%A _Antti Karttunen_, Sep 09 2018