欧拉工程，第二十一题，计算10000以下亲和数的和。 - 程序语言 - 其它

版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

本帖产生 3 个程序强帖，点击这里进行查看

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huycwork

金虫 (著名写手)

程序强帖: 22
应助: 0 (幼儿园)
金币: 953
散金: 663
红花: 8
沙发: 13
帖子: 1080
在线: 264.1小时
虫号: 1257243
注册: 2011-04-06
专业: 金融学

★ ★ ★
dubo(金币+1): 多谢交流 2011-06-04 14:42:08
余泽成(金币+2, 程序强帖+1): 鼓励交流！ 2011-06-04 19:31:23

C++代码：

CODE:

#include 

enum {BUFSZ = 10000};

size_t eular21(){

        size_t buf[BUFSZ];

        memset(buf, 0, sizeof buf);

        for(size_t i = 1; i < BUFSZ; ++i){

                for(size_t j = i+i; j < BUFSZ; j+=i){

                        buf[j] += i;

                }

        }

        size_t d, s = 0;

        for(size_t i = 2; i < BUFSZ; ++i){

                d = buf[i];

                if(i == buf[d]){

                        if(i != d)

                                s += i;

                }

        }

        return s;

}

int main(){

        std::cout<
}

赞一下(2人)

回复此楼

漩涡的中心有一块空地，空空的。

2楼2011-06-02 22:34:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

Project Euler 50 欧拉工程 50 题已经有12人回复
fluent的欧拉-拉格朗日法求助已经有15人回复
欧拉气固两相流模拟已经有5人回复
欧拉模型已经有5人回复
Project Euler 48 欧拉工程 48 题已经有30人回复
Project Euler 45 欧拉工程 45 题已经有7人回复
Euler Project Q17. 欧拉工程第十七题已经有4人回复
Euler Project Q13 欧拉工程第十三题已经有20人回复
Euler Project Q12 欧拉工程第十二题已经有23人回复
Euler Project Q8. 欧拉工程第八题已经有4人回复
Euler Project Q7. 欧拉工程第七题已经有14人回复
【求助】欧拉模型中颗粒相granular temperature定义已经有6人回复
【求助】fluent模拟气固流化床采用欧拉模型并行计算出现问题已经有12人回复
【求助】fluent模拟两段流化床采用欧拉和DPM模型问题已经有11人回复
【求助】关于欧拉－拉格朗日方程（Euler-Lagrange equation）【已解决】已经有11人回复
【求助】有谁知道这个是什么个离散方式啊，傅里叶，欧拉....？已经有3人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

libralibra

至尊木虫 (著名写手)

骠骑将军

程序强帖: 40
应助: 817 (博后)
金币: 12914.1
红花: 64
帖子: 2238
在线: 287.3小时
虫号: 696514
注册: 2009-02-05
专业: 计算机软件

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
dubo(金币+1): 多谢交流 2011-06-04 14:43:53
余泽成(金币+2, 程序强帖+1): 鼓励交流！ 2011-06-04 19:31:45

看笨拙的matlab吧

CODE:

 function result = euler21()

tic;

result = [];

for i=1:10000

    if d(i)~=i && d(d(i))==i

        result = [result,i];

    end

end

result = sum(unique(result));

toc;

end

%% Let d(n) denote the sum of proper divisors of n (numbers less than n which divide evenly into n).

% For example, the proper divisors of 220 are 1, 2, 4, 5, 10, 11, 20, 22, 44, 55 and 110; therefore d(220) = 284.

% sub function to compute d(n), called by 21 and 23

function s = d(n)

s = 0;

for i=1:n-1

    if mod(n,i)==0

        s = s+i;

    end

end

end

赞一下(3人)

回复此楼

matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com

4楼2011-06-03 16:36:21

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

dubo

金虫 (著名写手)

程序强帖: 4
应助: 23 (小学生)
贵宾: 0.779
金币: 569.2
散金: 3220
红花: 31
沙发: 1
帖子: 1821
在线: 349.5小时
虫号: 559371
注册: 2008-05-17
性别: GG
专业: 高分子组装与超分子结构

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
ben_ladeng(金币+1): 奖励下下 2011-06-04 17:52:28
余泽成(金币+2, 程序强帖+1): 鼓励交流！ 2011-06-04 19:32:24

引用回帖:

Originally posted by wangww2011 at 2011-06-03 13:22:19:
结果

CODE:

31626

elapsed time=0.020000 seconds.

c代码
[code]
#include <stdio.h>
#include <math.h>
#include <time.h>

#define TIMERSTART clock_t start_time,stop ...

// play_time_calu.cpp : 定义控制台应用程序的入口点。
//

//#include "stdafx.h"

#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 sumdivisors(int n){
  int i,sum=1,sqrtn=sqrt(float(n));//sqrtn=sqrt(n) 不正确

  for(i=2;i if(n%i==0)sum+=i+n/i;
  }
  if(sqrtn*sqrtn==n)sum-=sqrtn;
  return sum;
}

int euler21(int n){
  int i,sum=0,tmp;
  for(i=3;i tmp=sumdivisors(i);
if(tmp!=i&&tmp    sum+=i;
}
  }

  return sum;
}

int main(void){
//int i；没必要
TIMERSTART;

printf("%d\n",euler21(10000));

TIMERSTOP;

  return 0;
}

赞一下(3人)

回复此楼

8楼2011-06-04 14:55:25

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 huycwork 的主题更新

返回列表