登录
的行总和A286796型.
6

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