poj-2635（同余模运算）

大致题意：

给定一个大数K，K是两个大素数的乘积的值。

再给定一个int内的数L

问这两个大素数中最小的一个是否小于L，如果小于则输出这个素数。

  

高精度求模+同余模定理

 

1、  Char格式读入K。把K转成千进制Kt，同时变为int型。

把数字往大进制转换能够加快运算效率。若用十进制则耗费很多时间，会TLE。

千进制的性质与十进制相似。

例如，把K=1234567890转成千进制，就变成了：Kt=[  1][234][567][890]。

为了方便处理，我的程序是按“局部有序，全局倒序”模式存放Kt

即Kt=[890][567][234][1  ]  (一个中括号代表一个数组元素)

2、  素数打表，把10^6内的素数全部预打表，在求模时则枚举到小于L为止。

注意打表不能只打到100W，要保证素数表中最大的素数必须大于10^6，否则当L=100W且K为GOOD时，会因为数组越界而RE，这是因为越界后prime都是负无穷的数，枚举的while(prime[pMin]<L)循环会陷入死循环

3、  高精度求模。

主要利用Kt数组和同余模定理。

例如要验证123是否被3整除，只需求模124%3

但当123是一个大数时，就不能直接求，只能通过同余模定理对大数“分块”间接求模

具体做法是：

先求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<iostream>
#include<string.h>
using namespace std;
const int N=1001000;
char c[105];
int f[35];
int prime[1000000];
int l;
int pn;
void setprime()/*素数组打表*/ 
{
    prime[pn++]=2;
    for(int i=3;i<N;i+=2)//奇偶法
    {
        int j;
        for(j=0;prime[j]*prime[j]<=i;j++)
        {
            if(i%prime[j]==0) break;
        }
        if(prime[j]*prime[j]>i){
            prime[pn++]=i;
        }
    }
}
bool mod(int m,int k)//高精度K对p求模，因数检查(整除),又学习了
{
    int ans=0;
    for(int i=k-1;i>=0;i--)
    {
        ans=(ans*1000+f[i])%m;
    }
    if(ans) return false;
    return true;
}
int main()
{
    setprime();
    while(cin>>c>>l)
    {
        if(c[0]=='0'&&l==0) break;
        pn=0;
        memset(f,0,sizeof(f));
        int k=0;
        for(int i=strlen(c)-1;i>=0;i-=3)//把K转换为千进制Kt，其中Kt局部顺序，全局倒序    
        {                                //如K=1234567=[  1][234][567] ，则Kt=[567][234][1  ]
            if(i>=2)
            f[k]+=100*(c[i-2]-48)+10*(c[i-1]-48)+c[i]-48;
            else{
                    int s=1;
                for(int j=i;j>=0;j--)
                {
                    f[k]+=(c[j]-48)*s;
                    s*=10;
                }
            }
                k++;
        }
        int s=0;
        while(prime[s]<l) //枚举prime中比L小的素数
        {
            if(mod(prime[s],k)) {
                cout<<"BAD "<<prime[s]<<endl;
                break;
            }
            s++;
        }
        if(prime[s]>=l) cout<<"GOOD"<<endl;
    }
}