OFFSET
1,2
MATHEMATICA
NestList[# + Count[Divisors[#], _?OddQ] &, 1, 62] (* Michael De Vlieger, Oct 25 2021 *)
PROG
(PARI) f(n) = sumdiv(n, d, d%2); \\ A001227
a(n) = if (n==1, 1, my(x=a(n-1)); x + f(x)); \\ Michel Marcus, Oct 26 2021
(Python)
from math import prod
from itertools import islice
from sympy import factorint
def A348623gen(): # generator of terms
n = 1
yield n
while True:
n += prod(q+1 for p, q in factorint(n).items() if p > 2)
yield n
A348623_list = list(islice(A348623gen(), 20)) # Chai Wah Wu, Dec 13 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
Patryk Kisieniowski, Oct 25 2021
STATUS
approved