%我#21 2021年3月12日22:24:46
%S 1,3,6,11,20,34,56,9114322033449873210641528217130584269,
%电话:591081241108815034202642715436189479886332483176108780,
%电话:141672183776237499305812392406501856639781813108103054130192816405722061850
%N将2*N+1划分为不等于0,+-4,+-6,+-10,16的部分(mod 32)。
%C Ramanujan theta函数:f(q)(见A121373)、phi(q)
%H G.C.Greubel,n表,n=0..1000时的a(n)</a>
%H Michael Somos,《Ramanujan theta函数简介》</a>
%H Eric Weistein的《数学世界》,<a href=“http://mathworld.wolfram.com/RamanujanThetaFunctions.html“>Ramanujan Theta函数</a>
%F(phi(q^2)/phi(-q)-1)/(2*q)的q次幂展开式,其中phi()是Ramanujan theta函数。
%周期16序列的F欧拉变换[3,0,1,2,1,1,2,3,0,3,2,1,2,1,1,0,3,0…]。
%F 2*a(n)=A208850(n+1)。a(n)=A185083(n+1)。
%e 1+3*q+6*q^2+11*q^3+20*q^4+34*q^5+56*q^6+91*q^7+143*q^8+。。。
%e a(2)=6,因为2*2+1=5=3+2=3+1+1=2+2+1=2+1+1=1+1+1+1+1 6种方式。
%t a[n_]:=级数系数[(椭圆θ[3,0,q^2]/椭圆θ[3],0,-q]-1)/(2*q),{q,0,n}];表[a[n],{n,0,50}](*_G.C.Greubel_,2018年3月5日*)
%o(PARI){a(n)=局部(a);如果(n<0,0,n=2*n+2;a=x*o(x^n);polcoeff((eta(x^2+a)^3/(eta
%Y参考A185083、A208850。
%K nonn公司
%0、2
%A _迈克尔·索莫斯,2012年3月2日