poj2635

最新推荐文章于 2021-02-21 10:14:30 发布

Jiahe Xu

最新推荐文章于 2021-02-21 10:14:30 发布

阅读量486

点赞数

分类专栏： poj

本文链接：https://blog.csdn.net/EOD_realize/article/details/37929189

版权

poj 专栏收录该内容

42 篇文章 0 订阅

订阅专栏

大整数求余，注意是多位的高精度。还有要严格小于给的数。

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
#define maxn 1100000
#define hh printf("orz\n");
#define sys system("pause");
char st[100000];
int prime[maxn],g[10000];
int num,k;
bool vis[maxn];
void init()
{
    for(int i=2;i<=1000100;i++)
    {
        if(!vis[i])
        {
            num++;
            prime[num]=i;
        }
        for(int j=i*2;j<=1000100;j+=i)
        {
            vis[j]=1;
        }
    }
}
bool judge(const int p,const int len)
{
    int sq=0;
    for(int i=len-1;i>=0;i--)  //千进制K是逆序存放
        sq=(sq*1000+g[i])%p;  //同余模定理

    if(!sq)   //K被整除
        return true;
    return false;
}
int main()
{
    init();
    bool flag;
    while(scanf("%s%d",st,&k)&&k)
    {
        memset(g,0,sizeof(g));
        int len=strlen(st);
        for(int i=0;i<len;i++)
        {
            int pKt=(len+2-i)/3-1;
            g[pKt]=g[pKt]*10+(st[i]-'0');
        }
        len=(len+2)/3;
        /*
        for(int i=0;i<len;i++)
        {
            printf("%d ",g[i]);
        }*/
        flag=true;
        for(int i=1;prime[i]<k;i++)
        {
            if(judge(prime[i],len))
            {
                printf("BAD %d\n",prime[i]);
                flag=false;
                break;
            }
        }
        if(flag)
        {
            printf("GOOD\n");
        }
    }
    return 0;
}