The Embarrassed Cryptographer POJ - 2635 同余模+高精度处理 +线性欧拉筛（每n位一起处理）

题意：给出两数乘积K（1e100） 和 一个数L(1e6)  问有没有小于L（不能等于）的素数是K的因数

思路：把数K切割 用1000进制表示   由同余模公式知   k%x=(a*1000%x+b*1000*1000%x+c*1000*1000*1000%x....)

   a b c等为 相应位置的三位数  这样切割可以减少模的次数 防止超时

 

#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
#include<iostream>
using namespace std;
const int maxn=1e6+1000;
int primes[maxn];
int vis[maxn];
int cnt;
void init(){
     cnt=0;
    for(int i=2;i<maxn;i++){
        if(!vis[i])primes[cnt++]=i;
        for(int j=0;j<cnt&&i*primes[j]<maxn;j++){
            vis[i*primes[j]]=1;
            if(i%primes[j]==0)break;
        }
    }
}
char k[1000];
int kt[1000];
int l;
int lenkt;
bool check(int prime){
    int ans=0;
    for(int i=lenkt-1;i>=0;i--){
        ans=(ans*1000+kt[i])%prime;
    }
    if(ans==0)return 1;
    return 0;
}
int main(){
    init();
    while(scanf("%s%d",k,&l)==2&&l&&k[0]!='0'){
        int len=strlen(k);
        memset(kt,0,sizeof(kt));
        for(int i=0;i<len;i++){
            int temp=(len-i+2)/3-1;//从后往前3个3个数 +2向上取整
            kt[temp]=kt[temp]*10+k[i]-'0';
        }
         lenkt=(len+2)/3;//+2是向上取整
         int flag=0;
         for(int i=0;i<cnt&&primes[i]<l;i++){
             if(check(primes[i])){
                 printf("BAD %d\n",primes[i]);
                 flag=1;
                 break;
             }

         }
         if(!flag){
             printf("GOOD\n");
         }
    }
    return 0;
}