免费注册 查看新帖 |

Chinaunix

  平台 论坛 博客 文库
123
最近访问板块 发新帖
楼主: 大狗狗
打印 上一主题 下一主题

求个算法。 [复制链接]

论坛徽章:
0
21 [报告]
发表于 2009-01-07 13:46 |只看该作者
不知道perl运行要多久,我用C++,很快,大约1秒给出结果,和那位老兄结果一致。

#include <cstdlib>
#include <iostream>
#include <cmath>

using namespace std;

long F(unsigned long n)
{
  unsigned long sum = 0;
  for(unsigned long i=1; i<=n; ++i)
    sum += i;
  
  return sum;
}

int CountYueShu(unsigned long sum)
{
  unsigned long sq = sqrt(sum);
  int count = 0;
  
  for(unsigned long i=1; i<sq; ++i)
    if(sum % i == 0)
      count ++;
  
  count = count * 2;
  
  if(sum % sq == 0)
    count ++;
  
  return count;
}

int main(int argc, char *argv[])
{
    unsigned long sum=0;
    unsigned long n=1;
   
    while(1)
    {
      sum = F(n);
      if( CountYueShu(sum) >= 500 )
      {
        cout<< "The n is " <<n <<endl;
        cout<< "The F(n) is " << sum <<endl;
        break;
      }
      ++n;
    }
   
   
    system("PAUSE");
    return EXIT_SUCCESS;
}

论坛徽章:
0
22 [报告]
发表于 2009-01-07 13:52 |只看该作者
在我的机器上Perl跑了10多秒,而C++是1秒到2秒

论坛徽章:
0
23 [报告]
发表于 2009-01-07 15:48 |只看该作者
my $t = eval { map { my $t = $_; my $i = scalar (grep {$_}(map { $t % $_ ? 0 : 1 } (1 .. $t))); if ($i >= 500) {return $t; }} (map { ($_ + 1) * $_ / 2 } (1 .. 100000))};

map {if ((($_ + 1) * $_ / 2) == $t) {print $_-1; exit;}} (1 .. 100000);

---------------------

玩玩
没运行,只是测试过 50
您需要登录后才可以回帖 登录 | 注册

本版积分规则 发表回复

  

北京盛拓优讯信息技术有限公司. 版权所有 京ICP备16024965号-6 北京市公安局海淀分局网监中心备案编号:11010802020122 niuxiaotong@pcpop.com 17352615567
未成年举报专区
中国互联网协会会员  联系我们:huangweiwei@itpub.net
感谢所有关心和支持过ChinaUnix的朋友们 转载本站内容请注明原作者名及出处

清除 Cookies - ChinaUnix - Archiver - WAP - TOP