F： =程序（d）局部A0、A4、Res、q、r；
Res:=NULL；
q： =（10^（d+1）-1）/9；
对于组合中的A0:-powerset（{$1..d-1}）do
对于组合A4:-powerset（{$1..d-1}减去A0）do
r： =q-加法（10^i，i=A0）+3*加法（10 ^i，i=A4）；
如果isprime（r）和isprime
Res:=Res，r
fi（菲涅耳）
od od；
雷斯
结束进程：
排序（[seq（F（d），d=1..7）]）#罗伯特·伊斯雷尔2018年5月3日
