|
MATHEMATICA
|
maxExponent = 50; a[0] = 0; a[1] = 1;
a[n_] := Module[{}, aparts = List @@ (Product[1/(1 - x^j), {j, 1, n}] // Apart); cc = aparts + O[x]^maxExponent // CoefficientList[#, x]&; f[k_] = Total[FindSequenceFunction[#, k]& /@ cc]; f[6^n-n(n+1)/2+1] // Round];
Table[an = a[n]; Print[n, " ", an]; an, {n, 0, 10}] (* Jean-François Alcover, Nov 15 2018 *)
|