nn=10; Flatten[{1, Table[Subscript[c, j]*j!, {j, 1, nn}]/.Solve[Table[SeriesCoefficient[x/(x*(1+Sum[Subscript[c, j]*x^j, {j, 1, nn}])-1), {x, 0, k}]==SeriesCoefficient[Tan[-Sum[Subscript[c, j]*x^j, {j, 1, nn}]], {x, 0, k}], {k, 0, nn}]]}] (*
Vaclav Kotesovec, Jun 15 2013 *)