%I#31 2024年5月22日17:27:30
%S 2,3,5,7,13,17,41,73,89,97137193257281313409449521569577,
%电话:6176416737617699299531031049115312491289140916011657,
%电话:169717211801191320812113215322972441259327133137325733613449
%N个素数,其二进制表示不是较小素数的二进制表示的串联(允许前导0)。
%C类似于A090422,但允许在任何素数的表示中使用前导零。例如,碱基2中的19是10011,可以写成(10)(011),因此不出现在该序列中(但出现在A090422中)。
%C经验上,从a(6)=17开始,a(n)==1(mod 8)。-_雨果·普福尔特纳,2021年3月6日
%C这个观察结果源于这样一个事实:与前6个素数相对应的正则表达式(0*10+0*11+0*101+0*111+0*101/0*1101)*有一个补语,它只包含1、01、一些以0结尾的单词和一些以001结尾的单词_杰弗里·沙利特,2021年3月7日
%H Robert Israel,n的表,n的a(n)=1..10000</a>
%p CSP:=proc(n)选项记住;局部g;
%p g:=proc(k)局部v;v: =n模2^k;isprime(floor(n/2^k))和(isprim(v)或CSP(v))end proc;
%价格(g,[2美元..ilog2(n)])
%p端程序:
%p CSP(0):=假:
%p删除(CSP,[seq(ithprime(i),i=1..1000)]);#_罗伯特·伊斯雷尔,2024年5月22日
%o(Python)
%o来自sympy import isprime,primerange
%o定义正常(p):
%o b=箱(p)[2:]
%o对于范围(2,len(b)-1)中的i:
%o如果是i素数(int(b[:i],2)):
%o如果isprime(int(b[i:],2))或不正常(int(b[i:]2)):返回False
%o返回True
%o定义aupto(lim):如果确定(p),返回[p代表素数范围(2,lim+1)中的p]
%o打印(aupto(3449))#_Michael S.Branicky_,2021年3月7日
%Y参考A090422。
%K nonn,基础
%O 1,1号机组
%A _杰弗里·沙利特_,2021年3月7日
|