POJ 2635 The Embarrassed Cryptographer(千进制 高精度取模)

POJ 2635

题目大意

给定一个K和L，K是由两个素数组成，如果这两个素数存在比L小(严格小于,因为这个原因WA了两次)的则输出最小的素数,都大于L输出“GOOD”$(4\le K\le 10^{100},2\le L\le 10^{6})​$

分析

L比较小直接枚举L中的素数p，可以打表，然后高精度取模判断K是否是p的整数倍
大数求模可以用千进制+同余来做。
千进制是位了提高效率。比如123456789转换成千进制就是123|456|789
利用同余定理对数进行分块取余。

引用http://blog.csdn.net/lyy289065406/article/details/6648530的例子
以10进制数124%3为例
具体做法是：
先求1%3 = 1
再求（1*10+2）%3 = 0
再求 （0*10+4）% 3 = 1
那么就间接得到124%3=1，这是显然正确的
而且不难发现， （1*10+2）*10+4 = 124
这是在10进制下的做法，千进制也同理，*10改为*1000就可以了

代码

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
using namespace std;
#define uint unsigned int
const int Range=1100000;
int pt[200000];//素数表
int ptcnt=0;
void PrimeTable()
{
      bool bj[Range+1];
      memset(bj,0,sizeof(bj));
      for(int i=2;i<=Range;i++)
      {
            if(bj[i]==0)
            {
                  pt[++ptcnt]=i;
                  for(int j=i*2;j<=Range;j+=i)bj[j]=1;
            }
      }
}
void ToQian(string K,int qian[])//将string K转化成千进制存到数组qian中
{
       int a[105];//先将K存到数组a中
       int lenK=K.size();
       for(int i=lenK-1;i>=0;i--) a[i+1]=K[i]-'0';
       int lenq=ceil((double)lenK/(double)3);
       qian[0]=0;
       int t=0;
       int flag=0;
       for(int i=1;i<=lenK;i++)
       {
             t=t*10+a[i];
             if(flag==0 && i==lenK%3)
             {
                    qian[lenq-(++qian[0])+1]=t;
                    t=0;
                    flag=1;
             }
             else if(i%3==lenK%3)
             {
                   qian[lenq-(++qian[0])+1]=t;
                   t=0;
             }
       }
}
int Mod(int qian[],int p)//千进制取模，返回结果
{
      int t=0;
      for(int i=qian[0];i>=1;i--)
      {
            t=(t*1000+qian[i])%p;
      }
      return t;
}
void Work(int qian[],int L)
{
      for(int i=1;i<=ptcnt;i++)
      {
            if(pt[i]>=L)break;
            if(Mod(qian,pt[i])==0)
            {
                  cout<<"BAD "<<pt[i]<<endl;
                  return ;
            }
      }
      cout<<"GOOD"<<endl;
}
int main()
{
    PrimeTable();
    string  K;
    int L;
    int qian[50];//千进制,qian[0]表示位数

    while(cin>>K>>L )
    {
         if(K.size()==1 && K[0]=='0' && L==0 )break;
         ToQian(K,qian);
         Work(qian,L);
    }
}