5

我写了一个整数分解函数，但是在弄乱了它之后，我意识到它有几个数字的问题......

>>> pFactors(99) # it does work for numbers with multiple of one prime factor
[3, 3, 11]
>>> pFactors(999) # however, sometimes, it doesn't
[3, 37] # the actual prime factorization of 999 is [3, 3, 3, 37]. 
>>> pFactors(88)
[2, 11]
>>> pFactors(888)
[2, 3, 37]

我的代码有什么问题？

def pFactors(n):
   """Finds the prime factors of 'n'"""
   from primes import isPrime
   from math import sqrt
   pFact, limit, check, num = [], int(round(sqrt(n), 2)) + 1, 2, n
   if isPrime(n):
      return [n]
   for check in range(2, limit):
      if isPrime(check) and num % check == 0:
         pFact.append(check)
         num /= check
         if isPrime(num):
            pFact.append(num)
            break
   pFact = sorted(pFact)
   return pFact
4

2 回答 2

10

像这样修改：

def pFactors(n): 
        """Finds the prime factors of 'n'""" 
        from math import sqrt 
        pFact, limit, check, num = [], int(sqrt(n)) + 1, 2, n 
        if n == 1: return [1] 
        for check in range(2, limit): 
             while num % check == 0: 
                pFact.append(check) 
                num /= check 
        if num > 1: 
          pFact.append(num) 
        return pFact 

for i in range(1,1000):
        print pFactors(i)

尽管我喜欢您最初编写的代码，但有几点：

  1. 您不需要 isPrime。原因是任何在限制范围内除以 num 的素数也将是除以 num 的任何复合的最小除数，因此当您除掉这些素数时，您将阻止它们组成的复合被发现为范围后面的除数，只剩下主要因素。

  2. 您不需要对数组进行排序，它已经通过check升序排序。

  3. 添加的while循环确保只要它们继续除num，就可以正确找到重复因子。

  4. 您可以使用同余过滤掉小于限制的所有数字的 2/3 以检查为除数，您知道怎么做吗？

上面结果的最后几行是：

[11, 89]
[2, 2, 5, 7, 7]
[3, 3, 109]
[2, 491]
[983]
[2, 2, 2, 3, 41]
[5, 197]
[2, 17, 29]
[3, 7, 47]
[2, 2, 13, 19]
[23, 43]
[2, 3, 3, 5, 11]
[991]
[2, 2, 2, 2, 2, 31]
[3, 331]
[2, 7, 71]
[5, 199]
[2, 2, 3, 83]
[997]
[2, 499]
[3, 3, 3, 37]
于 2013-01-27T19:05:24.177 回答
1

在 999 示例中，将其除以 3 后，结果 (333) 也可除以 3，即得到 111，您也可以除以 3（得到 37）。但是在您的代码中，您只需将其除一次 - 您需要做的是，一旦您找到一个除以您当前数字的素数，就是继续除以该数字，直到它不再可能。

于 2013-01-27T18:48:36.077 回答