filter：=proc（n）局部i，j，count，d，S，x，y；
如果不是isprime（n），则返回false fi；
d： =ilog10（n）+1；
计数：=0；S： ={}；
对于i从0到d-1 do
x： =地板（n/10^i）；
对于j从i到d-1 do
y： =x模10^（j-i+1）；
如果不是成员（y，S）和isprime（y），则计数：=计数+1；S： =S并集{y}；如果计数>2*d，则返回真fi-fi
od od；
假
结束进程：
选择（过滤器，[seq（i，i=1..10^6，2）]）#罗伯特·伊斯雷尔2020年11月11日
