%I#9 2021年3月12日22:24:48
%S 1,0,1,1,2,1,3,2,5,3,7,5,11,7,15,11,22,15,30,22,31,56,43,77,58101,
%电话:801351061771422321872992463888319495415634532803683,
%电话:1017869127711031605139020001751249221893087233819
%N f(x^3,x^21)/f(-x^2,-x^4)的展开式,其中f(,)是Ramanujan广义θ函数。
%C Ramanujan theta函数:f(q)(见A121373)、phi。
%H G.C.Greubel,n的表格,n=0..2500时的a(n)</a>
%H Michael Somos,《Ramanujan theta函数简介》</a>
%H Eric Weisstein的数学世界,<a href=“http://mathworld.wolfram.com/RamanujanThetaFunctions.html“>Ramanujan Theta函数</a>
%周期48序列的F Euler变换[0,1,1,1,0,0,0,0,1,0,1,0,1,0,1,0,1,1,1,0,1,1,1,0,0,0,0。
%F a(n)=-A143067(2*n+3)。
%e.G.f.=1+x^2+x^3+2*x^4+x^5+3*x^6+2*x^7+5*x^8+3*x^9+。。。
%e G.f=q^77+q^173+q^221+2*q^269+q^317+3*q^365+2*qq^413+。。。
%t a[n_]:=系列系数[QPochhammer[-x^3,x^24]QPochharmer[-x^21,x^24]QPoch hammer[x^24]/QPochhamer[x^2],{x,0,n}];
%o(PARI){a(n)=如果(n<0,0,a=x*o(x^n));polceoff(subst(prod(k=1,n\3,1-x^k*[1,1,0,0,0,0,1][k%8+1],1+x*o;
%Y参考A143067。
%K nonn公司
%0、5
%A _迈克尔·索莫斯,2015年9月10日
|