OFFSET
1,2
FORMULA
a(n) = Sum_{k=1..n} A127136(n,k)*k.
MATHEMATICA
Needs["Combinatorica`"];
nn=30; s[n_, k_]:=s[n, k]=a[n+1-k]+If[n<2 k, 0, s[n-k, k]]; a[1]=1; a[n_]:=a[n]=Sum[a[i] s[n-1, i] i, {i, 1, n-1}]/(n-1); rt=Table[a[i], {i, 1, nn}]; c=Drop[Apply[Plus, Table[Take[CoefficientList[CycleIndex[CyclicGroup[n], s]/.Table[s[j]->Table[Sum[rt[[i]] x^(k*i), {i, 1, nn}], {k, 1, nn}][[j]], {j, 1, nn}], x], nn], {n, 1, 30}]], 1]; Drop[CoefficientList[Series[D[Product[1/(1-y x^i)^c[[i]], {i, 1, nn-1}], y]/.y->1, {x, 0, nn}], x], 1] (* after code given by Robert A. Russell in A000081 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Geoffrey Critzer, Oct 13 2012
STATUS
approved