Почему мне не удается Project Euler #10?
Вопрос в том, чтобы найти сумму всех простых чисел ниже 2 миллионов.
Я в значительной степени сделал вещь Sieve of Erastothenes, и программа ниже, кажется, работает для небольшого числа, то есть определите LIMIT, поскольку 10L выдает 17 как ответ.
Я отправил 1179908154 в качестве ответа, выданного следующей программой, и он был неверным.
Пожалуйста, помогите указать на проблему. Благодарю.
#include <stdio.h>
#define LIMIT 2000000L
int i[LIMIT];
int main()
{
unsigned long int n = 0, k, sum = 0L;
for(n = 0; n < LIMIT; n++)
i[n] = 1;
i[0] = 0;
i[1] = 0;
unsigned long int p = 2L;
while (p*p < LIMIT)
{
k = 2L;
while (p*k < LIMIT)
{
i[p*k] = 0;
k++;
}
p++;
}
for(n = 0; n < LIMIT; n++)
if (i[n] == 1)
{
sum += n;
}
printf("%lu\n",sum);
return 0;
}
3 ответа
Вы вычисляете простые числа правильно, но сумма слишком велика (более 2^32) и не помещается в 32-битную длину без знака. Вы можете использовать 64-битное число (long long
на некоторых компиляторах) чтобы это исправить.
Кажется, ваша логика верна, но вы путаетесь с типами данных и их диапазонами. Проверьте, работает ли это или нет:
#include <stdio.h>
#define LIMIT 2000000
int i[LIMIT];
int main()
{
long long int n = 0, k, sum = 0;
for(n = 0; n < LIMIT; n++)
i[n] = 1;
i[0] = 0;
i[1] = 0;
long long int p = 2;
while (p*p < LIMIT)
{
k = 2;
while (p*k <LIMIT)
{
i[p*k] = 0;
k++;
}
p++;
}
for(n = 0; n < LIMIT; n++)
if (i[n] == 1)
{
sum += n;
}
printf("%lld\n",sum);
return 0;
}
Output :142913828922
Вы также можете обнаружить, что вам также нужно использовать ключ компилятора -std=c99. Я сделал с gcc (GCC) 3.4.5 (mingw-vista special r3).
т.е.
gcc -Wall -std=c99 -o problem10 problem10.c