%I#16 2016年2月16日17:39:40

%S 2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,61,67,79,83,97101107，

%电话1131631811971992292332731317331433439457569859

%N个素数不在3X3个素数幻方的中心。

%C“魔和”总是中央入口的三倍。

%C没有其他<5000的术语。

%C没有其他<100000的条款_罗伯特·伊斯雷尔（Robert Israel），2016年2月16日

%H<a href=“/index/Mag#magic”>为与幻方相关的序列索引条目</a>

%p N:=10000:#获取所有术语<=N

%p p:=选择（isprime，{seq（p，p=3..2*N，2）}）：

%p计数：=1：

%p A[计数]：=2：

%当p[ic]<=N do

%p c:=p[ic]；

%pV:=地图（`-`，p[ic+1..-1]，c）相交地图（t->c-t，p[1..ic-1]）；

%p nv:=nops（V）；

%p VV:={seq（seq（V[j]-V[i]，j=i+1..nv），i=1..nv-1）}与V相交；

%p nvv:=nops（VV）；

%发现p：=false；

%p表示ia从1到nvv，而未找到do

%p a：=VV[ia]；

%p表示ib从ia+1到nvv，而VV[ib]<c-a do

%p b：=VV[ib]；

%如果b<>2*a和{c-a-b，c-a+b，c-b+a，c+a+b}子集p，则

%发现p：=true；

%p断点

%功率因数

%日期

%日期：

%p如果找不到，则

%p计数：=计数+1；

%p A[计数]：=c；

%功率因数

%日期：

%p序列（A[i]，i=1..计数）；#_罗伯特·伊斯雷尔（Robert Israel），2016年2月16日

%Y参考A073473。A024351中给出了中心为59的幻方。

%K nonn公司

%O 1,1号机组

%2002年8月25日，A·达维德·威尔逊