N： =30:#对于最多2*N个数字的术语
R： ={1}：T：={1，9}：
对于从2到N的d do
T： =选择（T->转换（转换（T^2 mod 10^d，base，10），set）子集{0，1，4，9}，映射（T->（T，T+10^（d-1），T+4*10^；
R： =R联合选择（t->转换（转换（t^2，base，10），set）子集{0，1，4，9}，t）；
日期：
R2:=映射（t->t^2，R）：
Res:=映射（t->seq（t*10^（2*i），i=0..（2*N-ilog10（t）-1）/2），R2）并集｛0｝：
排序（convert（Res，list））#罗伯特·伊斯雷尔2024年8月26日
