OFFSET
1,4
LINKS
FORMULA
Multiplicative with a(p^e) = A000078(e+3).
MATHEMATICA
t[e_] := t[e] = If[e < 5, 2^(e-1), t[e-1] + t[e-2] + t[e-3] + t[e-4]]; a[1] = 1; a[n_] := Times @@ t /@ Last @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Oct 23 2023 *)
PROG
(PARI) for(n=1, 87, print1(direuler(p=2, n, 1/(1 - X - X^2 - X^3 - X^4))[n], ", "))
CROSSREFS
KEYWORD
nonn,mult
AUTHOR
Ilya Gutkovskiy, Feb 16 2022
STATUS
approved