24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3386  |  回复: 23
本帖产生 7 个 程序强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

libralibra

至尊木虫 (著名写手)

骠骑将军

[交流] Euler Project Q12 欧拉工程第十二题已有9人参与

Question 12:
The sequence of triangle numbers is generated by adding the natural numbers. So the 7th triangle number would be 1 + 2 + 3 + 4 + 5 + 6 + 7 = 28. The first ten terms would be:

1, 3, 6, 10, 15, 21, 28, 36, 45, 55, ...

Let us list the factors of the first seven triangle numbers:

     1: 1
     3: 1,3
     6: 1,2,3,6
    10: 1,2,5,10
    15: 1,3,5,15
    21: 1,3,7,21
    28: 1,2,4,7,14,28

We can see that 28 is the first triangle number to have over five divisors.

What is the value of the first triangle number to have over five hundred divisors?

翻译:

自然数求和可生成三角数列.第七个三角数是1 + 2 + 3 + 4 + 5 + 6 + 7 = 28.前10个三角数列元素是:
1, 3, 6, 10, 15, 21, 28, 36, 45, 55, ...

下面列出前7个三角数及其整除数:
     1: 1
     3: 1,3
     6: 1,2,3,6
    10: 1,2,5,10
    15: 1,3,5,15
    21: 1,3,7,21
    28: 1,2,4,7,14,28
可以看出,28是第一个有超过5个整除数的三角数.

那么,第一个有超过500个整除数的三角数是多少?

[ Last edited by libralibra on 2011-5-17 at 05:41 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)

引用回帖:
Originally posted by libralibra at 2011-05-17 05:38:26:
对此题我相当无语,难度在20题以内,我也不知道具体程序跑了多久,反正我吃了个泡面回来,打印出结果了,应该超过10分钟了

c的,matlab对此表示鸭梨很大

[code] #include <stdio.h>
#include <stdlib.h ...

你又遍历!
漩涡的中心有一块空地,空空的。
3楼2011-05-17 11:10:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
余泽成(金币+3, 程序强帖+1): 好详细的分析! 2011-05-17 22:24:53
想到一个依赖组合计数的算法:
首先分析特殊情况,对于一个a1*a2*a3*a4*...*an=A酱紫的数,假定a1 如,2*3*5*7=210的因数数有:
2, 3, 5, 7
6, 10, 14,15, 21,35
30, 70,105,42
210
麻烦的问题是不完全是素数的情况,对于任意C=A*B,令A和B为素数积,当且仅当AB没有交集的时候C才是素数积,否则就带来可可怕的乘方,乘方会产生很多相同的因数。为了解决这个问题,需要分析素数的乘方能产生多少个相同解:
首先分析下素数乘方可能产生的相同解个数,假定某数可以分解出n个相同的素数因子:a1,a2,...an,这些因子的恼人之处就在于,不同的乘数可能产生相同的乘积,原本这些乘积是独立计数的,最终却需要剔除。一个简单的解决办法是,将n个数叠起来,使它们像这样:a1, a1^2, a1^3...a1^n,对这样的数列再次应用素数组合,但分析的时候要独立分析:这些数中只取一个的时候,产生独立解,其余的情况,都不能再计数,理由是,数列中的任意数的组合要么相同,要么不可能产生。因此,对于n个素数相同的情形,只能将其算一个素数的组合位置,但是要记入n种状况。最终的组合要借助计算机来完成,举例如下:
8*6 = 2*2*2*2*3,分离素数乘方表和素数表,得
乘方表:2,4,8,16,素数表2,3
在2的位置需要替换4种情况,当C21时,2本来只计数1,现在需要计数4,C22时2*3本来只有一个解,现在有4个解。对于2和3而言,原本的解像这样:
2, 3, 2*3
现在的解要将4个值替换到2的位置:
2,4,8,16,3,2*3,4*3,8*3,16*3
在更多情况下和这个类似,如果是24*6 = 2*2*2*3*2*3,分离出来的素数表仅有2,3,乘方表则有两个:2,4,8,16和3,9
原来的解是2,3,2*3,分两遍替换进去就成了这样:
替换2变成:2,4,8,16,3,2*3,4*3,8*3,16*3
再替换掉3变成:2,4,8,16,3,9,2*3,2*9,4*3,4*9,8*3,8*9,16*3,16*9
可以看到,由于无法知道哪个数是需要替换的,所以也无法直接构建公式,光组合算法还不成,得在组合算法之上再递归计数。
生成因数表和递归计数都不是啥问题,问题是组合算法,排列组合的算法的复杂度很高啊…
漩涡的中心有一块空地,空空的。
5楼2011-05-17 13:41:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
余泽成(金币+2): 鼓励交流! 2011-05-17 22:26:43
引用回帖:
Originally posted by sudo at 2011-05-17 12:54:30:
抛砖引玉一个思路:

[code]
1. 建立素数表prime[N],N取一个比较大的值
2. 对x=(1+i)*i/2进行质因数分解:
    x=k_1^m_1 * k_2^m_2 ... * k_n^m_n
    其中^表示指数,k_1 ... k_n表示素数
3. 计算x因子 ...

素数表倒是不需要估计,俺给出的这个素数算法可以动态定制素数表,在参数里面传递函数指针进去就可以随意处理每个素数。不过这个算法有点小bug,会漏掉一个数的测试,修改offset为3的倍数,BUFSZ修改为6的倍数即可解决这个问题。http://muchong.com/bbs/viewthread.php?tid=3192387&fpage=2#pid1151737

[ Last edited by huycwork on 2011-5-17 at 14:31 ]
漩涡的中心有一块空地,空空的。
6楼2011-05-17 13:54:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
余泽成(金币+1): 鼓励交流! 2011-05-17 22:29:39
引用回帖:
Originally posted by sudo at 2011-05-17 14:18:50:
嗯,有空我拜读一下~

话说你5楼的算法复杂化了...

像4楼的算法就直接应用排列组合里面的乘法定理就出来了,不需要那么麻烦的

先以你举的210=2*3*5*7为例:

因数一共有16个:
1,[ ...

不错不错,还得打基础呃~~~~
这样一来,复杂度就主要看如何分离因数表了。
结果又回到了素数算法呃……
漩涡的中心有一块空地,空空的。
9楼2011-05-17 14:40:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
余泽成(金币+1): 鼓励参与交流! 2011-05-19 19:16:36
引用回帖:
Originally posted by holmescn at 2011-05-18 17:52:06:
虽然有了因数分解这个方法,可以快速的知道一个三角数有多少个整除数,但好像还只胡一个一个算,一个一个分解。这太慢了。

显然第n个三角数就是n(n+1)/2。令这个数为a则可以解得n=(sqrt(8a+1)-1)/2

这就要求 ...

辗转相除法可以计算gcd,计算完之后两边再单独分解即可。
漩涡的中心有一块空地,空空的。
18楼2011-05-18 19:48:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 libralibra 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见