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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler Project Q12 欧拉工程第十二题
52/1返回列表
查看: 3391  |  回复: 23
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 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
1楼2011-05-17 05:30:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangww2011

木虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+2, 程序强帖+1): 谢谢参与交流！ 2011-05-17 22:30:26
楼上的分析的很多阿
还是直接写代码吧
最简单就是遍历了，其实也不是很慢，2秒多一点
优化一下也没有快多少，0.1s
先看运行结果
CODE:
triangle_num=76576500
slow version, elapsed time=2.060000 seconds.
triangle_num=76576500
normal version, elapsed time=0.090000 seconds.

代码：
CODE:
#include
#include
#include

#define TIMERSTART clock_t start_time,stop_time;double elapsed_time;start_time = clock();
#define TIMERSTOP stop_time = clock();elapsed_time=(double)(stop_time-start_time)/CLOCKS_PER_SEC;printf("elapsed time=%f seconds.\n",elapsed_time);

int euler12(int a, long n){
  static int last_a=0;
  static int count=1;
  int tmp_count=0;
  int i=0;
  
  for(i=a;i     if(n%i==0){
      if(i>last_a){
         last_a=i;
         tmp_count=count;
         count=2;
         return tmp_count*euler12(i,n/i);
       }
      
       count++;      
       return euler12(i,n/i);
    }
  }


  tmp_count=count;
  count=1;
  last_a=0;

  if (i*i-n==0){
    return 3*tmp_count;
  } else {
    return 2*tmp_count;
  }
}



int slow(int n){
  int i=0;
  int count=0;
  for(i=2;i<=sqrt(n);i++){
    if(n%i==0){
      count++;
     }
  }

  if (i*i-n==0) {
    return 2*count+1;
  }

  return 2*(count+1);
}


int main(void){
  long triangle_num=1;
  int i=2;
  int divisor_num=500;
  

  TIMERSTART;

  do {
    triangle_num+=i;
    i+=1;
  }while(slow(triangle_num)
  printf("triangle_num=%ld\n slow version, ",triangle_num);
  TIMERSTOP;



  start_time = clock();  
  triangle_num=1;
  i=2;

  do {
    triangle_num+=i;
    i+=1;
  }while(euler12(2,triangle_num)
  printf("triangle_num=%ld\n normal version, ",triangle_num);
  TIMERSTOP;

  return 0;
}
赞一下(2人) 回复此楼
11楼2011-05-17 16:04:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答

libralibra

至尊木虫 (著名写手)

骠骑将军
★ ★ ★
ben_ladeng(金币+3, 程序强帖+1): 欢迎继续 2011-05-17 06:43:27
对此题我相当无语,难度在20题以内,我也不知道具体程序跑了多久,反正我吃了个泡面回来,打印出结果了,应该超过10分钟了

c的,matlab对此表示鸭梨很大
CODE:
#include
#include

int main(int args, char* argv[])
{
        long result=1,trinum=0,curnum=1;
        int n=0,i=0,big=0;
        int stop = 500;                //停止要求的除数个数,改成5可以测试28
        while(n         {
                n = 0;
                trinum = 0;
                // 计算当前三角数
                for(i=1;i<=curnum;i++)
                        trinum += i;

                // 从1开始除,计算除数个数
                for(i=1;i<=trinum;i++)
                {
                        if(trinum%i==0) // 如果整除,n+1
                                n += 1;
                }
               
                if(n>=stop) //如果超过stop个,保存结果跳出循环,这里是500个
                {
                        result = trinum;
                        break;
                }

                //如果没超过stop个,记录当前最大的个数并输出
                //这个if可以不要,因为我运行以为电脑死机了
                //所以加了这个判断,找到更接近stop的除数个数时输出,让我有点盼头
                if(n>big)  
                {
                        big = n;
                        printf("Current number: %i [%d]\n",trinum,big);
                }

                //自加自然数增加1,trinum = sum(1:curnum)
                curnum++;
        }
       
        // 打印结果
        printf("%i\n",result);
       
        system("PAUSE");
        return 0;
        }

结果
CODE:
% ans =
%       76576500

[ Last edited by libralibra on 2011-5-17 at 18:14 ]
赞一下(1人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
2楼2011-05-17 05:38:26
 已阅   回复此楼   关注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的回帖

sudo

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+3, 程序强帖+1): 谢谢参与交流！ 2011-05-17 22:10:21
抛砖引玉一个思路：
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因子总数：
    divisors = (m_1 + 1) * (m_2 + 1) * ... * (m_n + 1)
4. 判断divisors是否超过500，如果没有，跳到第2步，尝试下一个x；如果已经超过500，则输出x

有素数表之后，这个算法应该比较快，因为质因数分解的时候，只要从小到大判断对素数的整除性，然后顺便统计质因数的指数，同时如果能整除，试除用的x可以缩小为
CODE:
x/=prime[i]

再继续作测试，这样，算法很快就收敛了

问题在于事先对prime[N]中的N的估计怎样才合理，恩...

[ Last edited by sudo on 2011-5-17 at 12:56 ]
赞一下(2人) 回复此楼
4楼2011-05-17 12:54:30
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 24 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭