A346100 - OEIS

%I #6 Jul 19 2021 01:26:11

%S 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,

%T 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,

%U 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

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

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

%o (PARI)

%o 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); };

%o A319626(n) = (n / gcd(n, A064989(n)));

%o A346100(n) = isprimepower(gcd(n, A064989(A319626(A324886(n))))); \\ Rest of program given in A324886.

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

%Y Cf. A346090 (positions of 0's).

%K nonn

%O 1,4

%A _Antti Karttunen_, Jul 07 2021