OFFSET
1
COMMENTS
LINKS
FORMULA
a(n) <= A377874(n).
EXAMPLE
A276085(174) = 223092873 = 3^4 * 1063 * 2591, which has a divisor of the form p^e, with e >= p, thus a(174) = 0.
PROG
(PARI)
A276085(n) = { my(f = factor(n), pr=1, i=1, s=0); for(k=1, #f~, while(i <= primepi(f[k, 1])-1, pr *= prime(i); i++); s += f[k, 2]*pr); (s); };
A377868(n) = if(isprime(n), 1, my(x=A276085(n), pp); forprime(p=2, , pp = p^p; if(!(x%pp), return(0)); if(pp > x, return(1))));
(PARI)
A359550(n) = { my(pp); forprime(p=2, , pp = p^p; if(!(n%pp), return(0)); if(pp > n, return(1))); };
A377868(n) = if(1==n, 0, if(isprime(n), 1, my(f = factor(n), pr=1, i=1, s=0); if(2==vecsum(f[, 2]), return(2==f[1, 1]), for(k=1, #f~, while(i <= primepi(f[k, 1])-1, pr *= prime(i); i++); s += f[k, 2]*pr); A359550(s)))); \\ Antti Karttunen, Nov 17 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Nov 10 2024
STATUS
approved