%I#9 2019年5月30日22:54:11

%S 1,1,2,3,6,14,31,7417642510552637667111699543507112102，

%电话：29034775583619767175189664136733143613954395795807254610593，

%电话：6783859181162593148481779101299991452334921821516939712649412532323626836402935184793512776850017033925611355476815244236777239029896764310264271633640605450739353335433772014350935223715549303915708822614992221809540249

%N G.f.：A（x）=exp（和{N>=1}x ^N/N*（1+x ^N*A（x）^（2*N））/（1+x^N*A（x）^N））。

%H Paul D.Hanna，n表，n=0..300时的a（n）</a>

%电子表格：A（x）=1+x+x^2+2*x^3+3*x^4+6*x^5+14*x^6+31*x^7+74*x^8+176*x^9+425*x^10+1055*x^11+2637*x^12+6671*x*^13+16995*x^14+43507*x^15+。。。

%e这样

%e log（A（x））=x*（1+x*A（x）^2）/（1+x*A（x（1+x^5*A（x）^10）/。。。

%e明确表示，

%e对数（A（x））=x+x^2/2+4*x^3/3+5*x^4/4+16*x^5/5+46*x*x^6/6+113*x^7/7+317*x^8/8+823*x^9/9+2206*x^10/10+6051*x^11+16418*x^12/12+45007*x^13/13+123033*x^14/14+336244*x^15+15+。。。

%o（PARI）{a（n）=my（a=[1]）；对于（i=1，n，a=Vec（exp（总和（m=1，#a+1，x^m/m*（1+x^m*Ser（a）^（2*m））/（1+x^m*Ser（a，^m））））；a[n+1]}

%o表示（n=0,40，打印1（a（n），“，”）

%K nonn公司

%0、4

%A _保罗·D·汉纳，2019年3月29日