login
A380477
a(n) = 1 if n is not a prime, and A380459(n) has no divisors of form p^p, otherwise 0.
3
1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0
OFFSET
1
FORMULA
a(n) = A380467(n) - A010051(n) = A359550(A380459(n)) - A010051(n).
PROG
(PARI)
A359550(n) = { my(pp); forprime(p=2, , pp = p^p; if(!(n%pp), return(0)); if(pp > n, return(1))); };
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
A349394(n) = { my(p=0, e); if((e=isprimepower(n, &p)), p^(e-1), 0); };
A380459(n) = { my(m=1); fordiv(n, d, m *= A276086(d)^A349394(n/d)); (m); };
A380477(n) = (!isprime(n) && A359550(A380459(n)));
CROSSREFS
Characteristic function of A380478.
Sequence in context: A132918 A358750 A205809 * A353370 A355940 A348033
KEYWORD
nonn
AUTHOR
Antti Karttunen, Feb 03 2025
STATUS
approved