A346100 a(n) = A100995(gcd(n, A064989(A319626(A324886(n))))).

0, 1, 1, 2, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 1, 1, 0, 1, 1, 0, 0, 0, 1, 1, 1, 2, 1, 1, 1, 1, 1, 0, 1, 1, 0, 2, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1, 2, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 1, 1, 3

Offset: 1

Author

Antti Karttunen, Jul 07 2021

Keywords

• nonn

Crossrefs

Cf. A064989, A100995, A319626, A324886, A346087, A346095, A346096, A346098, A346099.

Cf. A346090 (positions of 0's).

Programs

PARI

A064989(n) = { my(f = factor(n)); if((n>1 && f[1,1]==2), f[1,2] = 0); for (i=1, #f~, f[i,1] = precprime(f[i,1]-1)); factorback(f); };
A319626(n) = (n / gcd(n, A064989(n)));
A346100(n) = isprimepower(gcd(n, A064989(A319626(A324886(n))))); \\ Rest of program given in A324886.

Formula

a(n) = A100995(A346099(n)) = A100995(gcd(n, A064989(A319626(A324886(n))))).