OFFSET
0,2
LINKS
G. C. Greubel, Table of n, a(n) for n = 0..1000
MAPLE
b:= proc(n) option remember; `if`(n=0, 1, add(
26*(-d)^(n/d+1), d=numtheory[divisors](n)))
end:
a:= proc(n) option remember; `if`(n=0, 1,
add(b(j)*a(n-j), j=1..n)/n)
end:
seq(a(n), n=0..30); # Alois P. Heinz, Jul 18 2018
MATHEMATICA
With[{nmax=50}, CoefficientList[Series[Product[(1+m*q^m)^26, {m, 1, nmax}], {q, 0, nmax}], q]] (* G. C. Greubel, Jul 18 2018 *)
PROG
(PARI) m=50; q='q+O('q^m); Vec(prod(n=1, m, (1+n*q^n)^26)) \\ G. C. Greubel, Jul 18 2018
(Magma) Coefficients(&*[(1+m*x^m)^26:m in [1..40]])[1..50] where x is PolynomialRing(Integers()).1; // G. C. Greubel, Jul 18 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved