L2=N[Log[10,2],50];L5=N[Log[10,5],50];k=1;Do[While[Take[RealDigits[10^ FractionalPart[L2*k]][[1],n]!=取[RealDigits[10^分数部分[L5*k]][[1],n],k++];打印[k],{n,1,10}]
L2=N[Log[10,2],50];L5=N[Log[10,5],50];k=1;Do[While[Take[RealDigits[10^ FractionalPart[L2*k]][[1],n]!=取[RealDigits[10^分数部分[L5*k]][[1],n],k++];打印[k],{n,1,7}]
f[n_,k_]:={楼层[10^(k-1+n[分数部分[n Log[5]/Log[10]],20])]整平@块[{$MaxExtraPrecision=\[Infinity]},块[{l=Denominator/@Convergents[Log10[2],1000]},数组[k\[Function]l[[压扁@位置[f[#/2,k]&/@l,{x_,x_},{1},1]]]/2,20]]](*J.Mulder(jasper.Mulder(AT)plant.nl),2010年2月3日*)
(*备用程序*)
n=100$MaxExtraPrecision=n;安=
连续分数[Log10[5/2],n];数据=
分母/@
压扁[桌子[
FromContinuedFraction[Join[ans[[1;;p-1]],{#}]&/@
范围[1,ans[[p]]],{p,2,n}]];溶胶=
选择[表[{k,a=N[分数部分[{k Log10[2],k Log10[5]}],N];
10^a,b=实际数字[10^a][[All,1]];
LengthWhile[Range[Length[b[[1]]],b[[1,#]]==b[[2,#]]&],
10年。{-1,1},实际数字[10^a.{-1,1}][[-1]]},{k,数据}],
绝对值[#[[-2]]<1&];
acc=关联[{}];s=溶胶[[全部,{1,3}]];对于[i=1,
i<长度[s],i++,
如果[Lookup[acc,s[[i,2]],0]==0,
acc[s[[i,2]]]=s[[i,1]]];最终=
Rest[Sort[Normal[acc]]/。规则->列表;
bcc=关联[{}];对于[i=Max[Keys[acc]],i>=Min[Keys[Cac]],i--,
j=i;而[Lookup[acc,j,0]==0&&j<Max[Keys[acc]],j++];
bcc[i]=符合[j];j=i;而[bcc[j]>=bcc[j+1],j++];
bcc[i]=最小值[bcc[i],bcc[j]]];英国广播公司=
Rest[Sort[Normal[Reverse[bcc]]]/。规则->列表(*王贤文,2023年6月2日*)
|