OFFSET
1,2
COMMENTS
1 is the only odd term in this sequence.
LINKS
Amiram Eldar, Table of n, a(n) for n = 1..10000
MATHEMATICA
f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[AllTrue[e, OddQ], Times @@ (e + 1), Nothing]]; f[1] = 1; Array[f, 150]
PROG
(PARI) lista(max) = for(k = 1, max, my(e = factor(k)[, 2], isexpodd = 1); for(i = 1, #e, if(!(e[i] % 2), isexpodd = 0; break)); if(isexpodd, print1(vecprod(apply(x -> x+1, e)), ", ")));
(Python)
from math import prod
from itertools import count, islice
from sympy import factorint
def A366438_gen(): # generator of terms
for n in count(1):
f = factorint(n).values()
if all(e&1 for e in f):
yield prod(e+1 for e in f)
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Amiram Eldar, Oct 10 2023
STATUS
approved