revdigs:=proc（n）局部L，i；五十： =换算（n，基数，10）；添加（L[-i]*10^（i-1），i=1..nops（L））结束进程：
isemirp1：=进程（n）局部r；
如果不是isprime（n），则返回false fi；
r： =转速（n）；
r>n与isprime（r）
结束进程：
E： =选择（isemirp1，[seq（seq（i*10^d+j，j=1..10^d-1，2），i=[1,3,7,9]），d=1..5）]）：
五： =排序（映射（t->t+revdigs（t），E））：
N： =无（V）：
W： =矢量（16）：
i： =1：
而i<N do
对于从1到N-i的j，当V[i+j]=V[i]做od时：
如果j<=16且W[j]=0，则W[j]：=V[i]fi；
i： =i+j；
日期：
转换（W，列表）；
