题目链接:http://poj.org/problem?id=2635
题目描述:
这题的意思是给定一个很大的整数K,问其所有的质因数中有没有比L小的,有就输出。
主要思想是,先打一张素数表,然后用那些小于L的素数去试,这里注意的是用1000进制对K进行处理,用100进制会超时。
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std ;
const int MAXM = 1000000 ;
bool prim[MAXM+10] ;
int ck[105] ; //记录千分后的数
//选素数
void getPrim(){
for( int i = 0; i < MAXM+5; i++ ){
prim[i] = true ;
}
for( int i = 2; i < MAXM+5; i++ ){
if( prim[i] ){
for( int j = 2*i; j < MAXM+5; j+=i ){
prim[j] = false ;
}
}
}
}
//高精度取余,当存在一个比L小的因数时返回true,否则返回false
bool highMod(int lenth,int t){
int res = 0 ;
for( int i = lenth-1; i >= 0; i-- ){
res = (res*1000+ck[i]) % t ;
}
if( res == 0 ) return true ;
else
return false ;
}
//main
int main(){
int L, num ;
char k[1000] ;
bool judge = false ;
getPrim() ;
while( cin >> k >> L ){
if( L == 0 ) break ;
int len = strlen(k) ;
judge = false ;
for( int i = 0; i < 100; i++ ){
ck[i] = 0 ;
}
//千分化
for( int i = 0; i < len; i++ ){
int ord = (len+2-i) / 3 - 1 ;
ck[ord] = ck[ord]*10 + (k[i]-'0') ;
}
len = (len+2) / 3 ;
for( int i = 2; i < L; i++ ){
if( prim[i] == true ){
if( highMod(len,i) ){
cout << "BAD " << i << endl ;
judge = true ;
break ;
}
}
}
if( !judge ) cout << "GOOD\n" ;
}
return 0 ;
}