%I#29 2022年9月8日08:46:22
%S 2,5,3,4,3,8,3,2,4,7,2,5,9,4,9,6,7,2,4,21,5,7,3,5,3,11,
%电话:4,5,7,8,2,6,23,2,5,6,14,3,16,3,2,3,14,2,4,15,17,7,4,2,11,18,2,4,47,
%U 14,5,6,4,2,7,3,2,3,13,3,5,15,4,8,6,9,2,4,11,6,5,22,4
%N a(N)是最小正整数k,因此k^N以2开头。
%C对于n>1,取整数d>-n log_10(3^(1/n)-2^(1/1n))。
%C然后10^(d/n)>1/(3^(1/n)-2^(1/1n))
%所以(3*10^d)^
%C,因此a(n)<=上限(2*10^d)^(1/n))。
%C特别是,a(n)总是存在的。
%H Robert Israel,n的表,a(n)表示n=1..10000</a>
%F a(n)=A067443(n)^(1/n)。
%F A000030(a(n)^n)=2。
%e a(5)=3,因为3^5=243以2开头,而1^5=1和2^5=32不以2开头。
%p f:=进程(n)局部x,y;
%从2do到x的p
%p y:=x^n;
%p如果地板(y/10^ilog10(y))=2,则返回x fi
%日期
%p端程序:
%p映射(f,[1..100]美元);
%o(PARI)a(n)=表示(k=1,oo,if(数字(k^n)[1]==2,返回(k)))
%o(Python)
%o n=1
%o当n<100时:
%o k,s=2,str(2**n)
%o whiles s[0]!="2":
%o k=k+1
%o s=str(k**n)
%o打印(n,k)
%o n=n+1#_A.H.M.Smeets_,2019年8月14日
%o(岩浆)m:=1;溶胶:=[];对于[1..100]中的n,do k:=2;反向(Intseq(k^n))[1]ne 2 do k:=k+1;结束while;溶胶[m]:=k;m: =m+1;结束;溶胶;//_Marius A.Burtea,2019年8月15日
%Y参考A000030、A067443、A309707。
%K nonn,基础
%O 1,1号机组
%2019年8月14日,以色列罗贝尔特
|