#欧拉考虑了n^2+n+41的素数值；
#这是2.93 GHz机器上76秒的计算
h:=n^2+n+41；
y:=c*a*z^2+（（（（d+2）*（1/3））*c-2）*a/d+1）*z+（（（367*d^2+d+1）*；
y2：=子（n=y，h）；
y3:=系数（y2）；
#注意，y是4个变量中的表达式。
#在合成函数后，进行代数因式分解
#可以在y3中观察到。只要y3是一个整数，它就会
#是复合的。这是因为y3因素和这两个因素
#是大于1的整数。
最大值：=6000；
答：={}：
对于n到maxn do
克：=n^2+n+41：
如果isprime（g）=false，则A:=`union`（A，{n}）结束，如果：
结束do：
#现在，A集包含表单的组合值
#n^2+n+41小于maxn。
c:=1:a:=1:d:=1:z:=-1:p:=41:
q:=c*a*z^2+（（（（d+2）*（1/3））*c-2）*a/d+1）*z+（（（367*d^2+d+1）*（1/9）
A2：=A：
而q＜maxn do
而`和`（q<maxn，d<100）确实
while q＜maxn做while
q<最大do
A2:=`减号`（A2，{q}）；
A2：=“减去”（A2，{c*a*z^2+（（（（d+2）*（1/3））*c-2）*a/d+1）*z+（（（（367*d^2+d+1）*（1/9））*c^2-（（d+2）*（1/3））*c+1）*a/d^2-（（（（d-1）*（1/3））*c+1）/d）/c}）；
z:=z+1；
A2:=`减号`（A2，{c*a*z^2-（（（（d+2）*（1/3））*c-2）*a/d+1）*（1*z）+；q:=c*a*z^2+（（（（d+2）*（1/3））*c-2）*a/d+1）*z+（（（367*d^2+d+1）*（1/9）
结束do；
a:=a+1；z:=-1；
q:=c*a*z^2+（（（（d+2）*（1/3））*c-2）*a/d+1）*z+（（（367*d^2+d+1）*（1/9）
结束do；
d：=d+1：a：=1：
q:=c*a*z^2+（（（（d+2）*（1/3））*c-2）*a/d+1）*z+（（（367*d^2+d+1）*（1/9）
结束do：
c：=c+1：d：=1：
q:=c*a*z^2+（（（（d+2）*（1/3））*c-2）*a/d+1）*z+（（（367*d^2+d+1）*（1/9）
结束do：
A2；
