原题地址:https://www.luogu.org/problem/P1217
题目描述
因为 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
说明
时空限制:1000ms,128M
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;//(处理回文数...)
}
}
}
思路:模拟题,容易TLE,我这里用了快速读入,埃氏筛法,以及把不符合下列条件的数都跳过,所以时间减少了一大半,只需要500多ms。
-
如果一个数是质数,那么除了2,他一定是奇数,所以回文质数一定是奇数。
-
所有偶数位的回文数除了11都不是质数,所以回文质数一定是奇数位。因为所有偶数位的回文数都是11的倍数,所以除了11,他们都不是质数
代码如下:
#include <iostream>
#include <cstdio>
#include <cmath>
using namespace std;
const int N=9999999;
int prime[N];
inline int read(){ //快速读入
char c=getchar();int num=0;
for(;!isdigit(c);c=getchar());
for(;isdigit(c);c=getchar())
num=num*10+c-'0';
return num;
}
bool digit(int x){ //所有偶数位的回文数除了11都不是质数
if((x>=1000&&x<=9999)||(x>=100000&&x<=999999))
return 0;
if(x>=10&&x<=99&&x!=11)
return 0;
return 1;
}
bool huiwen(int x){ //判断是否为回文数
int y=x,num=0;
while(y)
{
num=num*10+y%10;
y/=10;
}
if(num==x)
return 1;
else
return 0;
}
int main(){
for(int i=0;i<=N;i++) //先打表,把所有素数标记为1,合数标记为0
prime[i]=1;
prime[0]=prime[1]=0;
for(int i=2;i<=N;i++){
if(!prime[i])
continue;
for(int j=i*2;j<=N;j+=i)
prime[j]=0;
}
int a,b;
a=read(),b=read();
if(a%2==0) //寻找第一个奇数
a++;
b=min(9999999,b); //题目所给范围最大回文质数是9999999
for(int i=a;i<=b;i+=2){ //遍历所有奇数即可
if(!digit(i)) //跳过位数不符合条件的数
continue;
if(!huiwen(i)) //跳过不是回文数的数
continue;
if(!prime[i]) //跳过不是质数的数
continue;
printf("%d\n",i); //剩下的就都符合条件
}
return 0;
}