%I#10 2017年5月26日22:18:32

%S 1,1,6,49542727811382420178813984293486539142220486717908，

%电话5248163121061446387659442675958040851613423937511765492，

%电话：4485155272443939687315865571912065180326592432249691301719190232900865260992371261809616145742392898851628641607005553456012247259844

%N行总和为A286798。

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

%t最大值=22；y0[x_，t]=1；y1[x_，t]=0；对于[n=1，n<=最大值，n++，y1[x_，t_]=1+x y0[x，t]^2+3 t x ^3 y0[x，t]*2 D[y0[x]，t]，x]+x^2（2 y0[x，t]D[y0[x，t]，x]+t简化；y0[x_，t]=y1[x，t]]；

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

%t a[n_]：=系数列表[P[n，t]，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，而（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 A286798_ser（N，t=t）={

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

%o}；

%o Vec（A286798_ser（21,1））

%Y参考A286798。

%K nonn公司

%0、3

%A Gheorghe Coserea，2017年5月21日