（Python）
从同情进口到犯罪；从数学导入sqrt，ceil
定义neib（m）：
如果m==1：返回[3，5，7，9，11，13]
如果m==3：返回[17,19,5,13,15]
L=[m，对于范围（6）中的i）]；n=int（ceil（（3+sqrt（6*m+3））/6））；x=6*n*n；y=12*n
a0=x-18*n+15；a1=x-16*n+11；a2=x-14*n+9
a3=x-y+7；a4=x-10*n+5；a5=x-8*n+3；a6=x-6*n+1
p＝0如果m＝＝a0 else 1如果m＞a0并且m＜a1 else 2如果m＝＝a1 else 3如果m＞a1并且m＜a2 else 4如果m＝＝a2 else 5如果m＞a2并且m＜a3 else 6如果m＝＝＝a3 else 7如果m＞a3并且m＜a4 else 8如果m＝＝＝a4 else 9如果m＞a4并且m＜a5 else 10如果m＝＝＝＝＝a5 else 11如果m＞a5并且m＜a6 else 12
如果p<=4，则L[0]+=y-10；如果p<=6，则-2；如果p>=9，则y+16；否则2
如果p<=1，则L[1]+=2；如果p<=6，则y-8；如果p>=8，则-2；如果p=8，则y-14
如果p<=1，L[2]+=-y+24；如果p<=3，则为2；如果p<0=8，则为y-6；如果p>=10，则为-2；否则为-y+12
如果p==0，则L[3]+=-2；否则-y+22，如果p<=3，则为2；如果p<=5，则为2-4，如果p>=10，则为-2
L[4]+=y-14，如果p==0，否则-2，如果p<=2，否则-y+20，如果p>=5，否则2，如果p<0=7，否则y-2
如果p<=2，则L[5]+=y-12；如果p<=4，则为-2；如果p=7，则为y+18；如果p>=9，则为2
返回L
对于范围（11500）内的i：
m=2*i-1
如果isprime（m）==1：
L1=[neib（m）[j]，对于范围（6）中的j
如果总和（L1中k的i素数（k））==0：打印（m）
