terms = 22; A[_] = 1; Do[A[x_] = Sum[j! x^j/Product[(1 - k x A[x]), {k, 1, j}], {j, 0, i}] + O[x]^i, {i, 1, terms}]; CoefficientList[A[x], x]
terms = 22; A[_] = 1; Do[A[x_] = Sum[x^j Sum[k! StirlingS2[j, k] A[x]^(j - k), {k, 0, j}], {j, 0, i}] + O[x]^i, {i, 1, terms}]; CoefficientList[A[x], x]