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,4

Links

Table of n, a(n) for n=1..108.

Formula

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

Prog

(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.

Crossrefs

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

Cf. A346090 (positions of 0's).

Sequence in context: A005091 A353455 A276516 * A253638 A337586 A345006

Adjacent sequences: A346097 A346098 A346099 * A346101 A346102 A346103

Keyword

nonn

Author

Antti Karttunen, Jul 07 2021

Status

approved