甲级PAT 1015 Reversible Primes（大量素数经典方法+二分查找）

A reversible prime in any number system is a prime whose "reverse" in that number system is also a prime. For example in the decimal system 73 is a reversible prime because its reverse 37 is also a prime.

Now given any two positive integers N (< 10^5^) and D (1 < D <= 10), you are supposed to tell if N is a reversible prime with radix D.

Input Specification:

The input file consists of several test cases. Each case occupies a line which contains two integers N and D. The input is finished by a negative N.

Output Specification:

For each test case, print in one line "Yes" if N is a reversible prime with radix D, or "No" if not.

Sample Input:

73 10
23 2
23 10
-2

Sample Output:

Yes
Yes
No

题目要求： 

输入多组数据，最后以一个负数结尾。每组数据，第一个整数N（N<10^5》），第二个表示进制，判断N以及N转化进制后的逆序再转化为十进制所得到的数，是否均为素数，若均为素数，则输出Yes,否则，输出No

解题思路：

这里当然先是找到100000以内所有的素数，但是这里并不能用传统的两重循环来找，因为对大量数据循环肯定会超时的。

采用经典的求解素数的方法。用数组isprime[i]作为整数I的标志,若I为素数，isprime = 1，反正为0。先将数组Isprime初始化为1。然后从2开始一直遍历到maxsize，若isprime[i]是素数，则所有在范围以内i的整数倍都不是素数，修改isprime[i的倍数]的值为0。这样用一重循环就找到了所有的素数。然后再循环一遍用prime从下标0开始存储这些素数。

要判断一个数是否为素数，则利用二分查找的方式，查找该数是否在prime数组内

对于一个数N按进制D转化之后逆转再转化为十进制。逆序转换只需要用一个数组trans从下标0开始存储N每次除以D的余数即可，所得到从下标0开始就是逆序后的值。然后将trans数组转换为十进制。

完整代码：

#include<iostream>
#include<string>
using namespace std;
#define maxsize 100001

int isprime[maxsize]; //素数标志位，为1表示是素数，为0表示不是素数 
int prime[maxsize/2]; //存储maxsize以内的所有素数 
int count=0; //用于存储素数的下标，从0开始，一共count+1个素数 

//将整数n转化为d进制，然后逆序转化为十进制，返回结果，这里区分c++原本函数reserve 
int reserve1(int n,int d){
	int trans[14];
	int j,i=0,x=0;
	while(n>0){  
		trans[i++] = n%d; //trans数组中存储的直接就是逆序的值 
		n /= d;
	}
	for(j=0;j<i;j++){
		x = x*d +trans[j]; //转化为十进制 
	}
	return x;
}

//利用二分查找检查x是否为素数，是返回1，不是返回0 
int check(int x){
	int l,r,mid;
	l=0;r=count;
	while(l<=r){
		mid = (l+r)/2;
		if(prime[mid]>x){
			r = mid - 1;
		}else if(prime[mid]<x){
			l = mid +1;
		}else{
			return 1;
		} 
	}
	if(l>r) return 0;
}

int main(){
	int i,j,N,D,temp;
	//传统的方法查找100000以内的素数肯定会超时，这里用经典的求素数方法 
	for(i=0;i<maxsize;i++) isprime[i] = 1;
	for(i=2;i<maxsize/2;i++){
		if(isprime[i]){
			j = 2;
			while(i*j<maxsize){ //素数的整数倍一定不是素数 
				isprime[i*j] = 0;
				j++;
			}
		}
	}
	for(i=2;i<maxsize;i++){
		if(isprime[i]) prime[count++] = i;
	}
	while(scanf("%d",&N)&&N>=0){
		scanf("%d",&D);
		temp = reserve1(N,D);
		if(check(N) && check(temp)){
			cout<<"Yes"<<endl;
		}else{
			cout<<"No"<<endl;
			
		}
	}
	return 0; 
} 

输入的写法也可以是while(scanf("%d %d",&N,&D)>1)，因为scanf返回的是成功输入的数据个数。