2,3,5,7,11,13,.... 是素数序列。 类似:7,37,67,97,127,1577,37,67,97,127,157 这样完全由素数组成的等差数列,叫等差素数数列。
上边的数列公差为 3030,长度为 66。
20042004 年,格林与华人陶哲轩合作证明了:存在任意长度的素数等差数列。 这是数论领域一项惊人的成果!
有这一理论为基础,请你借助手中的计算机,满怀信心地搜索:
长度为 1010 的等差素数列,其公差最小值是多少?
一开始直接暴力求解,直接跑不出来,参考答案之后,发现他寻找素数的时候不是返回这个数,而是返回T/F(1/0),是素数这个位置是1,不是这个位置是0,巧思。
def prime(x):
if(x==1 or x==2):
return 0
else:
for i in range(2,x):
if(x%i==0):
return 0
if(i*i>x):
break
return 1
X=[prime(i) for i in range(1,5000)]
for d in range(2,500): #公差d开始循环
for i in range(5000): #首项i
for j in range(10): #循环10个数
try:
if(X[i+d*j]==0):
break
else:
if(j==9):
print(d)
else:
pass
except:
pass
几个注意点:
①判断素数少判断一半的方法:i*i>x 就不用看了
②19行 i+d*j别忘了这个乘以j 不然就一直是一个数字在循环了
③ try except 是为了防止X[i+d*j]出现越界