[USACO1.5] 回文质数 Prime Palindromes
题目描述
因为 151 151 151 既是一个质数又是一个回文数(从左到右和从右到左是看一样的),所以 151 151 151 是回文质数。
写一个程序来找出范围 [ a , b ] ( 5 ≤ a < b ≤ 100 , 000 , 000 ) [a,b] (5 \le a < b \le 100,000,000) [a,b](5≤a<b≤100,000,000)(一亿)间的所有回文质数。
输入格式
第一行输入两个正整数 a a a 和 b b b。
输出格式
输出一个回文质数的列表,一行一个。
样例 #1
样例输入 #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 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<iostream>
#include<algorithm>
#include<cstring>
using namespace std;
const int N = 1e8+10;
int primes[N],cnt;
bool st[N];
int l,r;
void init(int x){
for(int i=2;i<=x;i++){
if(!st[i])primes[cnt++]=i;
for(int j=0;i*primes[j]<=x;j++){
st[i*primes[j]]=true;
if(i%primes[j]==0)break;
}
}
}
bool check(int x){//判断回文数
int t=x,ans=0;
while(t){
ans=ans*10+t%10;
t/=10;
}
if(ans==x)return true;
return false;
}
int main(){
cin>>l>>r;
if(r<l){
return 0;
}
init(N);
if(l%2==0)l++;
for(int i=l;i<=r;i+=2){
if(!st[i]&&check(i)){
cout<<i<<endl;
}
}
return 0;
}