%I#51 2022年7月20日15:19:47
%S 1,4,5,7,10,17,24,37,45,55,70,77107137150271364115716562004,
%电话:21263033348936454336659772791269013840201082169328888,
%电话:84155102930
%N个数m,使得2^m反转为素数。
%C a(35)>105000.-_乔瓦尼·雷斯塔,2013年2月22日
%C From _Bernard Schott_,2022年1月30日:(开始)
%如果m是偶数项,那么u=m/2是A350441的项,这是因为2^m=4^(m/2)。事实上,A350441的项是这个序列中偶数项的一半。
%如果m是3的倍数,那么k=m/3是A350442的倍数,这是因为2^m=8^(m/3)。第一个例子:m=24,45,150,1656。。。以及相应的k=8,15,50552。。。(结束)
%e4是一个项,因为2^4反转为61和素数。
%p与(numtheory):myarray:=[]:对于从1到4000的n,执行it1:=转换(2^n,base,10):it2:=和(10^(nops(it1)-i)*it1[i],i=1..nos(it1)):如果isprime(it2),则打印f(`%d,`,n)fi:od:
%t Do[If[PrimeQ[FromDigits[Reverse[IntegerDigits[2^n]]],打印[n]],{n,20000}](*_Robert G.Wilson v_,2005年1月29日*)
%t选择[Range[4400],PrimeQ[IntegerReverse[2^#]]&](*需要Mathematica版本10或更高版本*)(*程序生成序列的前25个项;要生成更多项,请增加Range常量,但程序将需要更长的运行时间。*)(*Harvey P.Dale_,2020年8月5日*)
%o(Python)
%o来自症状输入isprime
%o k,m,A057708_列表=1,2,[]
%o而k<=10**3:
%o如果isprime(int(str(m)[::-1])):
%o A057708_list.append(k)
%o k+=1
%o m*=2#_恰瓦乌,2021年3月9日
%o(PARI)isok(m)=i素数(来自数字(Vecrev(数字(2^m)))
%Y参考A004094,A341713。
%Y数字m,使k^m反转为素数:此序列(k=2)、A350441(k=4)、A058993(k=5)、A058 994(k=7)、A350 442(k=8)、A05 8995(k=13)。
%K基,nonn,更多
%O 1,2号机组
%A _G.L.Honaker,Jr.,2000年10月23日
%E更多来自Chris Nash(Chris_Nash(AT)prodigy.net)的条款,2000年10月25日
%E Robert G.Wilson v_,2001年1月29日的另外两个条款
%E 2004年8月5日,来自F arideh Firoozbakht的另外3个术语
%E a(33)-a(34)摘自Giovanni Resta,2013年2月22日
|