Clear[t, n, k, i, p]
p = 2
coeff = {p^0, -p^1, p^1, p^2, p^3, p^4, p^5, p^6, p^7, p^8, p^9, p^10,
p^11, p^12, p^13, p^14, p^15, p^16};
nn = Length[coeff]
cc = Range[nn]*0 + 1;
Monitor[Do[Clear[t];
t[n_, 1] := t[n, 1] = cc[[n]];
t[n_, k_] :=
t[n, k] =
If[n >= k,
Sum[t[n - i, k - 1], {i, 1, k - 1}] -
Sum[t[n - i, k], {i, 1, k - 1}], 0];
A4 = Table[Table[t[n, k], {k, 1, nn}], {n, 1, nn}];
A5 = A4[[1 ;; nn - 1]];
A5 = Prepend[A5, ConstantArray[0, nn]];
cc = Total[
Table[coeff[[n]]*MatrixPower[A5, n - 1][[All, 1]], {n, 1,