
MATHEMATICA

n = 23; p[s_, 1] := c[s]; p[0, n_] := c[0]^n; p[s_, n_] := p[s, n] = Expand[Sum[c[k] p[s  k, n  1], {k, 0, s}]]; d[n_, 0] := Sum[c[k] c[0]^k, {k, 0, n}]; d[n_, i_] := Sum[c[k] p[i, k], {k, 1, n}]; a = c[0] /. FindRoot[Table[d[n, k] k! == 1, {k, 0, n}], Table[{c[k], 2^k}, {k, 0, n}], WorkingPrecision > 30]; First[RealDigits[a, 10, 6]]
