N： =200:#使用所有（p，q）和p*q<N
素数：=选择（isprime，[$2.floor（N/2）]）：
A： ={}：
对于i从1到nops（素数）do
p： =素数[i]；
Qs:=选择（q->q<N/p，[seq（素数[j]，j=1..i-1）]）；
A： =并集{seq（（2^（p*q）-1）/（2^p-1）/（2 ^q-1），q=Qs）}；
日期：
#枫树12及以上
选择（isprime，A）；
#或在早期的Maple版本中
排序（[select（isprime，A）；#_）[]]）[]；
