Python实现欧拉计划27题

来源:互联网 发布:mac怎么下photoshop 编辑:程序博客网 时间:2024/06/04 22:30
import math


from time import clock
time1=clock()


def IsPrime(x):  
    if x<2:  
        return False  
    for i in range(2,int(math.sqrt(x))+1):  
        if x%i==0:  
            return False  
    return True


def func(a,b):  
    k=0  
    while True:  
        if IsPrime(k*k+a*k+b):  
            k+=1  
        else:  
            break  
    return k-1  
  
  
maxa,maxb=0,0  
num=0  
for j in range(-1000,1000):  
    if IsPrime(j):    #巧妙,缩短了三秒的时间。k(k+a)+b b肯定为质数
        for i in range(-1000,1000):  
            temp=func(i,j)  
            if temp>num:  
                num=temp  
                maxa,maxb=i,j
print(maxa*maxb)


time2=clock()
time=time2-time1
print 'it takes %f seconds'%time
0 0