题目描述
因为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
说明
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的回文数:
#include<iostream>
using namespace std;
bool isprime(int x){ //判断质数的函数,应该不用我多解释
if(x<2)return false;
if(x==2)return true;
for(int i=2;i*i<=x;i++){
if(x%i==0)return false;
}
return true;
}
int l,r,x,y; //l,r是题目范围,x,y是他们的长度
int a[10]; //储存创造出来的数
bool b=true; //用来作为是否跳出的判断
void go(int x,int k){ //k是总位数,x是目前枚举到的位置
if(x==(k+1)/2){ //(k+1)/2是k位的中间数,可以自己算一下
for(int i=k;i>x;i--){ //把剩下的位数完善
a[i]=a[k-i+1];
}
int shu=0; //将数组转化成数
for(int i=1;i<=k;i++){
shu=shu*10+a[i];
}
if(shu<l)return; //小于l就跳过
if(shu>r){ //大于r就跳出
b=false;
return;
}
if(isprime(shu))cout<<shu<<endl; //是质数就输出
return; //这个一定不能漏,要不然会一直搜下去(本人是亲身体验过的受害者)
}
int i;
if(x)i=0;
else i=1; //最高位不能是0
for(i=i;i<=9;i++){ //dfs(深度优先搜索)下一位,没学过可以去看下P1706
if(b==false)return; //如果已经出现大于r的数,就跳出函数
a[x+1]=i;
go(x+1,k);
}
return;
}
int weishu(int j){ //计算一个数的长度
int b=j,count=0;
while(b>0){
b/=10;
count++;
}
return count;
}
int main(){
cin>>l>>r;
x=weishu(l);
y=weishu(r);
for(int i=x;i<=y;i++){
if(i==1){ //一位数特判
if(l<=5&&r>=5)cout<<5<<endl;
if(l<=7&&r>=7)cout<<7<<endl;
continue;
}
if(i==2){ //两位数特判
if(l<=11&&r>=11)cout<<11<<endl;
continue;
}
if(i%2==0)continue; //偶数位就跳过
if(i==9)break; //九位直接跳出
b=true; //每次搜索要先重置这个变量
go(0,i); //从第0位开始搜
}
return 0;
}
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;//(处理回文数...)
}
}
}