%I#39 2021年1月24日10:29:57
%S 1,3,5,4,2,6,8,7,9,11,10,12,13,14,16,17,15,18,20,19,21,23,22,24,25,26,
%电话28,27,29,31,32,30,33,35,34,36,38,37,39,41,40,42,43,44,46,45,47,48,50,
%U 49、51、53、52、54、56、55、57、58、59、60、61、62、63、65、64、66、67、68、70、71、69、72
%N a(N)+a(N+1)决不是素数;在词典学上最早的这种不同的正整数序列。
%C基本相同的序列见A253074(并证明该序列是置换)。
%C序列A253074的定义方式相同,但以0开头。这恰好从下一个项开始产生相同的序列。这是序列族中的情况(M,N)=(2,0),其中M个连续项在其成对和中产生N个素数,请参阅wiki页面了解其他示例。-_M.F.Hasler,2019年11月26日
%H Robert Israel,n的表,n=1..10000的a(n)</a>
%H M.F.Hasler,邻项质数和,OEIS wiki,2019年11月23日
%H<a href=“/index/Per#IntegerPermutation”>自然数排列序列的索引项</a>
%F a(n)=A253074(n+1)(只要A25307S(1)=0)_M.F.Hasler,2019年11月26日
%e a(3)=5,因为已经使用了1和3,3+2=5和3+4=7都是质数,而3+5=8不是质数。
%p编号:=1000;#得到n到n的[n]
%pA:={1};
%pa[1]:=1;
%p代表n从2到n do
%p mA:=最大值(A);
%p R:={$1..mA}减去A;
%R do中x的p
%p如果不是isprime(a[n-1]+x),则
%p a[n]:=x;
%p断点
%功率因数
%日期:
%p如果未赋值(a[n]),则
%p代表mA+1 do中的x
%p如果不是i素数(a[n-1]+x),则
%p a[n]:=x;
%p断点
%功率因数
%日期
%p fi;
%p A:=并集{x};
%日期:
%p序列(a[n],n=1..n);#_罗伯特·伊斯雷尔(Robert Israel),2014年6月3日
%tf[s_]:=块[{k=1,a=s[[-1]]},而[Or[MemberQ[s,k],素数Q[a+k]],k++];追加[s,k]];嵌套[f,{1},121](*_Zak Seidov_,2009年10月21日*)
%tα={1};z=范围[2202];z=补码[z,a];当[Length[z]>1时,如果[!PrimeQ[z[[1]]+Last[a]],AppendTo[a,z[[1]]],如果[!PrimeQ[z[[2]]+Last[a]],AppendTo[a,z[[2]]],AppendTo[a,z[[3]]];z=补码[z,a]];打印[a](*速度明显加快*)(*_Vladimir Joseph Stephan Orlovsky_,2011年5月3日*)
%o(哈斯克尔)
%o导入数据。列表(删除)
%o a055266 n=a055266_列表!!(n-1)
%o a055266_list=1:f 1[2..]其中
%o f u vs=g vs其中
%o g(w:ws)|a010051'(u+w)==0=w:f w(删除wvs)
%o |否则=g ws
%o--_Reinhard Zumkeller_,2015年1月14日
%o(PARI)v=[1];n=1;while(n<100,if(!i素数(n+v[#v])&&!vecsearch(vecsort(v),n),v=concat(v,n);n=0);n++);v\\_Derek Orr_,2015年6月8日
%o(PARI)A055266_upto(n=99,u=1,u,a)={向量(n,n,n=u;while(bittest(u,n-u)||isprime(a+n),n++);if(n>u,u+=1<<(n-u),u>>=-u+u+=估值(u+2,2));a=n)+if(默认(调试),打印([u]))))}\\可选参数允许调整计算。如果debug>0,则在末尾打印最少未使用的数字。-_M.F.Hasler,2019年11月25日
%Y参见A055265、A010051、A249920(反向)、A203069、A253074。
%K容易,不是
%O 1,2号机组
%2000年5月9日,安利底特律
%E由_Zak Seidov修正,2009年10月21日
%E Name由M.F.Hasler_编辑,2019年11月26日