OFFSET
0,3
MATHEMATICA
a[n_] := Coefficient[Series[Product[(x^k + 1 + 1/x^k)^k, {k, 1, n}], {x, 0, 0}], x, 0]; Array[a, 14, 0] (* Amiram Eldar, Dec 24 2021 *)
PROG
(PARI) a(n) = polcoef(prod(k=1, n, (x^k+1+1/x^k)^k), 0);
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Dec 23 2021
STATUS
approved