题目描述
因为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的回文数:
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<bits/stdc++.h>
using namespace std;
bool prime[100000005];
void init(int b)
{
// 用埃氏筛法生成质数表
//初始化,默认全部都是质数
memset(prime, true, sizeof(prime));
prime[1]=false;//1不是质数
int n=sqrt(b);
for(int i=2;i<=n;i++)
if(prime[i])
{
//质数的倍数绝对不是质数,把所有质数的倍数全部设为false
for (int j=2;j<=b/i;j++)//i*j<=b
prime[i*j]=false;
}
}
bool Palin(int x)//判断是不是回文数
{
int t=0,y=x;
while(y)
{
t*=10;
t+=y%10;
y/=10;
}
if(t==x) return true;
else return false;
}
int main()
{
int a,b;
cin>>a>>b;
if (b>=10000000) b=9999999;//一定要加这个,不然最后一个测试点过不了
init(b);
for(int i=a;i<=b;i++)
if(prime[i]&&Palin(i))
cout<<i<<endl;
return 0;
}