Guest User

PrimeFactorComparison

a guest
Aug 15th, 2018
165
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
Python 1.48 KB | None | 0 0
  1. import numpy as np
  2. import math
  3.  
  4.  
  5. def primeFactors(n):
  6.     list=[]
  7.     while n % 2 == 0:
  8.         #print 2,
  9.         n = n / 2
  10.          
  11.     for i in range(3,int(math.sqrt(n))+1,2):
  12.  
  13.         while n % i== 0:
  14.            # print i,
  15.             n = n / i
  16.             list.append(i)
  17.  
  18.     if n > 2:
  19.       #  print n
  20.         list.append(n)
  21.    
  22.     return list
  23.  
  24.  
  25. #digit prime count
  26. def size_of_primefactors(n):
  27.     Investigate = False
  28.    
  29.     size = len(str(i))
  30.     if size>=3:
  31.         Investigate = True
  32.        
  33.     return Investigate
  34.  
  35. def expirement():
  36.     RGL_one=[]
  37.     for i in range(294):
  38.         RGL_one.append(np.random.randint(4,846))
  39.  
  40.     nine_RGL_one=[]
  41.     for i in range(294):
  42.         nine_RGL_one.append(RGL_one[i]*19)
  43.  
  44.     RGL_two=[]
  45.     sums=0
  46.     for i in range(294):
  47.         sums+=np.random.randint(1,141)
  48.         RGL_two.append(sums)
  49.  
  50.     lv=np.sum(RGL_two)
  51.     nv=np.sum(nine_RGL_one)
  52.     lvp=primeFactors(lv)
  53.     nvp=primeFactors(nv)
  54.    
  55.     SamePrime=[0]
  56.     SamePrimeFN=False
  57.  
  58.     for i in range(len(lvp)):
  59.         for j in range(len(nvp)):
  60.             Size= True
  61.             Size=size_of_primefactors(lvp[i])
  62.             if (lvp[i]==nvp[j]):
  63.                     #print "Same PrimeF, %3d"%(lvp[i])
  64.                 SamePrime.append(lvp[i])
  65.                
  66.    
  67.     return SamePrime
  68.    
  69. for i in range(1900000):
  70.     SamePrime=expirement()
  71.     for j in SamePrime:
  72.         if j>10000:
  73.             print SamePrime
Advertisement
Add Comment
Please, Sign In to add comment