HDU - 2012 素数打表

最新推荐文章于 2024-01-15 19:12:17 发布

原创最新推荐文章于 2024-01-15 19:12:17 发布 · 198 阅读

0 ·

CC 4.0 BY-SA版权

本文介绍了一种快速判断素数并计算指定范围内素数个数的方法，通过预先生成素数表和素数个数求和表，有效避免了在大规模数据输入下普通素数判断方法的超时问题。

注意题目给定时间限制：1000ms

而且数据的输入规模为多组输入（很庞大）

因此使用普通的素数判断的方法做题妥妥的超时，所以我们使用打表的方法。

以下给出本题代码：

#include <iostream>
#include <cstring>
#include <ctime>

using namespace std;

const int maxn = 3005;

bool prime[maxn];
int num[90];  //n=index-39

int caculate(int n)   //公式写成函数形式方便书写
{
    return n*n+n+41;
}

void get_prime()     //素数打表
{
    memset(prime, true, sizeof(prime));
    prime[0] = false;
    prime[1] = false;
    for (int i = 2; i <= 1300; i++) {
        for (int j = 2; i*j < 2600; j++) {
            prime[i*j] = false;
        }
    }
}

void get_num()      //素数个数求和表
{
    memset(num, 0, sizeof(num));
    if (prime[caculate(-39)] == true) num[0] = 1;
    for (int i = 1; i < 90; i++) {
        if (prime[caculate(i-39)] == true) num[i] = num[i-1] + 1;
        else num[i] = num[i-1];
    }
}

int main()
{
    ios::sync_with_stdio(false);     //加快输入速度
    int x, y;

    get_prime();
    get_num();

    while (cin >> x >> y) {
        if (!x && !y) break;
        int prime_num = prime[caculate(x)]==true ? num[y+39]-num[x+39]+1 : num[y+39]-num[x+39];  //判断

        if (y-x+1 == prime_num) cout << "OK" << endl;
        else cout << "Sorry" << endl;
    }

    return 0;
}