%I#21 2022年9月8日08:45:05
%S 1,0,0,0,1,0,5,0,9,1,3,1,0,1,11,2,15,18,4,16,9,20,12,20,
%电话:14,22,0,8,9,3,26,15,6,13,11,20,32,7,12,31,39,46,6,49,7,10,2,5,44,
%U 11,32,41,49,21,40,17,49,62,44,13,25,67,41,57,27,13,24,35,25,43,25,27,29
%N zeta(10)的阶乘展开:zeta(100)=Sum_{N>0}a(N)/N!。
%H G.C.Greubel,n表,n=1..100000的a(n)(文森佐·利班迪的术语1..150)
%t t=泽塔[10];s={};Do[n=楼层[t*i!];t-=不适用!;附加到[s,n],{i,1,30}];s(*_Amiram Eldar_,2018年11月25日*)
%t与[{b=Zeta[10]},表[If[n==1,Floor[b],Floor[n!*b]-n*Floor[(n-1)!*b],{n,1,100}]](*_G.C.Greubel_,2018年11月26日*)
%o(PARI)向量(30,n,if(n>1,t=t%1*n,t=zeta(10))\1)\\增加实精度(例如\p500)以计算更多项。-_M.F.Hasler,2018年11月25日
%o(PARI)默认值(realprecision,500);对于(n=1,80,print1(如果(n==1,floor(zeta(10)),floor*zeta(10)),“,”)\\_G.C.Greubel_,2018年11月26日
%o(Magma)SetDefaultRealField(RealFild(250));五十: =黎曼泽塔(RiemannZeta);[楼层(Evaluate(L,10))]cat[楼层(Factorial(n)*Evaluation(L,0))-n*楼层(Factor((n-1))*Eval(L,100)):n in[2..80]];//_G.C.Greubel,2018年11月26日
%o(鼠尾草)
%o定义A068459(n):
%o如果(n==1):返回楼层(zeta(10))
%o else:返回展开(floor(factorial(n)*zeta(10))-n*floor(阶乘(n-1)*zeta10))
%o【A068459(n)代表n in(1..80)】#_G.C.Greubel_,2018年11月26日
%Y参考A007514。
%K nonn公司
%O 1,7
%2002年3月10日,A _贝尼特·克洛伊特
%E Keywords cons,由M.F.Hasler轻松删除,2018年11月25日
|