Project Euler 7-10

7>判斷第10001個素數是多少？

from math import sqrt def isprime(number): flag=True begin=2 end=sqrt(number) while begin<=end: if number%begin==0: flag=False break begin+=1 return flag primes=[] begin=2 while len(primes)<=10001: if isprime(begin): primes.append(begin) print begin begin+=1

8>求1000個數中連乘積最大的

s=""" 73167176531330624919225119674426574742355349194934 96983520312774506326239578318016984801869478851843 85861560789112949495459501737958331952853208805511 12540698747158523863050715693290963295227443043557 66896648950445244523161731856403098711121722383113 62229893423380308135336276614282806444486645238749 30358907296290491560440772390713810515859307960866 70172427121883998797908792274921901699720888093776 65727333001053367881220235421809751254540594752243 52584907711670556013604839586446706324415722155397 53697817977846174064955149290862569321978468622482 83972241375657056057490261407972968652414535100474 82166370484403199890008895243450658541227588666881 16427171479924442928230863465674813919123162824586 17866458359124566529476545682848912883142607690042 24219022671055626321111109370544217506941658960408 07198403850962455444362981230987879927244284909188 84580156166097919133875499200524063689912560717606 05886116467109405077541002256983155200055935729725 71636269561882670428252483600823257530420752963450 """ s=s.replace('/n','') numbers=[int(i) for i in s] max=0 i=1 while i<996: tmp=numbers[i-1]*numbers[i]*numbers[i+1]*numbers[i+2]*numbers[i+3] if tmp>max: max=tmp i+=1 print max

9>求平方數使得a,b,c之和爲1000,且c平方=a平方+b平方

a,b=1,1 for ib in range(1,1000): flag=False for ja in range(1,i): if ib**2+ja**2==(1000-ib-ja)**2: print ib*ja*(1000-ib-ja) flag=True break if flag==True: break

10>計算200萬以下素數的和,isprime函數如上

i,sum=2,0 while i<=2000000: if isprime(i): sum+=i i++