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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第三题：寻找600851475143的最大质因子
191/2返回列表12下一页
查看: 4342  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

holmescn

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

FactorInteger[600851475143]

当然3楼的结果是对的。
不过，好像因为是线性查找，效率才不高。还有，干什么不从大到小找呢？那样快很快的。
赞一下(2人) 回复此楼
4楼2011-05-10 14:40:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

微尘、梦想

木虫 (知名作家)
★ ★
小木虫(金币+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的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 谢谢参与！ 2011-05-10 19:08:53
引用回帖:
Originally posted by holmescn at 2011-05-10 14:40:09:
其实有个作弊的解法哈哈，用Mathematica直接

FactorInteger[600851475143]

当然3楼的结果是对的。
不过，好像因为是线性查找，效率才不高。还有，干什么不从大到小找呢？那样快很快的。

CODE:
sqrt(600851475143)
ans =
          775146.099224527

结果是6857,从小到大快,只检测了6800多个数
从大到小,要检测77万多个数,时间就长了.呵呵
赞一下(2人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
5楼2011-05-10 16:50:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
微尘、梦想(金币+2): 谢谢参与！ 2011-05-10 19:08:27
引用回帖:
Originally posted by libralibra at 2011-05-10 16:50:34:
CODE:
sqrt(600851475143)
ans =
          775146.099224527

结果是6857,从小到大快,只检测了6800多个数
从大到小,要检测77万多个数,时间就长了.呵呵

如果从小到大检测的话

意味着需要验证求出的质因数是否是最大（不然怎么知道是最大质因数而没有更大的呢？），如果这么做会浪费更多时间

不如从大到小判定了...


PS：
路过...话说看到标题里面的分类是【其他】....我还以为是不熟悉的领域呢...原来是编程题啊....=,=|||会不会也有别人有同样的感觉然后就没打开帖子看看？

[ Last edited by sudo on 2011-5-10 at 18:50 ]
赞一下(2人) 回复此楼
6楼2011-05-10 18:47:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+1): 鼓励交流！ 2011-05-11 23:02:19
sudo所言极是,呵呵,
从大到小,找到直接就break,的确循环次数少
赞一下(2人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
7楼2011-05-10 21:09:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+2): 鼓励交流！ 2011-05-12 19:06:10
咳，其实这个问题相当有现实意义了...

看雪的密码学小组一直在研究这个....

目前貌似100位以下的整数的最快方法是二次筛法（咳在一本数论书里面说是115位），然后以上的目前最快的方法是数域筛法

谁来挑战一下二次筛法？
赞一下(2人) 回复此楼
8楼2011-05-12 08:35:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★
余泽成(金币+2, 程序强帖+1): 辛苦了！ 2011-05-12 19:08:14
写了一个python版的，不知道算不算筛选的。
CODE:
import timeit
from math import sqrt

def euler3():
    n = 600851475143

    primes = range(2, int(sqrt(n)))

    while n > 1:
        prime = primes[0]
        if n % prime == 0:
            print prime
            n = n / prime
        primes = [x for x in primes if x % prime != 0]


t = timeit.Timer("euler3.euler3()", "import euler3")
print t.timeit(1)

在我的电脑上用时不到20秒
不过，应该还能优化。因为用Mathematica不到1秒

PS:这个sqrt(n)假设好像有问题啊。比如本题的数分解为71, 839, 1471, 6857，如果是后两个数的积是10086647，开根号是3175.9，显然找不到后面一个数了啊。

[ Last edited by holmescn on 2011-5-12 at 11:15 ]
赞一下(1人) 回复此楼
9楼2011-05-12 10:14:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★
余泽成(金币+2): 呵呵！ 2011-05-12 19:08:26
非常的悲催啊，同样的算法，在matlab里只用不到7秒
CODE:
n = 600851475143;
tic;
primes = 2:round(sqrt(n));
while n > 1
    prime = primes(1);
    if mod(n, prime) == 0
        disp(num2str(prime));
        n = n / prime;
    end
    primes(find(mod(primes, prime)==0)) = [];
end
toc
赞一下(1人) 回复此楼
10楼2011-05-12 10:31:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 holmescn 的主题更新
191/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊 +7 Tide man 2026-08-10 8/400 2026-08-11 15:49 by 铁帽子农民
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] FileCode能看出啥？ +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +10 yufeiwaner 2026-08-09 11/550 2026-08-10 18:47 by yufeiwaner
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读？ +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭