%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