N： =10000:#获取所有术语<=N
P： =选择（isprime，{seq（P，P=3..2*N，2）}）：
计数：=1:
A[计数]：=2：
对于从1开始的ic，而P[ic]<=N do
c： =P[ic]；
五： =地图（`-`，P[ic+1..-1]，c）相交地图（t->c-t，P[1..ic-1]）；
nv:=nops（V）；
VV:={seq（seq（V[j]-V[i]，j=i+1..nv），i=1..nv-1）}与V相交；
nvv:=nops（VV）；
发现：=错误；
对于从1到nvv的ia，当找不到时，执行do
a： =VV[ia]；
当VV[ib]<c-a do时，ib从ia+1到nvv
b： =VV[ib]；
如果b<>2*a和{c-a-b，c-a+b，c-b+a，c+a+b}子集P，则
发现：=true；
打破
fi（菲涅耳）
操作系统
日期：
如果找不到，则
计数：=计数+1；
A[计数]：=c；
fi（菲涅耳）
日期：
seq（A[i]，i=1..计数）#罗伯特·伊斯雷尔2016年2月16日
