hdoj 素数回文 1431 (模拟)

素数回文

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 16372    Accepted Submission(s): 3621

Problem Description

xiaoou33对既是素数又是回文的数特别感兴趣。比如说151既是素数又是个回文。现在xiaoou333想要你帮助他找出某个范围内的素数回文数，请你写个程序找出 a 跟b 之间满足条件的数。(5 <= a < b <= 100,000,000); 

 

Input

这里有许多组数据，每组包括两组数据a跟b。

 

Output

对每一组数据,按从小到大输出a，b之间所有满足条件的素数回文数（包括a跟b）每组数据之后空一行。

 

Sample Input

  
  

   
   5 500
  
  

 

Sample Output

  
  

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

   
   思路：
  
  
  
  

   
   1、因为是素数，所以，它的最后一位只能是1,3,7,9；
  
  
  
  

   
   2、因为是回文素数，所以前面一半与后面一半数字对称（只用模拟前面一半就行了）；
  
  
  
  

   
   根据这两个规律就可以解决此题：
  
  
  
  

   
   
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<iostream>
#define INF 0x3f3f3f3f
#define ull unsigned long long
#define ll long long
#define IN __int64
#define N 10010
#define M 100000000
using namespace std;
int a[N];
int b[5]={1,3,7,9};
int c[10]={1000000,100000,10000,1000,100,10};
int kk;
void getp()
{
	int i,j,k,l;
	for(i=0;i<4;i++)
	{
		for(j=0;j<10;j++)/7位数的 
		{
			for(k=0;k<10;k++)
			{
				for(l=0;l<10;l++)
				{					
					int flag=0;
					int s=b[i]*c[0]+j*c[1]+k*c[2]+l*c[3]+k*c[4]+j*c[5]+b[i];
					for(int x=2;x<=sqrt(s);x++)
					{
						if(s%x==0)
						{
							flag=1;
							break;
						}
					}
					if(!flag)
						a[kk++]=s;
				}
			}
		}
		for(j=0;j<10;j++)//5、6位数的 
		{
			for(k=0;k<10;k++)
			{
				int flag=0;
				int s=b[i]*c[1]+j*c[2]+k*c[3]+k*c[4]+j*c[5]+b[i];
				for(int x=2;x<=sqrt(s);x++)
				{
					if(s%x==0)
					{
						flag=1;
						break;
					}
				}
				if(!flag)
					a[kk++]=s;
					
				flag=0;
				s=b[i]*c[2]+j*c[3]+k*c[4]+j*c[5]+b[i];
				for(int x=2;x<=sqrt(s);x++)
				{
					if(s%x==0)
					{
						flag=1;
						break;
					}
				}
				if(!flag)
					a[kk++]=s;
			}
		}
		for(j=0;j<10;j++)//3、4位数的 
		{
			int flag=0;
			int s=b[i]*c[3]+j*c[4]+j*c[5]+b[i];
			for(int x=2;x<=sqrt(s);x++)
			{
				if(s%x==0)
				{
					flag=1;
					break;
				}
			}
			if(!flag)
				a[kk++]=s;
				
			flag=0;
			s=b[i]*c[4]+j*c[5]+b[i];
			for(int x=2;x<=sqrt(s);x++)
			{
				if(s%x==0)
				{
					flag=1;
					break;
				}
			}
			if(!flag)
				a[kk++]=s;
		}
	}
	a[kk++]=11;a[kk++]=7;a[kk++]=5;//1、2位数的 
	sort(a,a+kk);
}
int main()
{
	getp();
	int l,r;
	int v=0;
	while(scanf("%d%d",&l,&r)!=EOF)
	{
		int x=lower_bound(a,a+kk,l)-a;
		int y=lower_bound(a,a+kk,r)-a;
		for(int i=x;i<y;i++)
		{
			printf("%d\n",a[i]);
		}
		printf("\n");
	}
	return 0;
}