%I#17 2016年5月14日13:46:04
%S 0,0,0,1,1,2,2,3,4,5,6,7,8,9,10,11,11,12,13,14,15,16,16,17,17,
%电话:18,19,20,21,21,212,22,23,24,25,26,27,28,28,29,29,30,31,32,33,34,
%U 34,35,36,37,38,38,39,40,41,42,43,44,45,46,47,48,49,50
%N与非赢素数的累积数有关。
%H Robert Israel和Paul Nanninga,n表,a(n)表示n=0..100000
%取奇素数序列op(n);设置a(0)=0;如果op(n+1)-op(n)=2a(n+1。
%p N:=100:#得到a(0)到a(N)
%p p:=[seq(ithprime(i),i=2..N+2)]:
%p列表工具:-部分和([0,seq(`if`(p[i]-p[i-1]=2,0,1),i=2..N+1)]);#_罗伯特·伊斯雷尔,2016年5月13日
%t累加@表[If[Prime@n-Prime[n-1]==2,0,1],{n,2,120}]-1(*_Michael De Vlieger_,2016年5月13日,在_Robert Israel_*之后)
%o(PARI)op(n)=素数(n+1);
%o lista(nn)={my(x=0);对于(n=1,nn,print1(x,“,”);如果((op(n+1)-op(n))>2,x++););}\\_Michel Marcus_,2016年5月13日
%Y参见A065091(奇数素数),A007510(非双素数)。
%K容易,不是
%0、6
%A Rudi Huysmans(鲁迪·休斯曼(AT)pandora.be),2002年9月18日
|