蟒

Question

完美数字是其适当除数之和的数字   完全等于数字。例如，正确的总和   28的除数是1 + 2 + 4 + 7 + 14 = 28，这意味着28   是一个完美的数字。

如果正确的除数之和为，则数字n称为不足   小于n，如果此总和超过n，则称为丰富。

由于12是最小的数，1 + 2 + 3 + 4 + 6 = 16，所以   最小的数字，可以写成两个丰富数字的总和   通过数学分析，可以显示所有整数   大于28123可写为两个数字的总和。   但是，通过分析不能进一步降低该上限   即使知道最大数量是不可能的   表示为两个有限数的总和小于此限制。

求出不能写成的所有正整数的总和   两个数字的总和。

问题的答案是4179871，我的程序显示3797954。

首先，我创建了一个函数来填充数字丰富的[]，其中包含28124以下的所有数字。这非常好用，因为我搜索了大量的数字，它们与我的数组完全匹配。

其次，我有另一个数字为1-28123的数组，我假设它们中的所有数字“不能写成两个数字的总和。”这些都写入数组hold []。

最后，我将可以写成两个数字之和的数字除去，通过在[[]中添加所有数字，并使用[]中的所有数字，并将hold []的值设置为0。（保持[丰富[0到n] +丰富[0到n]] = 0）在hold []中添加所有剩余的数字，我只得到3797954

我知道这个程序效率不高，因为它为所有数量增加了所有数量，但它应该可以正常工作。怎么了？

#include <iostream>
#include <cmath>

using namespace std;

int hold[28124];
int abundant[7000]; //hold abundant numbers, there are only 6919 abundant numbers below 28123

bool abundance(int x){  //returns true if x is abundant
    int counter = 1;    //holds "proper divisors" of numbers, default by 1 because every int is divisible by 1
    for (int i = 2; i < sqrt(x); i++){   //finds all divisors 2 - sqrt(n)
        if (x % i == 0){
            counter += i;
            counter += x / i;
        }
    }
    int y = sqrt(x);   
    if (x % y == 0){   //adds sqrt(n) if its modulus to n is 0
            counter += sqrt(x);
        }
    if (counter > x){
        return true;
    } else {
        return false;
    }
}

int main()
{
    int counter = 0;
    for (int i = 0; i < 28124; i++){ //assumes every number cannot be written as the sum of two abundant numbers,
        hold[i] = i;                 //goes up to 28123 because "it can be shown that all integers greater
    }                                //than 28123 can be written as the sum of two abundant numbers." - project euler
    for (int j = 10; j < 28124; j++){ 
        if (abundance(j) == true){  //copies all abundant numbers up to 28123 to abundant[]
            abundant[counter] = j;
            counter++;
        }
    }

    for (int m = 0; m < counter; m++){  //adds all numbers in abundant[], with all numbers in abundant[]
        for (int n = 0; n < counter; n++){
            if (abundant[m]+abundant[n] < 28124){
                hold[abundant[m]+abundant[n]] = 0; //sum of the abundant numbers in hold[] is set to 0
            } //hold[] now holds all numbers that cannot be written as the sum of 2 abundant numbers
        }
    }
    int counter2 = 0;
    for (int x = 0; x < 28124; x++){
        counter2 += hold[x];
    }
    cout << counter2 << endl;

}

Answer 1

问题出在你的abundance函数中，特别是这部分：

int y = sqrt(x);   
if (x % y == 0){   //adds sqrt(n) if its modulus to n is 0
    counter += sqrt(x);
}

x % (int)sqrt(x) == 0并不意味着sqrt(x)是一个整数。一个简单的反例是2. sqrt(2)约为1.414，或仅1为整数。但是2 % 1 == 0，即使它不是平方根。

因此，要修复代码，请将该部分更改为：

int y = sqrt(x);   
if (y * y == x){   //adds sqrt(n) if sqrt(n) is an integer
    counter += y;
}

你会得到正确的结果。

Answer 2

您可以查看显示的代码here - 只是因为它得到了正确的答案。如果不复制那里所做的事情，您可以确认问题是在您的丰富数字生成器中还是在其他部分中。它可能会帮助你弄清楚你哪里出错了。

Answer 3

my_set = set([i for i in range(1, 28123)])
print ('original sum',sum(my_set))
my_list = list(my_set)
my_l1 =[]
for k in range(12, int(my_list[-1]/2+1)):
    a = 0
    for j in range(1, int(k/2+1)):
        if k%j == 0:
            a += j

    if a > k:
        my_l1.append(k)
        my_set.remove(k*2)
#Calculating the sum of all the numbers which can be written as the sum of two abundant numbers   
l = 0
my_l2 = set([])
for d in my_l1:
    l += 1
    k = l
    while k < len(my_l1):
        a_s = d + my_l1[k]
        my_l2.add(a_s)
        k += 1
my_set.difference_update(my_l2)
print ('sum of all abbundant numbers:',sum(my_set))

这是我的问题23项目Euler的代码，这有什么问题？我现在不关心运行时，我只想要正确答案。

项目Euler＃23，在程序中找不到问题

3 个答案:

蟒