Question

自2天以来，我一直在努力解决SPOJ的PRIME1问题。我设法编写了一个对大数字有效的程序，但是我不知道为什么它会在开始范围内（通常是1-11）显示如此混乱的数字。我正在使用Eratosthenes的分段筛。 SPOJ link

问题出在这里

输入

输入以一行中的测试用例数t（t <= 10）开始。在接下来的t行中的每行中，都有两个数字m和n（1 <= m <= n <= 1000000000，n-m <= 100000），以空格隔开。

输出

对于每个测试用例，请打印所有素数p，以使m <= p <= n，每行一个数字，测试用例用空行分隔。

我的代码：

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

#define MAX         1000000000
#define MAX_SQRT    sqrt(MAX)

vector<bool> is_prime(MAX_SQRT, true);

void simpleSieve(int limit, vector<int>& primes)
{
    for(int p=2; p*p<=limit; p++) {
        if(is_prime[p] == true) {
            for(int i=p*p; i<=limit; i+=p)
                is_prime[i] = false;
        }
    }
    for(int i=2; i<=limit; i++) {
        if(is_prime[i]) primes.push_back(i);
    }
}

void segmentedSieve(int left, int right)
{
    int range = floor(sqrt(right)) + 1;
    vector<int> primes;
    simpleSieve(range, primes);

    int n = right - left + 1;

    bool sieve[n];
    memset(sieve, true, sizeof(sieve));

    for (int i = 0; i<primes.size(); i++) {

        int low = floor(left/primes[i]) * primes[i];
        if(low < left)
            low += primes[i];

        for(int a=low; a<=right; a+=primes[i]) {
            sieve[a - left] = false;
        }
    }
    for(int i=left; i<=right; i++)
        if(sieve[i - left])   cout << i << "\n";
}

int main()
{
    int t;
    cin >> t;
    while(t--) {
        int left, right;
        cin >> left >> right;
        segmentedSieve(left, right);
    }
    return 0;
}

我已经尝试对开头部分进行硬编码，但是根据输入的不同，它会有所不同。

Answer 1

首先打印大于质数的质数，这些质数是使用simpleSieve生成的。
另外，1不是素数，因此我们在打印其他素数时会跳过它：

for (int i = 0; i<primes.size(); i++) {
    if (primes[i] >= left)
        cout << primes[i] << "\n";
}

for(int i=left; i<=right; i++){
    if (i == 1){
        continue;
    }
    if(sieve[i - left])   cout << i << "\n";
}
cout << "\n";  //empty line before next testcase starts

素数生成器不适用于小数

1 个答案: