你另外帖子里面的代码：
  
可以改成下面这样，能少算几次：
flag = False
#五个素数的乘积-1
for il in range(0,1025):
     for i2 in range(i1+1,1026):
         m2 = data2[i1][0] * data2[i2][0]
         for i3 in range(i2+1,1027):
             m3 = m2 * data2[i3][0]
             for i4 in range(i3+1,1028):
                 m4 = m3 * data2[i4][0]
                 for i5 in range(i4+1,1029) :
                     五个素数的乘积 =  m4 * data2[i5][0]#print(五个素数的乘积)


计算的库可以试试 gmpy2。


【 在 feng321 的大作中提到: 】
: 我有1029个不相同的素数，然后  要从中取出若干个素数（无重复），把它们乘起来构成一个大整数 N,  要求N的每一个素因子 p 都满足条件
: p-1 整除 N-1，且  p+1 整除 N+1
: （只要计算奇数个素因子的乘积，再验算就行了）
: ...................
--
FROM 115.171.23.*