所有可能的因素对

时间:2012-08-14 04:04:05

标签: python math

我正在尝试编写一个程序,它接受带有一个前导整数和任意数量的尾随零的数字,然后打印出两个因子的所有可能组合。

100

因子是2 ^ 2,5 ^ 2

所以程序会打印出来:

(2,50),(4,25),(5,20)

600

因子是2 ^ 3,3,5 ^ 2

(2300),(4150),(8.75),(3200),(5120),(25,24),(6100),(12,50),(15,40),(30, 20),(60,10)

......我想这就是全部吗?是吗?我可以使用程序检查......

import itertools

facts=[[2,2,2],[3],[5,5]]
for n in itertools.product(*facts)
    print(n)

我发现我使用的不正确但这是我第一次尝试它。

这只给了(2,3,5)十次。

我想要像(2)*(2,3,5,5)和(2,2)*(3,5.5)之类的东西......

4 个答案:

答案 0 :(得分:2)

这是我将如何做到的:

def factors(n):
  # Fill this in

def factor_pairs(n):
  for i in factors(n):  # You need to write the factor() function
    yield i, n / i

if __name__ == '__main__':
  n = input('Enter an integer: ')

  for i, j in factor_pairs(n):
    print i, j

我不会完全为你编码,但你明白了。

答案 1 :(得分:2)

To generate all factors of a number given its prime factors

#!/usr/bin/env python
import itertools, operator

def all_factors(prime_dict):
    series = [[p**e for e in range(maxe+1)] for p, maxe in prime_dict.items()]
    for multipliers in itertools.product(*series):
        yield reduce(operator.mul, multipliers)

实施例

prime_dict = {2:3, 3:1, 5:2}
L = sorted(all_factors(prime_dict))
number_of_divisors = reduce(lambda prod, e: prod*(e+1), prime_dict.values(),1)
assert len(L) == number_of_divisors
# -> [1, 2, 3, 4, 5, 6, 8, 10, 12, 15, 20, 24,
#     25, 30, 40, 50, 60, 75, 100, 120, 150, 200, 300, 600]

生产配对:

n, isodd = divmod(len(L), 2)
print(zip(L[:n], reversed(L[n + isodd:])))
if isodd: # number is perfect square
   print((L[n], L[n]))

输出

[(1, 600), (2, 300), (3, 200), (4, 150), (5, 120), (6, 100),
 (8, 75), (10, 60), (12, 50), (15, 40), (20, 30), (24, 25)]

适用于少数人。您可以使用它来测试您的解决方案,该解决方案可以考虑您的数字的特殊形式:x00000...

答案 2 :(得分:0)

我认为这会做你想要的:

n = input('Enter a number? ')
factors = []

for i in range(int(sqrt(n))):
  if n % i == 0:
    factors.append((i, n / i))

根据因子的定义,您必须检查的最大数量取决于数字的平方根,因此如果您可以对此进行编码,则应设置。

如果你可以编码,你应该设置。

答案 3 :(得分:0)

你可以把它全部放在列表理解中

import math
n = 600 # or whatever...
[(x, n/x) for x in range(1, int(math.sqrt(n))+1) if n % x == 0]

返回:

[(1, 600), (2, 300), (3, 200), (4, 150), (5, 120), (6, 100), (8, 75), (10, 60), (12, 50), (15, 40), (20, 30), (24, 25)]

如果你不想要(1,600)只使用范围(2,int(math.sqrt(n))+ 1)。