鱼C论坛

 找回密码
 立即注册
查看: 2556|回复: 1

欧拉计划第三题疑问

[复制链接]
发表于 2017-3-17 17:13:07 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
求600851475143的最大质因数?

n=600851475143
a=2
while n!=a:
    if n%a==0:
        n=int(n/a)
    a+=1
    if n==a:
        print("最大质数为:",a)

结果为:6857

但是如果求20的最大质因数,为什么就不对了呢?
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2017-3-17 19:31:52 | 显示全部楼层
import math


#def is_prime(x):
#    """ 判断x是否是质数 """
#    x_sqrt = int(math.sqrt(x))
#    for x_ in range(x_sqrt, 1, -1):
#        if x % x_ == 0:
#            return False
#    return True
#    
#    
#def get_factor(x):
#    """ 获取x的质数因子 """
#    if is_prime(x):
#        return [x]
#    x_sqrt = int(math.sqrt(x))
#    for x_ in range(x_sqrt, 1, -1):
#        if x % x_ == 0:
#            tmp = get_factor(x_)
#            tmp.extend(get_factor(int(x / x_)))
#            return tmp


def is_prime(x):
    """ 判断x是否是质数 """
    if x == 2:
        return True
    x_sqrt = int(math.sqrt(x))
    l = [2]
    l.extend(range(3, x_sqrt + 1, 2))
    for i in l:
        if x % i == 0:
            return False
    return True


def get_factor(x):
    """ 获取x的质数因子 """

    def factor_div(xx):
        """ 
        递归求取xx的质数因子
        若 xx = i * j,xx的质数因子是i和j的质数因子集合
        """
        if is_prime(xx):
            factor.add(xx)
            return factor
        xx_sqrt = int(math.sqrt(xx))
        for i in range(2, xx_sqrt + 1):
            if xx % i == 0:
                if i not in factor:
                    factor_div(i)
                j = int(xx / i)
                if j not in factor:
                    factor_div(j)
                break
        return factor

    factor = set()
    factor_div(x)
    if x in factor:
        factor.remove(x)
    return factor


print(max(get_factor(600851475143)))

这样改
然后解释下你的,如果按照你的
n = 20
那么余数为0
不会输出任何东西
你的程序有问题
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2025-2-11 15:59

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表