%I#21 2023年1月18日19:21:37
%S 1110110011100111101110111110101111100011100011111100101010011,
%电话:110101010101111101101101111111111111110111100110101011001,
%U 10100001001010100010110101101110000011100111011011100011
%N只包含数字0和1的回文素数。(A020449中的回文术语。)
%H Jon E.Schoenfield,n的表,n的a(n)=1..15185(所有术语<10^38;T.D.Noe的术语1..200,Chai Wah Wu的术语201..2385)
%e a(3)=100111001,因为100111001是仅由1和0组成的第三个回文素数。
%t PalQ[n_]:=起始数字[反转@整数位数[n] ]==n;删除重复项[Flatten[Table[Select[FromDigits/@Tuples[{0,1},n],PrimeQ[#]&&PalQ[#]&],{n,15}]](*_Jayanta Basu_,2013年5月11日*)
%t选择[FromDigits/@Tuples[{0,1},15],PalindromeQ[#]&&PrimeQ[#]&](*哈维·P·戴尔,2023年1月18日*)
%o(Python)
%o来自sympy import isprime
%o A100580_列表=[11]
%o对于范围(2,2**16)内的i:
%o。。。。s=格式(i,'b')
%o。。。。x=整数(s+s[-2::-1])
%o。。。。如果是质数(x):
%o。。。。。。。。A100580_list.append(x)#_Chai Wah Wu_,2015年1月6日
%Y参考A020449、A002385、A000040、A002113。
%K基,nonn
%O 1,1
%A Chuck Seggelin(seqfan(AT)plasteddragon.com),2004年11月29日
|