%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日