n=30;(*来自Rains和Sloane的算法*)
S4[f_,h,x_]:=f[h,x]^4/24+f[h,x]^2 f[h,x^2]/4+f[h,x]f[h,x^3]/3+f[h,x^2]^2/8+f[h,x^4]/4;
S5[f,h,x_]:=f[h,x]^5/120+f[h、x]^3f[h和x^2]/12+f[h,x]*2f[h;x^3]/6+f[h[h,x]f[h;
T[-1,z_]:=1;T[h_,z_]:=T[h,z]=表[z^k,{k,0,n}]。取[系数表[z^(n+1)+1+S4[T,h-1,z]z,z],n+1];
求和[Take[CoefficientList[z^(n+1)+S5[T,h-1,z]z-S5[T,h-2,z]z-(T[h-1,z]-T[h-2,z])(T[h-1,z]-1),z],n+1],{h,1,n/2}]+PadRight[{0,1},n+1)(*罗伯特·拉塞尔2018年9月15日*)