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[3], s] /.
Table[s[i] -> Sum[b[[k]] x^(k*i), {k, 1, nn}], {i, 1, 3}], {x, 0,
