OFFSET
1
LINKS
FORMULA
a(n) = [n==1 or A373147(n)==0], where [ ] is the Iverson bracket.
PROG
(PARI)
A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*prod(i=1, primepi(f[k, 1]-1), prime(i))); };
(PARI) A373489(n) = if(n<=2, !!n, if(isprime(n), 0, my(f=factor(n), ad=n*sum(i=1, #f~, f[i, 2]/f[i, 1]), pr=1, q=1, s=0); for(k=1, #f~, forprime(p = nextprime(1+q), f[k, 1]-1, pr *= p; if(pr > ad, return(0)); q = p); s += f[k, 2]*pr; if(s > ad, return(0))); !(ad%s))); \\ Antti Karttunen, Jan 29 2026
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jun 09 2024
STATUS
approved
