OFFSET
0,6
COMMENTS
a(n) = index of the least non-unitary prime divisor of A276087(n) or 0 if no such prime-divisor exists.
LINKS
PROG
(PARI)
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A328828(n) = { my(i=1, p=2); while(n, if((n%p)>1, return(i)); i++; n = n\p; p = nextprime(1+p)); (0); };
(PARI)
A277885(n) = if(1==n, 0, my(f=factor(n)); for(i=1, #f~, if(f[i, 2]>1, return(primepi(f[i, 1])))); (0));
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 29 2019
STATUS
approved
