python蓝桥杯训练营一 等差素数列

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]出现越界

你可能感兴趣的:(python蓝桥杯训练营,python,蓝桥杯)