OFFSET
1,3
COMMENTS
a(n) = A198772(n) for n <= 32. - Reinhard Zumkeller, Oct 30 2011
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
PROG
(PARI) prim(f)=for(i=1, #f~, if(f[i, 1]%3!=1 && (f[i, 1]!=3 || f[i, 2]>1), return(factorback(f)==0))); 1
imprim(f)=my(t); for(i=1, #f~, if(f[i, 1]%3<2 && f[i, 2]>1, t=1); if(f[i, 1]%3==2, if(f[i, 2]%2, return(0), t=1))); t
is(n)=my(f=factor(n)); prim(f)+imprim(f)==1 \\ Charles R Greathouse IV, Nov 04 2015
CROSSREFS
After the initial 0, differs from A329963 next time at a(63) = 196, term which is not present in the latter.
KEYWORD
nonn
AUTHOR
EXTENSIONS
Extended by Ray Chandler, Jan 29 2009
Data section further extended up to a(71), to better differentiate from nearby sequences, by Antti Karttunen, Jul 04 2024
STATUS
approved