テートブリテン

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

#!/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