Advertisement
Guest User

euler-pythontr

a guest
Jan 29th, 2015
175
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 0.82 KB | None | 0 0
  1. class xseq:
  2.     def isprime(self, n):
  3.         if n == 2 or n == 3: return True
  4.         if n < 2 or n%2 == 0: return False
  5.         if n < 9: return True
  6.         if n%3 == 0: return False
  7.         r = int(n**0.5)
  8.         f = 5
  9.         while f <= r:
  10.             if n%f == 0: return False
  11.             if n%(f+2) == 0: return False
  12.             f +=6
  13.         return True
  14.  
  15.     def primes(self, inset=0, offset=0):
  16.         list = []
  17.         for i in range(inset, offset):
  18.             if self.isprime(i):
  19.                 list.append(i)
  20.         return list
  21.  
  22. ### Run ###
  23.  
  24. eul = 600851475143
  25. eul2 = 13195
  26. primelist = myseq.primes(1,int(eul/2))
  27. cont = []
  28.  
  29. for i in range(1, len(primelist)):
  30.     if (eul % primelist[i]) == 0:
  31.         eul = eul / primelist[i]
  32.         cont.append(primelist[i])
  33.         i -= 1
  34.  
  35. print(cont)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement