-1

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

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

在阅读了一些关于素数生成算法的信息后，我发现“阿特金斯筛子”是最快的素数生成算法。在本教程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;
}
4

1 回答 1

4

首先，阿特金筛子比埃拉托色尼筛子快是一个互联网神话。Atkin 筛在某些限制下更快，但 Eratosthenes 筛通常更快，并且是实际使用它来研究素数的数学家选择的筛。

其次，当筛子太大而无法放入内存时，您可以对筛子进行分段。有关详细信息，请参阅此答案。

于 2016-06-14T16:54:26.507 回答