(PARI)列表(lim)={
my(v=列表([5616]),t);
对于素数(p=2,log(lim)\log(8)+2,
列表输入(v,(4^p-1)<<p)
);
对于素数(p=3,log(2*lim)\log(27)+2,
列表输入(v,3^p*(9^p\2))
);
对于素数(p=3,log(lim)\log(32)+2,
列表输入(v,(4^p-1)*(2^p-1,<<(2*p))
);
对于素数(p=7,sqrtn(2*lim,3)+1,
如果(p%5>1&&p%5<4,listput(v,p^2\2*p))
);
vecsort(选择(n->n<=lim,Vec(v)))
};