%I#19 2017年5月26日22:18:24

%S 1,1,4,2,27,2248264,3028303610830,83823255768187461078，

%电话：59385996174041472046986576104017121832697694577881593664，

%电话：6268211220260189838270667422652636262499413103569882452962964342263000869747536857008706154096514160998750490967416896101551822350213865372020154257135428048205014786580878437729442000322560

%N行读取的三角形T（N，k）：公式部分中定义的多项式P_N（T）的系数。

%C行n>0包含楼层（2*（n+1）/3）术语。

%H Gheorghe Coserea，<a href=“/A286798/b286798.txt”>行n=0..123，扁平</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 y（x；t）=Sum_{n>=0}P_n（t）*x^n满足x^2*导数（y，x）=（1-y+x*y^2+2*x^2*t*y^3）/（t-（2+t）*y-3*x*t*y ^2），y（0；t）=1，其中P_n。

%F A000699（n+1）=T（n，0），A000108（n）=P_n（-1），A286799（n）=P_n。

%e A（x；t）=1+x+（4+2*t）*x^2+（27+22*t）*x^3+（248+264*t+30*t^2）*x*4+

%e三角形开始：

%电子[0][1][2][3][4][5]

%e[0]1；

%e[1]1；

%e[2]4，2；

%e[3]第27、22页；

%e[4]第248、264、30页；

%电子[5]2830，3610，830，8；

%e[6]38232、55768、18746、1078；

%电子[7]593859、961740、414720、46986、576；

%电子[8]10401712、18326976、9457788、1593664、62682、112；

%电子[9]202601898，382706674，226526362，49941310，3569882，45296；

%e[10]。。。

%t最大值=12；y0[x_，t]=1；y1[x_，t]=0；对于[n=1，n<=max，n++，y1[x_，t_]=1+x y0[x，t]^2+3 tx^3 y0[x，t]^2 D[y0[x，t]，x]+x^2（2 y0[x，t]D[y0[x，t]，x]+t（2 y0[x，t]^3-D[y0[x，t]，x]+y0[x，t]D[y0[x，t]，x]）+O[x]^n//正常//简化；y0[x_，t]=y1[x，t]]；

%t P[n_，t_]：=系数[y0[x，t]，x，n]；

%t行[n_]：=系数列表[P[n，t]，t]；

%t表[第[n]行，{n，0，max}]//Flatten（*_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；）；y0；

%o}；

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

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

%o}；

%o concat（应用（p->Vecrev（p），Vec（A286798_ser（12）））

%o\\测试：y=A286798_ser（50）；x^2*y'==（1-y+x*y^2+2*x^2*t*y^3）/（t-（2+t）*y-3*x*t*y*2）

%Y参见A286781、A286782、A286793、A28678、A28675。

%K nonn，标签

%0、3

%A _Gheorghe Coserea，2017年5月21日