Question

我解决了问题10 in project euler（但我的解决方案需要2000秒），所以我试图找到更快的解决方案。我发现了这个（执行时间0.1秒），但我不明白这个解决方案，我需要帮助。

#include <stdio.h>
#include <string.h>
#include <math.h>
using namespace std;
#define HIGHEST 2000000

char Prime[HIGHEST / 2];

int main(int argc, char** argv)
{
   unsigned int i, j;
   unsigned long long sum = 2;
   unsigned int total = 0;

   /* Set entire array to true (prime) */
   memset(Prime, 1, sizeof(char)*HIGHEST / 2);
   /* except for 1 */
   Prime[0] = 0;

   for(i = 3; i < HIGHEST; i += 2) {
      if(Prime[i / 2] == 1) {
         sum += (unsigned long long)i;
         total++;
      }
      for(j = (i+i+i)/2; j < HIGHEST/2; j += i) {
         Prime[j] = 0;
      }
   }
   printf("Sum: %llu (%d prime numbers)\n", sum, total);
   return 0;
}

Answer 1

解决方案是Eratosthenes筛子的一个有点奇怪的版本。看起来作家试图通过循环计数器i做一些偷偷摸摸的技巧来过早地进行优化。

我建议您在纸上完成循环的前几次迭代，以了解代码的作用。

在projecteuler中混淆问题10的解决方案

1 个答案: