blob: 7c1646bf24db7b3fa110925d3402391c14a4222b (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
import PIL, math
import numpy as np
#Problem 10 sum of primes
###INEFFICIENT###
#def isPrime(number):
# if (number <=1):
# return False
# for i in range(2,number):
# if(number%i==0):
# return False
# return True
#def prime_factors(number):
#primes = []
#for i in range(number, 2, -1):
#if(number%i==0):
#if (isPrime(i)):
#return i
#return primes
def prime_factors(number):
primes = []
if (number <=1):
return [1]
a = np.ones(number)
a[0]=0
a[1]=0#not prime
for i in range(2,int(math.ceil(math.sqrt(number)))):
if(a[i]==1):
j = i**2 #cross out all the multiples
while(j < number):
a[j] = False
j+=i
for k in range(2, number):
#is prime
if(a[k]==1):
primes.append(k)
return primes
out = prime_factors(2000000)
print str(sum(out))
|
