OFFSET
1,1
COMMENTS
LINKS
PROG
(PARI)
A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
isA359830(n) = !!(A048675(n)%3);
(Python)
from itertools import count, islice
from sympy import primepi, factorint
def A359830_gen(startvalue=1): # generator of terms >= startvalue
return filter(lambda n:sum(e<<primepi(p)-1 for p, e in factorint(n).items())%3, count(max(startvalue, 1)))
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Jan 17 2023
STATUS
approved