%I#12 2020年12月23日07:26:30

%S 0,0,0,1,2,2,2,2,7,8,3,9,7,9,9,8,6,14,3,14,18,9,14,37,26,16,28,24,13，

%电话：40,7,45,42,11,20,54,27,35,37,15,43,82,22,13,99,49,24,66,32,54,90,67，

%U 32123,90,86141,55,53142,49,36164,24,96114,4811413,92118133,83,79198,40号

%N a（N）是素数对p<q的（q mod p）之和，使得p+q=2*N。

%H Robert Israel，n的表，a（n）表示n=1..10000</a>

%F a（n）=和{k=1..n}（（2*n-k）modk）*c（k）*c（2*n-k），其中c是素特征（A010051）_韦斯利·伊万·赫特，2020年12月22日

%e a（9）=7，因为2*9=5+13=7+11，其中（13 mod 5）+（11 mod 7）=3+4=7。

%p N:=200:#对于a（1）。。a（N/2）

%p p：=选择（isprime，[seq（i，i=3..N，2）]）：

%p nP:=无（p）：

%p V：=矢量（N）：

%p代表i从1到nP do

%p代表j从i+1到nP do

%p v:=p[i]+p[j]；

%p如果v>N，则打破fi；

%pV[V]：=V[V]+（p[j]修改p[i]）

%操作说明：

%p[seq（V[2*i]，i=1..N/2）]；

%t表[Sum[Mod[2 n-k，k]（PrimePi[k]-PrimePi[k-1]）（PrimePi[2 n-k]-PriemPi[2 n-k-1]），{k，n}]，{n，60}]（*Wwesley Ivan Hurt_2020年12月22日*）

%Y参考A010051。

%K nonn，看

%O 1,4型

%A._J.M.Bergot_和_Robert Israel_，2020年12月20日