%I#15 2020年8月11日19:40:47
%S 1,2,11,14,25,28,30,37,39,41,43,46,48,51,54,57,60,64,66,71,76,78,80,
%电话82,84,90,95101103106110113115117127133140146152157,
%电话:160162165167170173179181185189196200203206209212215220
%(8/13)的贪婪力量:和{N>=1}(8/13,^a(N)=1。
%当0.5<x<1时,x的第n个贪婪幂是不导致幂级数sum_{k=1..n}x^a(k)超过单位的最小整数指数a(n)。
%C一个启发式论证表明,A(n)/n的极限是m-sum_{n=m..inf}log(1+x^n)/log(x)=3.8170308430…,其中x=8/13,m=floor(log(1-x)/log(x))=1_Paul D.Hanna,2002年11月16日
%当你达到sum_{n=1..59}(8/13)^a(n)时,这个和和1之间的差只有1.6*10^-47。
%H Robert Israel,<a href=“/A0777475/b0777475.txt”>n,a(n)表,n=1.-10000</a>
%F a(n)=sum_{k=1..n}楼层(g_k),其中g_1=1,g_{n+1}=log_x(x^frac(g_n)-x)(n>0),x=(8/13),frac(y)=y-楼层(y)。
%F a(n)似乎渐近于c*n,c约为3.7…-贝诺伊特-克洛伊特_
%e a(3)=11,因为(8/13)+。
%p V:=矢量(100):
%p V[1]:=1:温度:=1-8/13:
%p代表从2到100 do的n
%p V[n]:=-楼层(log[13/8](T));
%p T:=T-(8/13)^V[n];
%日期:
%p转换(V,列表);#_罗伯特·伊斯雷尔(Robert Israel),2020年8月11日
%t s=0;a={};做[如果[s+(8/13)^n<1,s=s+(8:13)^n;a=Append[a,n]],{n,1,250}];一
%t启发式极限[x_]:=(m=楼层[Log[x,1-x]])+1/24+Log[x、乘积[1+x^n,{n,1,m-1}]/DedekindEta[I Log[x]/-Pi]*DedekindEta[-I Log[x]/2/Pi]];N[启发式极限[8/13],20]
%Y参考A077468、A077469、A077470、A077471、A077472、A077473、A077474。
%K容易,不是
%O 1,2号机组
%A·保罗·D·汉纳,2002年11月6日
%E由_Robert G.Wilson v_编辑和扩展,2002年11月8日。2002年11月6日由_Benoit Cloitre_延长
|