Project Euler - Problem 21

亲和数求和算法

最新推荐文章于 2025-12-09 18:25:52 发布

原创最新推荐文章于 2025-12-09 18:25:52 发布 · 395 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#iterator #numbers #include #list #ios #测试

本文介绍了一种用于寻找并计算10000以内所有亲和数之和的C++算法实现。通过构建素数列表辅助判断，利用list存储因数进行计算，最终输出所有符合条件的亲和数及其总和。

#include <iostream> #include <fstream> #include <string> #include <time.h> #include <list> using namespace std; //计算一个数的因数并放入list中的因数和，并且把list内的因数清空 int divisorsum (int n); //判断一个数是否是亲和数 bool IFAmicableNumber (list<int>′ ,int n); int main() { double start = clock(),end(0); int dsum = 0;//某个数因数之和 int SumOfAll = 0;//所有10000以下亲和数之和 bool cheak = false; //判断是否为亲和数 //收集所有的素数 int N = 3; bool cheakN = true;//判断是否为素数 int a[] = {2,3}; list<int> prime(a,a+2); //cout.setf(ios::fixed); while(N != 10001) { for(list<int>::iterator iter= prime.begin();iter!=prime.end();++iter) { if(N % (*iter) == 0) { N+=2; cheakN = false; break; } } if(cheakN) { if(N < 3000) prime.insert(prime.end(),N); N+=2; } cheakN = true; } for(int n = 1; n<10001 ; n++) //被测试数 { if(IFAmicableNumber (prime ,n)) continue; dsum = divisorsum(n); if(n!=dsum && n == divisorsum(dsum)) { cout << "amicable numbers is " << n << " and " << dsum << endl; SumOfAll += n; prime.push_front(dsum); } } cout << endl << "the sum of all the amicable numbers under 10000 is " << SumOfAll << endl; end = clock(); double result = (end - start)/1000; cout.setf(ios::fixed); cout <<endl<< "3.00GHz's run time = " << result << "s."<<endl; return 0; } //计算一个数的因数并放入list int divisorsum (int n) { list<int> divisor; for(int d = 1; d<n ; d++) //因数判断 { if (0 == n % d) divisor.push_back(d); } int sum = 0; for(list<int>::iterator iter = divisor.begin(); iter!=divisor.end(); iter++) { sum += *iter; } divisor.clear(); return sum; } //判断一个数是否是亲和数 bool IFAmicableNumber (list<int>′ ,int n) { for(list<int>::iterator iter= prime.begin();iter!=prime.end();++iter) { if(*iter == n) return true; } return false; }