编辑
经核准的
a(n)~2^(n+5/2)*n^(n+2)/exp(n+2)-瓦茨拉夫·科特索维奇2022年3月8日
提出
最大值=20;y0[x_,t]=1;y1[x_,t]=0;对于[n=1,n<=max,n++,y1[x_,t_]=(1+x*(1+2*t+x*t^2)*y0[x,t]^2+t*(1-t)*x^2*y0[x,t]^3+2*x^2*y0[x,t]*D[y0[x,t],x])/(1+2*x*t)+O[x]^n//正常//简化;y0[x_,t]=y1[x,t]];
a[n_]:=系数列表[系列系数[y0[x,t]/(1-x*t*y0[x,t]),{x,0,n}],t]//总计;
表[a[n],{n,0,max-1}](*Jean-François Alcover公司2017年5月24日,改编自PARI*)
A286799型1986年7月_ser(N)=我的(q=A049464号_ser(N));q/(x-x*q);
Vec公司(A286799型1986年7月_ser(20))
(PARI)
A049464号_ser(N)={\\用于A049464号(0)=0
my(s=Ser(concat(1,向量(N+1,N,(2*N)/(2^n*n!)),g=(1/s-1/s^2)/x);
1-1/subst(g,'x,serreverse(x*g^2*s^2));
};
A286799型_ser(N)=我的(q=A049464号_ser(N));q/(x-x*q);
Vec公司(A286799型_ser(20))