P1217 [USACO1.5]回文质数 Prime Palindromes【暴力枚举、欧拉筛】

题目描述

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

写一个程序来找出范围 [a,b] (5 \le a < b \le 100,000,000)[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;//(处理回文数...)
         }
     }
 }

解题思路

这里，我们可以先打表保存1~b范围内的所有素数，然后从找出的素数中再判断哪些是回文。

这里可以先缩小规模，对于偶数位的回文串一定可以被11整除，所以偶数位的回文串一定不是素数。

证明：

观察偶数位的回文数，提取所有奇数位的数字，与提取所有偶数位的数字，正好是相反的顺序。
因此，偶数位数和等于奇数位数和，从而差等于0，而0能被11整除，因此这个回文数，可以被11整除

这样，对于b>=10000000的情况下我们只需要判断b==9999999即可。此外，这里我采用欧拉筛来打素数表（用埃氏筛也可以），其具体算法见：埃氏筛、欧拉筛

代码

#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;
#define MAXN 10000000

int prime[MAXN], p[MAXN] = {0};
int pNum = 0;
int a, b;

void Find_Prime()   //欧拉筛
{
    for(int i = 2;i<=b;i++)
    {
        if(p[i] == 0)
            prime[++pNum] = i;
        for(int j=1;j<=pNum&&i*prime[j]<=b;j++)  //筛去所有i的倍数
        {
			p[i*prime[j]] = 1;
            if(i%prime[j]==0)
				break;
		}
    }
}

int isPalindromes(int num)   //判断回文数，将要判断的数倒过来比较
{
	int tmp = num, num2 = 0;
	while (tmp)
	{
		num2 = num2 * 10 + tmp % 10;
		tmp /= 10;
	}
	if (num2 == num)
		return 1;
	else return 0;
}

int main()
{
	scanf("%d%d", &a, &b);
	if (b >= 10000000) b = 9999999;
	Find_Prime();    //先查找1~b以内的素数
	for (int i = a; i <= b; i++)
	{
		if (p[i] == 0 && isPalindromes(i))
			printf("%d\n", i);
	}
	return 0;
}