G[n_]:=模[{G},Do[G[x_]=1+x*(G[x]^3/6+G[x^2]*G[x]/2+G[x^3]/3)+O[x]^n//正规,{n}];克[x]];
T[n_,k_]:=模[{T=G[n],G},T=x*((T^2+(T/.x->x^2))/2);g[e_]=(正规[t+O[x]^商[n,e]]/.x->x^e)+O[x]^n//正规;系数[(和[EulerPhi[d]*g[d]^(k/d),{d,除数[k]}]/k+如果[OddQ[k],g[1]*g[2]^商[k,2],(g[1]^2+g[2])*g[2]^(k/2-1)/2])/2,x,n]];
a[n]:=T[n+5,5];
|