%I#22 2022年5月22日09:47:48
%S 1,0,2,1,4,1,5,0,8,1,11,7,8,4,4,4-11,13,1,6,15,13,8,12,22,25,14,9,
%电话:13,11,30,9,16,25,3,12,35,41,29,29,11,27,43,32,16,2,5,29,3,2,
%U 30,18,30,32,56,44,38,44,27,4
%N平方rt(2)=Sum_{N>=1}a(N)/N!的阶乘展开!,使用贪婪算法。
%H G.C.Greubel,n的表格,n=1..10000的a(n)</a>
%H<a href=“https://oeis.org/index/Fa#facbase“>阶乘基表示的索引条目</a>
%平方(2)=1+0/2!+2/3! + 1/4! + 4/5! + 4/6! + 1/7! + 5/8!+。。。
%p A009949:=程序(a,n)局部i,b,c;b:=a;c:=[楼层(b)];对于i从1到n-1,做b:=b-c[i]/i!;c:=[op(c),楼层(b*(i+1)!)];od;c;结束时间:
%t带[{b=Sqrt[2]},表[If[n==1,Floor[b],Floor[n!*b]-n*Floor[(n-1)!*b]],{n,1100}]](*_G.C.Greubel_,2018年12月10日*)
%o(PARI)默认值(realprecision,250);b=平方英尺(2);对于(n=1,80,打印1(如果(n==1,楼层(b),楼层(n!*b)-n*楼层((n-1)*b) ),“,”)\\_G.C.Greubel_,2018年12月10日
%o(PARI)默认值(realprecision,900);我的(t=sqrt(2));对于(n=1,80,t=t*n;打印1(楼层(t),“,”);t=压裂(t);\\_Joerg Arndt_2018年12月17日
%o(Magma)SetDefaultRealField(RealFild(250));[楼层(平方(2))]类别[楼层(阶乘(n)*Sqrt(2)_G.C.Greubel,2018年12月10日
%o(Sage)b=平方(2);
%o定义a(n):
%o如果(n==1):返回楼层(b)
%o else:返回展开(floor(factorial(n)*b)-n*floor(阶乘(n-1)*1))
%o【a(n)代表n in(1..80)】#_G.C.Greubel_,2018年12月10日
%Y参考A002193(十进制扩展),A040000(续分数)。
%Y参见A067881(平方码(3))、A068446(平方码)、A320839(平方码,7)。
%K nonn公司
%O 1,3
%A _N.J.A.Sloane,_比尔·戈斯珀_
|