回文质数---（质数的产生）

题目描述

因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。

写一个程序来找出范围[a,b](5 <= a < b <= 100,000,000)( 一亿)间的所有回文质数;

输入输出格式

输入格式：

 

第 1 行: 二个整数 a 和 b .

 

输出格式：

 

输出一个回文质数的列表，一行一个。

 

输入输出样例

输入样例#1： 复制

5 500

输出样例#1： 复制

5
7
11
101
131
151
181
191
313
353
373
383

说明

Hint 1: Generate the palindromes and see if they are prime.

提示 1: 找出所有的回文数再判断它们是不是质数（素数）.

Hint 2: Generate palindromes by combining digits properly. You might need more than one of the loops like below.

提示 2: 要产生正确的回文数，你可能需要几个像下面这样的循环。

题目翻译来自NOCOW。

USACO Training Section 1.5

产生长度为5的回文数:

for (d1 = 1; d1 <= 9; d1+=2) {    // 只有奇数才会是素数
     for (d2 = 0; d2 <= 9; d2++) {
         for (d3 = 0; d3 <= 9; d3++) {
           palindrome = 10000*d1 + 1000*d2 +100*d3 + 10*d2 + d1;//(处理回文数...)
         }
     }
 }

#include<stdio.h>
#define N_7 10000000
#define N_6 1000000
#define N_5 100000
#define N_4 10000
int isprime(int n){
    int i;
	for(i=2;i*i<=n;i++)
	if(n%i==0)
	return 0;
	return 1;
}

int main(){
    long long int d1,d2,d3,d4,a,b,num,t,n=0;
	scanf("%lld %lld",&a,&b);
	t=b;
	while(b>0){
	    b/=10;
		n++;
	}
	b=t;
	if(n>=1)
	  for(d1=5;d1<=9;d1+=2){
		if(d1<=b&&d1>=a&&isprime(d1))
		   printf("%lld\n",d1);
		else if(d1>b)
		         goto L;
	}
	if(n>=2)
	 for(d1=1;d1<=9;d1+=2){
	        if(d1*10+d1<=b&&d1*10+d1>=a&&isprime(d1*10+d1))
			     printf("%lld\n",d1*10+d1);
		     else if(d1*10+d1>b)
			      goto L;
	 }
     if(n>=3){
	    for(d1=1;d1<=9;d1+=2)
		for(d2=0;d2<=9;d2++){
		num=d1*100+d2*10+d1;
		if(num<=b&&num>=a&&isprime(num))
		    printf("%lld\n",num);
			else if(num>b)
			goto L;
		}
		for(d1=1;d1<=9;d1+=2)
		for(d2=0;d2<=9;d2++){
		num=d1*1000+d2*100+d2*10+d1;
		if(num<=b&&num>=a&&isprime(num))
		    printf("%lld\n",num);
			else if(num>b)
			goto L;
		}
		}
		if(n>=5){
		    for(d1=1;d1<=9;d1+=2)
			for(d2=0;d2<=9;d2++)
			  for(d3=0;d3<=9;d3++){
			     num=d1*N_4+d2*1000+d3*100+d2*10+d1;
				 if(num<=b&&num>=a&&isprime(num))
		    printf("%lld\n",num);
			else if(num>b)
			goto L;
			  }
			  for(d1=1;d1<=9;d1+=2)
			for(d2=0;d2<=9;d2++)
			  for(d3=0;d3<=9;d3++){
			     num=d1*N_5+d2*10000+d3*1000+d3*100+d2*10+d1;
				 if(num<=b&&num>=a&&isprime(num))
		    printf("%lld\n",num);
			else if(num>b)
			goto L;
			  }
			}
			if(n>=7){
			for(d1=1;d1<=9;d1+=2)
	         for(d2=0;d2<=9;d2++)
		       for(d3=0;d3<=9;d3++)
                for(d4=0;d4<=9;d4++){
				num=d1*N_6+d2*N_5+d3*N_4+d4*1000+d3*100+d2*10+d1;
				  if(num<=b&&num>=a&&isprime(num))
				  printf("%lld\n",num);
				  else if(num>b)
				    goto L;
				}
				for(d1=1;d1<=9;d1+=2)
	            for(d2=0;d2<=9;d2++)
		          for(d3=0;d3<=9;d3++)
                  for(d4=0;d4<=9;d4++){
				  num=d1*N_7+d2*N_6+d3*N_5+d4*N_4+d4*1000+d3*100+d2*10+d1;
				  if(num<=b&&num>=a&&isprime(num))
				  printf("%lld\n",num);
				  else if(num>b)
				    goto L;
				}
			}
	L:;
    return 0;	
}