POJ 2635 The Embarrassed Cryptographer(素数线性筛法)

Description
给出一个由两个素数相乘得到的大数s和一个整数l,如果两个素数中的较小值小于l则输出BAD和这个素数,否则输出GOOD
Input
多组用例,每组用例包括两个整数s和l,以文件尾结束输入
(4<=s<=10^100,2<=L<=10^6)
Output
对于每组用例,如果组成s的两个素数中的较小值小于l则输出BAD和这个素数,否则输出GOOD
Sample Input
143 10
143 20
667 20
667 30
2573 30
2573 40
0 0
Sample Output
GOOD
BAD 11
GOOD
BAD 23
GOOD
BAD 31
Solution
首先预处理出10^6以内的素数,这个用普通筛法会爆,所以要用素数线性筛法,筛完之后枚举l以内的素数,判断其是否能整除s,这个要用到大数取模,即先将s转化为千进制(万进制需用long long,因为有乘法会爆int),然后每次模拟取模过程判断s是否能整除当前素数,如果找不到满足条件的素数则是GOOD,找到了就是BAD,然后输出这个素数即可
Code

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
#define maxn 1111111
int mark[maxn],prime[maxn],res;//prime存储素数,从0开始,mark[i]表示i的最小素因数 
void get_prime()//素数线性筛法,得到n以内所有素数共res个 
{
    memset(mark,0,sizeof(mark));
    res=0;
    for(int i=2;i<maxn;i++)
    {
        if(!mark[i]) 
            mark[i]=prime[res++]=i;
        for(int j=0;j<res&&prime[j]*i<maxn;j++)
        {
            mark[i*prime[j]]=prime[j];
            if(i%prime[j]==0) 
                break;
        }
    }
} 
int main()
{
    get_prime();
    char s[111];
    int k[33],l,n;
    while(~scanf("%s%d",s,&l),s[0]!='0'&&l)
    {
        n=0;
        int i=strlen(s)-1,j,temp,mul,flag=1;
        while(i>=0)
        {
            j=0,temp=0,mul=1;
            while(j<3&&i>=0)
            {
                temp=temp+(s[i]-'0')*mul;
                i--,j++,mul*=10;
            }
            k[n++]=temp;
        }
        i=0;
        while(prime[i]<l)
        {
            temp=0;
            for(j=n-1;j>=0;j--)
                temp=(temp*1000+k[j])%prime[i];
            if(!temp)
            {
                flag=0;
                printf("BAD %d\n",prime[i]);
                break;
            }
            i++;
        }
        if(flag)printf("GOOD\n");
    }
    return 0;
}

你可能感兴趣的:(POJ 2635 The Embarrassed Cryptographer(素数线性筛法))