POJ2635-The Embarrassed Cryptographer-大整数素因子

计算一个大整数（10^100）中有没有一个小于L的素因子。这个大整数是两个素数的乘积。其实就是RSA加密。

做法是把大整数表示成千进位，用数组存储，然后一位一位地取模。

/*--------------------------------------------------------------------------------------*/
//        Helica's header
//        Second Editions
//        2015.11.7
//
#include <algorithm>
#include <iostream>
#include <cstring>
#include <ctype.h>
#include <cstdlib>
#include <cstdio>
#include <vector>
#include <string>
#include <queue>
#include <stack>
#include <cmath>
#include <set>
#include <map>

//debug function for a N*M array
#define debug_map(N,M,G) printf("\n");for(int i=0;i<(N);i++)\
{for(int j=;j<(M);j++){\
printf("%d",G[i][j]);}printf("\n");}
//debug function for int,float,double,etc.
#define debug_var(X) cout<<#X"="<<X<<endl;
/*--------------------------------------------------------------------------------------*/
using namespace std;

int N,M,T,L;
char save[];
int num[];

const int MAXN = ;
int prime[MAXN];

void getPrime()
{
    memset(prime,,sizeof prime);
    for(int i=;i<=MAXN;i++)
    {
        if(!prime[i]) prime[++prime[]] = i;
        for(int j=;j<=prime[]&&prime[j]<=MAXN/i;j++)
        {
            prime[prime[j]*i] = ;
            if(i%prime[j] == ) break;
        }
    }
}

bool mod(int p,int len)
{
    int r = ;
    for(int i=len-;i>=;i--)
        r = (r* + num[i])%p;

    if(!r) return false;
    else return true;
}

int main()
{
    getPrime();

    while(scanf("%s%d",save,&L) && L)
    {
        int flag = ;
        int len = strlen(save);

        memset(num,,sizeof num);

        for(int i=;i<len;i++)
        {
            int p = (len+-i)/-;
            num[p] = num[p]*+(save[i]-'');
        }
        len = (len+)/;

        for(int i=;prime[i]<L;i++)
        {
            int p = prime[i];
            //printf("p=%d\n",p);
            if(!mod(p,len))
            {
                printf("BAD %d\n",p);
                flag = false;
                break;
            }
        }
        if(flag) printf("GOOD\n");
    }
}