应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第三题：寻找600851475143的最大质因子
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
52/1返回列表
查看: 3649  |  回复: 18
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 7 个 程序强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

holmescn

金虫 (正式写手)

[交流] Euler 工程 第三题：寻找600851475143的最大质因子已有7人参与

昨天没有放出第三题，今天赶早补上。
前两个题目都比较简单了，只要会基本的数学和编程语言，就可以完成。
第三题就有点意思了。

第三题：寻找一个合数的最大质因数

对一个数（非质数）进行因数分解，比如13195＝5x7x13x29。最大的质因数是29.
那么 600851475143 怎么分解呢？最大的质因数又是多少？

[ Last edited by holmescn on 2011-5-12 at 15:06 ]
回复此楼

» 猜你喜欢
1楼2011-05-10 08:51:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feign_te

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
如果已经有一个合适大小的素数列，肯定是最快的。
单独为一个问题去求素数列不太经济。
方法和前面一样，n除遍奇数找出第一个质因数p1之后，用n/p1找出下一个质因数p2，可代码写出来感觉差别好大啊。
CODE:
clear
tic
n = 600851475143;
prime=zeros(fix(log2(p)),1);
k=0;
m=n;
p=fix(n^.5);
r=3;
    for i =r:2:p
        if mod(m,i)==0
            k=k+1;
            r=i;         
            prime(k)=i;
            m=n/i;
            p=fix(m^.5);
            continue
        end     
    end
prime=prime(find(prime~=0))
toc

运行结果
CODE:
prime =
          71
         839
        1471
        6857
Elapsed time is 0.021253 seconds.
赞一下 回复此楼
19楼2013-12-10 06:53:21
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

微尘、梦想

木虫 (知名作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xzhdty(金币+1): 谢谢微尘、梦想 斑斑 2011-05-10 15:10:46
余泽成(程序强帖+1): 2011-05-12 19:07:28
CODE:
#include "stdio.h"
void main(void)
{
    int i,a;
    printf("请输入一个整数：" );
    scanf("%d",&a);

    for(i=2;a!=1;i++)
        if(a%i==0)
        {
            a/=i;
            printf("%d\t",i);
            i--;
        }
        printf("\n" );
}

由于32位内存的限制，无法求出太大的数！

[ Last edited by 微尘、梦想 on 2011-5-12 at 16:19 ]
赞一下(3人) 回复此楼
任风云变幻，我笑对人生！
2楼2011-05-10 13:19:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
xzhdty(金币+2): 呵呵谢谢应助 2011-05-10 15:11:23
余泽成(程序强帖+1): 2011-05-12 19:07:37
此题很变态,matlab运行55s,
CODE:
function result = euler3()
tic;
result = 0;
n = 600851475143;
for i=3:sqrt(n)
    if isprime(i)==1 && mod(n,i)==0
        result = i;
    end
end
toc;
end

答案
CODE:
Elapsed time is 55.277722 seconds.
ans =
        6857

[ Last edited by libralibra on 2011-5-10 at 16:46 ]
赞一下(3人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
3楼2011-05-10 13:58:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★
余泽成(金币+2): 辛苦了！ 2011-05-11 23:01:56
其实有个作弊的解法哈哈，用Mathematica直接

FactorInteger[600851475143]

当然3楼的结果是对的。
不过，好像因为是线性查找，效率才不高。还有，干什么不从大到小找呢？那样快很快的。
赞一下(2人) 回复此楼
4楼2011-05-10 14:40:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭