%I#25 2025年6月12日15:52:57

%S 1,1,2,1,1,2,3,4,3,3,,3,6,7,8,10,9,9,9，9,11,13,16,20,22,25,28,27,28，

%电话：29,30,32,35,40,45,53,60,67,73,79,85,87,92,95,98105111120132145，

%电话：160178196212231247263280291305319334352371393

%N计算一组受限分区。

%H Vaclav Kotesovec，n=0..10000的n，a（n）表</a>

%F G.F.：求和{k>=0}x^（k^2）*Product_{j=1..k}（1+x^j）^2=1+x^1*（1+x）^2+x^4*（1++x）^2*（1+x^2）^2+。..+x^k^2*（1+x）^2*。..*（1+x^k）^2+。..

%F a（n）~phi^（3/2）*exp（Pi*sqrt（2*n/15））/（4*5^（1/4）*sqrt（n）），其中phi=A001622=（1+sqrt）/2是黄金比率。-_Vaclav Kotesovec_，2024年9月29日

%t取[系数表[Sum[x^（k^2）*乘积[1+x^i，{i，k}]^2，{k，0，7}]，x]，63]（*_Giovanni Resta_，2013年3月13日*）

%t nmax=100；p=1；s=1；做[p=展开[p*（1+x^k）*（1+x^k，*x^（2*k-1）]；p=取[p，Min[nmax+1，指数[p，x]+1，长度[p]]]；s+=p；，{k，1，平方[nmax]}]；取[系数表[s，x]，nmax+1]（*_Vaclav Kotesovec_，2024年10月9日*）

%Y参见A001622、A306734、A376542、A376580、A376812和A376813。

%K nonn公司

%0、3

%A _David S.Newman，2013年3月13日

%E a（14）-a（62）摘自Giovanni Resta，2013年3月13日