%I#10 2022年9月8日08:45:05
%S 1,1,0,3,2,2,2,3,6,6,8,11,12,7,6,13,7,3,2,2,9,20,9,16,11,0,12,13,
%电话:19,25,26,31,18,24,21,32,12,12,34,22,24,13,14,41,20,34,29,22,50,4,33,
%U 50,39,8,15,24,14,59,40,3,9,29,27,14,18,39,59,44,28,30,35,5,64,20,18
%N zeta(2)的阶乘展开:zeta(1)=Sum_{N>=1}a(N)/N!。
%H G.C.Greubel,n的表格,n=1..10000的a(n)</a>
%F a(n)=楼层(n!*zeta(2))-n*楼层((n-1)*zeta(2)),对于n>=2。
%t与[{b=Zeta[2]},表[If[n==1,Floor[b],Floor[n!*b]-n*Floor[(n-1)!*b],{n,1,100}]](*_G.C.Greubel_,2018年11月26日*)
%o(PARI)默认值(realprecision,250);b=zeta(2);对于(n=1,80,打印1(如果(n==1,楼层(b),楼层(n!*b)-n*楼层((n-1)*b) ),“,”)\\_G.C.Greubel_,2018年11月26日
%o(Magma)SetDefaultRealField(RealFild(250));L: =黎曼齐塔();[楼层(Evaluate(L,2))]cat[楼层(Factorial(n)*Evaluation(L,2中))-n*楼层(Factor(n-1))*Eval(L,二)):n in[2..80]];//_G.C.Greubel,2018年11月26日
%o(鼠尾草)
%o定义A067279(n):
%o如果(n==1):返回楼层(zeta(2))
%o else:返回展开(floor(factorial(n)*zeta(2))-n*floor(阶乘(n-1)*zeta2))
%o【A067279(n)表示n in(1..80)】#_G.C.Greubel_,2018年11月26日
%Y参见A067277(zeta(3))、A068447(zeta(4))、AO68454(zeta(5))、P068455(zeta)(6))、C068456(zeta。
%K容易,不是
%O 1,4型
%2002年3月10日,A _贝尼特·克洛伊特
%E a(1)由G.C.Greubel修正,2018年11月26日
|