テートブリテン

素晴らしい晴天(夕方からは雨だったけど）。テート・ブリテンへ。バスで5分乗り換えなしは魅力的な目的地。
財布を持っていくのを忘れてコーヒーも飲めなかった。

view plainprint?

1. #!/usr/bin/python  
2.   
3. #Project Euler Problem 27  
4. #Start 01/May/2010  
5. #End   01/May/2010  
6.   
7. import math  
8.   
9. def is_prm(i):  
10.     if i<0:  
11.         return 0  
12.     else:  
13.         sq_i=math.ceil(math.sqrt(i))  
14.         k=2  
15.         while k<=sq_i:  
16.             if i%k==0:  
17.                 return 0  
18.             k+=1  
19.         return 1  
20.   
21. def fun(a,b,n):  
22.     return math.pow(n,2)+a*n+b  
23.   
24. prm=[]  
25. prm.append(1)  
26. for i in range(2,1000):  
27.     k=is_prm(i)  
28.     if k==1:  
29.         prm.append(i)  
30. print prm  
31. mx=0  
32. for a in range(-1000,1001,1):  
33.     for b in prm:  
34. #        print a,b  
35.         n=0  
36.         while 1:  
37.             k=is_prm(fun(a,b,n))  
38.             if k==1:  
39.                 n+=1  
40.             else:  
41.                 break  
42.         if n>mx:  
43.             mx=n  
44.             print a,b,n,a*b  

#!/usr/bin/python

#Project Euler Problem 27
#Start 01/May/2010
#End   01/May/2010

import math

def is_prm(i):
    if i<0:
        return 0
    else:
        sq_i=math.ceil(math.sqrt(i))
        k=2
        while k<=sq_i:
            if i%k==0:
                return 0
            k+=1
        return 1

def fun(a,b,n):
    return math.pow(n,2)+a*n+b

prm=[]
prm.append(1)
for i in range(2,1000):
    k=is_prm(i)
    if k==1:
        prm.append(i)
print prm
mx=0
for a in range(-1000,1001,1):
    for b in prm:
#        print a,b
        n=0
        while 1:
            k=is_prm(fun(a,b,n))
            if k==1:
                n+=1
            else:
                break
        if n>mx:
            mx=n
            print a,b,n,a*b