数因式分解

Question

任务是将foldes数编写成素数函数。通过给定的数字'n‘，该函数应该返回元组p_i，c^i的列表，例如，如果输入为100，则输出为(2,2)，(5,2)。所以，下面是我尝试写它的方法：

def factor(n):
c = 1
pre_ans = list()
temp_n=n
for i in range(2,temp_n+1):
    if (is_prime(i) == True) and (temp_n % i == 0):
        for j in range (2,temp_n+1):
            if (temp_n % (i ** j) == 0):
                pre_ans.append((i,j))
                temp_n /= (i **j)
        pre_ans.append((i,c))
        temp_n /= i
print(pre_ans)

这是错误的，但我找不到一个错误

Answer 1

修正了密码。这是工作版本

def factor(n):
c = 1
pre_ans = list()
temp_n=n
for i in range(2, n // 2 + 1):
    if (is_prime(i) == True):
        k = 1
        while temp_n % (i ** k) == 0:
            if temp_n % (i ** (k + 1)) == 0:
                k += 1
            else:
                k += 1
                break
        if k > 1:
            pre_ans.append((i, k - 1))
return pre_ans

Answer 2

你的大致想法是可以的。但是，您的代码的以下部分存在一些小问题：

for j in range (2,temp_n+1):
    if (temp_n % (i ** j) == 0):
        pre_ans.append((i,j))
        temp_n /= (i **j)
pre_ans.append((i,c))
temp_n /= i

实际上，主要的问题是您需要在这个语句for j in range (2,temp_n+1)中向相反的方向迭代。如果你把它重写成

def factor(n):
    c = 1
    pre_ans = list()
    temp_n=n
    for i in range(2,temp_n+1):
        if (is_prime(i) == True) and (temp_n % i == 0):
            for j in range (temp_n+1, 0,-1):
                if (temp_n % (i ** j) == 0):
                    pre_ans.append((i,j))
                    temp_n /= (i **j)
    print(pre_ans)

看起来不错。整个代码也可以编写得更短一点：

from collections import Counter

def factor(n):
    lst = []
    for i in range(2, n+1):
        while n % i == 0:
            lst.append(i)
            n = n / i
    return Counter(lst).items()

print(factor(100))