|
|
A318966
|
|
Expansion of e.g.f. Product_{i>=1, j>=1, k>=1} 1/(1 - x^(i*j*k))^(1/(i*j*k)).
|
|
1
|
|
|
1, 1, 5, 21, 165, 1077, 11457, 103905, 1345257, 15834825, 237535389, 3372509709, 59235634125, 979573962429, 19224990899865, 366788042231193, 8019002662543953, 171360055378885905, 4132946756763614133, 97947895990285022085, 2576516749059849502581, 67124117357620005459141
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
E.g.f.: Product_{k>=1} 1/(1 - x^k)^(tau_3(k)/k), where tau_3 = A007425.
E.g.f.: exp(Sum_{k>=1} ( Sum_{d|k} Sum_{j|d} tau(j) ) * x^k/k), where tau = number of divisors (A000005).
|
|
MAPLE
|
a:=series(mul(mul(mul(1/(1-x^(i*j*k))^(1/(i*j*k)), k=1..21), j=1..50), i=1..50), x=0, 22): seq(n!*coeff(a, x, n), n=0..21); # Paolo P. Lava, Apr 02 2019
|
|
MATHEMATICA
|
nmax = 21; CoefficientList[Series[Product[Product[Product[1/(1 - x^(i j k))^(1/(i j k)), {i, 1, nmax}], {j, 1, nmax}], {k, 1, nmax} ], {x, 0, nmax}], x] Range[0, nmax]!
nmax = 21; CoefficientList[Series[Product[1/(1 - x^k)^(Sum[DivisorSigma[0, d], {d, Divisors[k]}]/k), {k, 1, nmax}], {x, 0, nmax}], x] Range[0, nmax]!
a[n_] := a[n] = (n - 1)! Sum[Sum[Sum[DivisorSigma[0, j], {j, Divisors[d]}], {d, Divisors[k]}] a[n - k]/(n - k)!, {k, 1, n}]; a[0] = 1; Table[a[n], {n, 0, 21}]
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|