HDU 2012 素数判定(素数)

HDU 2012 素数判定(素数)

http://acm.hdu.edu.cn/showproblem.php?pid=2012

题意：水题一枚

       对于表达式n^2+n+41，当n在（x,y）范围内取整数值时（包括x,y）(-39<=x<y<=50)，判定该表达式的值是否都为素数。

分析：

       本题的数据范围很小，求出表达式在-39到50内的所有可能值可以得到下面的数：

1523 1447 13731301 1231 1163 1097 1033 971  911

853  797 743  691  641 593  547  503 461  421

383  347 313  281  251 223  197  173 151  131

113  97  83   71   61  53   47   43  41   41

43   47  53   61   71  83   97   113 131  151

173  197 223  251  281 313  347  383 421  461

503  547 593  641  691 743  797  853 911  971

1033 1097 11631231 1301 1373 1447 1523 1601 1681

1763 1847 19332021 2111 2203 2297 2393 2491 2591

最大的数为2591，所以我们只要求出3000范围内所有素数，然后每次求出表达式范围值，判断都是否素数即可。

       当然本题更快的方法是预处理所有可能的数，看上面矩阵中哪些是素数，哪些不是素数，然后直接对于范围x，y输出结果。

AC代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=3000;

//生成素数表
bool flag[maxn+5];//flag[i]==true表i为素数
int prime[maxn+5];
int get_prime()
{
    for(int i=2;i<=maxn;i++)
    {
        if(prime[i]==0)
        {
            flag[i]=true;
            prime[++prime[0]]=i;
        }
        for(int j=1;j<=prime[0] && prime[j]<=maxn/i;j++)
        {
            prime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
    return prime[0];
}

//计算公式的值
int compute(int x)
{
    return x*x+x+41;
}

//判断范围内是否都为素数
bool check(int a,int b)
{
    for(int i=a;i<=b;i++)
        if(!flag[compute(i)]) return false;
    return true;
}

int main()
{
    get_prime();
    int a,b;
    while(scanf("%d%d",&a,&b)==2)
    {
        if(a==0&&b==0) break;
        printf("%s\n",check(a,b)?"OK":"Sorry");
    }
    return 0;
}