Question

我在spoj.com上遇到了这个问题

http://www.spoj.com/problems/PRIME1/

在阅读了关于素数生成算法的一些内容之后，我发现了＃34;阿特金斯的筛选＆＃34;是最快的素数生成算法。我已经在本教程http://www.geeksforgeeks.org/sieve-of-atkin/

的帮助下部分实现了算法

当我运行代码时，它给了我分段错误。调试后我才知道这段代码不是最优的，它使用了更多的内存，存储空间。有人能告诉我如何优化程序。这是我的代码片段。是的，它是不完整的。

#include <iostream>
#include <bits/stdc++.h>
using namespace std;

int sieveOfAtkin(long int limit)
{
    // if limit is either 2 or 3, print them
    if (limit >= 2)
    {
        cout << 2 << endl;
    }

    if (limit >= 3)
    {
        cout << 3 << endl;
    }


    bool sieve[limit];
    for (long int i=0; i<=limit; i++)
        sieve[i] = false;

    for (long int x=1; x<limit; x++)
    {
        for (long int y=1; y<limit; y++)
        {

            long int n = (4*x*x) + (y*y);
            if (n<=limit && (n%12==1 || n%12==5))
                sieve[n] ^= true;

            n = (3*x*x) + (y*y);
            if (n<=limit && (n%12==7))
                sieve[n] ^= true;

            n = (3*x*x) - (y*y);
            if (x>y && n<=limit && (n%12==11))
                sieve[n] ^= true;
        }
    }

    for (long int r=5; r*r<limit; r++)
    {
        if(sieve[r])
        {
            for(long int i=r*r; i<limit; i+=r*r)
                sieve[r] = false;
        }
    }

    for (long int i=5; i<limit; i++)
        if(sieve[i])
            cout << i << endl; 

    return 0;
}

int main(void)
{
    long int x;
    cin >> x;
    sieveOfAtkin(x);
    return 0;
}

Answer 1

首先，互联网神话是阿特金的筛子比Eratosthenes筛子更快。在某些限制条件下，阿特金的筛子速度更快，但是Eratosthenes的筛子总体上更快，并且是实际用它来研究质数的数学家选择的筛子。

其次，当筛子太大而不适合记忆时，你可以分割筛子。有关详细信息，请参阅this answer。

value需要1000000000个字节，超过max-value-size

1 个答案: