OFFSET
4,3
LINKS
Andrew Howroyd, Table of n, a(n) for n = 4..500
FORMULA
G.f.: A(x) = ( B(x)^4 + B(x^2)^2 + 2*B(x^4) )/4 where B(x) is the o.g.f. for A000081. - Geoffrey Critzer, Aug 09 2013
MATHEMATICA
nn = 20; f[x_] := Sum[a[n] x^n, {n, 0, nn}]; sol =
SolveAlways[
0 == Series[
f[x] - x Product[1/(1 - x^i)^a[i], {i, 1, nn}], {x, 0, nn}],
x]; b = Flatten[Table[a[n], {n, 1, nn}] /. sol]; CoefficientList[
Series[CycleIndex[CyclicGroup[4], s] /.
Table[s[i] -> Sum[b[[k]] x^(k*i), {k, 1, nn}], {i, 1, 4}], {x, 0,
nn}], x] (* Geoffrey Critzer, Aug 08 2013 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved