OFFSET
1,2
MATHEMATICA
g[p_, e_] := 5*e; a[1] = 1; a[n_] := Times @@ g @@@ FactorInteger[n]; Array[a, 100]
PROG
(PARI) for(n=1, 100, print1(direuler(p=2, n, (1+3*X+X^2)/(1-X)^2)[n], ", "))
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Vaclav Kotesovec, Feb 28 2023
STATUS
approved