poj解题报告——2635

这题特别好理解，就是代码坑啊。题意就是给一个数，这个数是两个大素数的积，再给出一个数，如果最小的素数比给的数大，打印GOOD，否则打印BAD和最小的素数。

这题用的方法是高精度求模+同余模定理。还有素数打表，把10^6内的素数全部预打表，在求模时则枚举到小于L为止。注意打表不能只打到100W，要保证素数表中最大的素数必须大于10^6。千进制的性质与十进制相似。例如，把K=1234567转成千进制，就变成了：Kt=[ 1][234][567]。为了方便处理，按“局部有序，全局倒序”的方式存放Kt中即Kt=[567][234][1 ] 。

代码如下

#include<stdio.h>

#include<string.h>

#include<stdlib.h>

#define Range 1000100

int Kt[10000];

int L;

int prime[Range+1];

void PrimeTable()

{

int pNum=0;

int i,j,flag;

prime[pNum++]=2;

for(i=3;i<=Range;i+=2)

{

flag=1;

for(j=0;prime[j]*prime[j]<=i;j++)

if(!(i%prime[j]))

{

flag=0;

break;

}

if(flag==1)

prime[pNum++]=i;

}

return;

}

int mod(const int* K,const int p,const int len)

{

int sq=0;

int i;

for(i=len-1;i>=0;i--)

sq=(sq*1000+K[i])%p;

if(!sq)

return 0;

return 1;

}

void main()

{

char K[10000];

int i,lenK,pKt,lenKt,pMin,flag;

PrimeTable();

while(scanf("%s %d",K,&L)!=EOF&&L!=0)

{

memset(Kt,0,sizeof(Kt));

lenK=strlen(K);

for(i=0;i<lenK;i++)

{

pKt=(lenK+2-i)/3-1;

Kt[pKt]=Kt[pKt]*10+(K[i]-'0');

}

lenKt=(lenK+2)/3;

flag=1;

pMin=0;

while(prime[pMin]<L)

{

if(mod(Kt,prime[pMin],lenKt)!=1)

{

flag=0;

printf("BAD %d\n",prime[pMin]);

break;

}

pMin++;

}

if(flag==1)

printf("GOOD\n");

}

}