%I#29 2022年9月8日08:44:46
%S 1,-1,-2,-1,-1,5,1,13,4,0,2,-8,-61,-31,13,-156,211223,92,91426,
%电话:972165141、-11264401294、-4684、-2755、-5748、-2414、-667910511,
%U-10048,-1936919635226291402776969,-199040193,-1067879525215767,-54322,-40882
%N产品的扩展{m>=1}(1-m*q^m)。
%C(9)是这个序列中唯一出现的0吗_罗伯特·伊斯雷尔(Robert Israel),2015年6月2日
%H Robert Israel,n的表,a(n)表示n=0..10000</a>
%p p:=倍数(1-m*q^m,m=1..100):
%p S:=系列(p,q,101):
%p seq(系数(S,q,j),j=0..100);#_罗伯特·伊斯雷尔(Robert Israel),2015年6月2日
%p#第二个Maple程序:
%p b:=proc(n,i)选项记忆`如果`(i*(i+1)/2<n,0,
%p`if`(n=0,1,b(n,i-1)-`if`
%p端:
%pa:=n->b(n$2):
%p序列(a(n),n=0..60);#_肖恩·艾尔文(Sean A.Irvine_)(以阿洛伊斯·海因茨(_Alois P.Heinz)命名),2019年5月19日
%t nmax=40;系数列表[系列[产品[1-k*x^k,{k,1,nmax}],{x,0,nmax{],x](*_Vaclav Kotesovec_,2015年12月15日*)
%t nmax=40;系数列表[系列[Exp[-Sum[PolyLog[-j,x^j]/j,{j,1,nmax}],{x,0,nmax{],x](*_Vaclav Kotesovec_,2015年12月15日*)
%t(*更高效的程序:*)nmax=50;poly=常量数组[0,nmax+1];聚[1]]=1;聚[2]]=-1;Do[Do[poly[[j+1]]-=k*poly[[j-k+1]],{j,nmax,k,-1}];,{k,2,nmax}];poly(*Vaclav Kotesovec_,2016年1月7日*)
%o(PARI)m=50;q='q+O('q^m);Vec(产品(n=1,m,(1-n*q^n))\\_G.C.Greubel_,2018年2月18日
%o(Magma)系数(&*[(1-m*x^m):[1.40]]中的m)[1.40],其中x是多项式环(整数()).1;//_G.C.Greubel,2018年2月18日
%Y参见A006906、A022629、A022693、A266964。
%K符号
%0、3
%A _N.J.A.斯隆_
|