%I#17 2022年3月8日14:17:24

%序号1,2,10,82898120181876263323682656076821424967394337908874，

%电话：8643725766262382554347123470307467263201822118247888976170，

%电话：739081808704195650261461161294004838429763820587771744516503838629686672749972852215849866941237056394386

%N行A286796的总和。

%H Gheorghe Coserea，n的表格，n=0..200的a（n）</a>

%H Luca G.Molinari，Nicola Manini，<a href=“https://arxiv.org/abs/cond-mat/0512342“>多体骨架图的枚举，arXiv:cond-mat/0512342[cond-mat.str-el]，2006。

%F a（n）=和{k=0..n}A286796（n，k）。

%F a（n）~2^（n+5/2）*n ^（n+2）/exp（n+2）。-_Vaclav Kotesovec_，2022年3月8日

%t最大值=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]]；

%t a[n_]：=系数列表[系列系数[y0[x，t]/（1-x*t*y0[x，t]），{x，0，n}]，t]//总计；

%t表[a[n]，{n，0，max-1}]（*_Jean-François Alcover_，2017年5月24日，改编自PARI*）

%o（PARI）

%o A286795_ser（N，t=t）={

%o我的（x='x+o（'x^N），y0=1，y1=0，N=1）；

%o while（n++，

%o y1=（1+x*（1+2*t+x*t^2）*y0^2+t*（1-t）*x^2*y0,3+2*x^2*y0*y0'）；

%o y1=y1/（1+2*x*t）；如果（y1==y0，break（））；y0=y1；)；年；

%o}；

%o A286796_ser（N，t=t）=我的（v=A286795_ser（N，t））；v/（1-x*t*v）；

%o Vec（A286796_ser（20,1））

%o（PARI）

%o A049464_ser（N）={\\对于A049464（0）=0

%o my（s=Ser（concat（1，向量（N+1，N，（2*N）！）！/（2^n*n！）），g=（1/s-1/s^2）/x）；

%o 1-1/subst（g，'x，serreverse（x*g^2*s^2））；

%o}；

%o A286797_ser（N）=我的（q=A049464_ser（N））；q/（x-x*q）；

%o Vec（A286797_ser（20））

%Y参考A286796。

%K nonn公司

%O 0,2

%A Gheorghe Coserea，2017年5月21日