poj 2635

http://poj.org/problem?id=2635

好久之前看过的一道题，不过一直没有什么思路， 然后突然想起来， 看了下解题报告， 真的很佩服他们

  10^100太大我们就改成用1000进制存的， 然后用同余模定理

#include <iostream>
#include <stdio.h>
#include <stdlib.h>
#include <cmath>
#include <string.h>


using namespace std;
int knum[1111];
char num[1111];
int prime[1000100];


void PrimeTable( )
{
    int cnt = 0;
    prime[cnt++] = 2;
    for( int i = 3; i < 1000100; i++)
    {
        bool flag = false;
        for( int j = 0; prime[j]*prime[j] <= i; j++)
        {
            if(i%prime[j] == 0)
              {
                  flag = true;
                   break;
              }
        }
        if(!flag)
          prime[cnt++] = i;
    }


}


bool check( int pos, int len)
{
    int mod = 0;
    for( int i = len-1; i >= 0; i--)
      mod = (mod * 1000 + knum[i])%prime[pos];
    if(mod)
      return false;
    return true;
}


void solve()
{
    int k;
    while(scanf("%s%d",&num, &k) != EOF && k)
    {
      int len = strlen(num);
      int bit = 0;
      for( int i = len - 1; i >= 0; i-=3)
      {
          int res = (num[i]-'0') ;
          if(i-1>=0)
            res += (num[i-1]-'0')*10 ;
          if(i-2>=0)
            res += (num[i-2]-'0')*100;
          knum[bit++] = res;
      }
      bool flag = false;
      int pos = 0;
      while(prime[pos] < k)
      {
          if(check(pos,bit))
           {
               printf("BAD %d\n",prime[pos]);
               flag = true;
               break;
           }
           pos++;
      }


     if(!flag)
       printf("GOOD\n");
    }
}
int main()
{
    PrimeTable();
    solve();
    return 0;
}