P1217 [USACO1.5]回文质数 Prime Palindromes(洛谷题解)

---

题目描述：因为 151 既是一个质数又是一个回文数（从左到右和从右到左是看一样的），所以 151 是回文质数。

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

题目提示：
提示 1: 找出所有的回文数再判断它们是不是质数（素数）.
提示 2: 要产生正确的回文数，你可能需要几个像下面这样的循环.

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;//(处理回文数...)
         }
     }
 }

---

一、问题分析

分析题目，题目中要求我们输出一个又是质数又是回文数的一个数字，大部分同志可能像我一样，首先想到的就是使用一个for循环，再嵌套一个while循环将for中的每一个数字倒置，并判断是否相等来计算回文数

代码如下

for (int i = a; i <= b; i++)
{
	int sum = 0, c = i;
	while (c > 0)	//数字倒置
	{
		sum = sum * 10 + c % 10;
		c /= 10;
	}
	if (sum == i)huiwen[k++] = i;
}

但这样做的问题是，在 $(5\le a<b\le 100,000,000)$的取值范围中，当数字变得特别大的时候，每一次计算回文数就要花费大量的时间，而且在后面判断数组中存储的数字是否是质数的时候也需要用到循环遍历，很容易就会被oj判定为超时。

文章作者尝试过使用提前将回文数存储在一个给定数组里面，再进行质数的计算，或者引入各种博客上dalao们写的快速判断质数的方法之后，还是发现会超时。

所以解决超时的问题关键在于回文数的判断上，而不在于质数。

二、解决方式

1.怎样解决回文数问题？

与前面的代码思路相反的是，我们可以利用回文数的特点“主动出击”，根据题目中的提示，在代码中主动生成回文数，并将其存储在回文数组中

值得注意的是，根据回文数的定义，在题目中所给的数值判断范围中，回文数只有1，3，5，7位，而不存在4，6，8位的回文数（因为回文数从左到右读和从右到左读是一样的，具有对称性）

2.考虑质数问题

考虑以下几点

• 偶数不是质数，在生成回文数的第一位和最后一位d1时，可以直接从1，3，5，7，9算起
• 质数n的判断范围从2到$\sqrt{n}$（具体可参考dalao写的文章）
  C++判断是否是素数（质数）

三、代码

根据上述分析，可得代码如下：

#include <iostream>
#include <math.h>
using namespace std;

int isPrime(int s){//判断质数 
	for(int i = 3;i<=sqrt(s);i+=2)
	//由于生成回文数时已经不考虑偶数，可直接从奇数作因子判断 
		if(s % i == 0) return 0;
	return 1;
}

int main(int argc, char** argv) {
	int a,b;
	cin>>a>>b;
	
	int h[100005] = {5,7,11};
	//一位的回文质数易知，可直接初始化 
	int k = 3;
	int d1,d2,d3,d4;
	//初始化四个整数用于遍历,只存在下面几种组合 
	//百万 十万 万 千 百  十  个 
	//              d1  d2 d1
	//        d1 d2 d3  d2 d1 
	//d1   d2 d3 d4 d3  d2 d1
	
	//看看输入b是否要求咱们要计算这么大的数
	//如果没到要求，咱就不算了 
	if(b>=100)	//d1 d2 d1
		for(d1 = 1; d1 <= 9; d1 += 2)	//偶数非质数，直接从奇数开始生成 
			for(d2 = 0; d2<= 9; d2++){
				int res = d1*100+d2*10+d1;
				if(isPrime(res))
					h[k++] = res;	//回文质数预存数组	
			}
	
	if(b>=10000)	//d1 d2 d3 d2 d1
		for(d1 = 1;d1<=9;d1+=2)
			for(d2=0;d2<=9;d2++)
				for(d3 = 0;d3<=9;d3++){
					int res = d1*10000+d2*1000+d3*100+d2*10+d1;
					if(isPrime(res))
						h[k++] = res;	//回文质数预存数组	
				}
	
	if(b>=1000000)	//d1 d2 d3 d4 d3 d2 d1
		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++){
						int res = d1*1000000+d2*100000+d3*10000+d4*1000+d3*100+d2*10+d1;
						if(isPrime(res))
							h[k++] = res;	//回文质数预存数组	
					}
					
	for(int j =0;j<k;j++){//判断数组中所存数字是不是在给定ab的范围内
	//k是数组内存储的回文质数的个数
		if(h[j]>=a&&h[j]<=b){
			cout<<h[j]<<endl;
		}
	}
			
	return 0;
}

四、调试与运行结果

input:

5 500

output:

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

洛谷oj：

---

写在最后

本文作者才疏学浅，若有纰漏，恳请斧正

本文章仅用于各位同志作为学习交流之用，不作任何商业用途，望悉知