%I#15 2023年3月31日05:13:21
%S 1,0,0,1,1,0,2,1,1,1,3,2,2,5,6,3,5,8,9,7,8,13,14,10,14,19,20,17,
%电话:20,29,30,26,32,42,45,41,47,63,64,60,70,88,91,87,99124128123143,
%电话:172179176200240246246279325338381440456461519590615
%N将N划分为不含奇数平方和偶数(非正方形)的部分的次数,奇数部分以单个重数出现,偶数部分没有限制。
%H G.C.Greubel,n表,n=0..5000的a(n)</a>
%H努里丁主席,<a href=“http://arxiv.org/abs/hep-th/0409011“>部分超对称的分区恒等式</a>,arXiv:hep-th/0409011,2004。
%F G.F.:产品{k>0}(1+x^(2*k-1))/(1-(-1)^k*x ^(k^2))。
%e a(16)=6,因为16=13+3=11+5=7+5+4=5+3+4=4+4+4。
%p系列(乘积((1+x^(2*k-))/(1-(-1)^k*x^)(k^2)),k=1..100),x=0100);
%t带[{m=80},系数列表[系列[产品[(1+x^(2*k-1))/(1-(-1)^k*x^(k^2)),{k,m+2}],{x,0,m}],x]](*_G.C.Greubel_,2023年3月28日*)
%o(岩浆)
%o m:=80;
%o f:=函数<x|(&*[(1+x^(2*k-1))/(1-(-1)^k*x^;
%o R<x>:=PowerSeriesRing(整数(),m);
%o系数(R!(f(x));//_G.C.Greubel,2023年3月28日
%o(SageMath)
%o m=80
%o定义f(x):返回(1,m+2)范围内k的乘积((1+x^(2*k-1))/(1-(-1)^k*x^
%o定义A100522_列表(前c):
%o P.<x>=PowerSeriesRing(ZZ,prec)
%o返回P(f(x)).list()
%o A100522_列表(m)#_G.C.格鲁贝尔,2023年3月28日
%Y参考A100926。
%K nonn公司
%O 0.8
%2004年11月25日,法庭主席
%E偏移由G.C.Greubel修正,2023年3月28日
|