SDIBT-OJ 2318 回文质数

Description

因为151即是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 号是回文质数。 写一个程序来找出范围[a,b](5 <= a < b <= 100,000,000)间的所有回文质数;

Input

第 1 行: 二个整数 a 和 b

Output

输出一个回文质数的列表，一行一个。

Sample Input

5 500

Sample Output

5
7
11
101
131
151
181
191
313
353
373
383

思路：

模拟题（自己在模拟题方面还是比较弱的啊）

求回文素数，如果枚举所有的数，先判断是不是素数，在判断是不是回文数，会很浪费时间

巧妙地方法是：先构造出回文数，然后判断是不是素数即可（回文的数量比素数的数量要少）

回文数是呈对称的状态，所以可以先构造出左半边，对称的构造出右半边，再判断是不是素数即可

 

回文质数满足的条件：

1、偶数位除了 11 是回文素数之外，其余的偶数位不存在回文数是素数，因为他们都能被 11 整除，所以不是素数

      换句话说：回文素数必须满足奇数位，以中间一位为对称轴

2、在奇数位回文数中，首位为2、4、6、8的数均不是素数。

     首位是它们，根据回文数的性质，末尾也是他们，但不是素数

3、分析奇数位的回文素数：

     （1）、个位数的奇数只有  2、3、5、7

                  偶数位只有 11,可以特判 输出这几个数

     （2）、3 位的回文素数   aba

                  直接模拟即可，  a 不能取   0 （首位），b 可以取 0～9

                  类似的模拟 5位 、7 位 的即可

                 9 位已经超出了范围

4、因 5 的任何末尾数是 5, 所以在奇数位回文数中，首位为 5  的数均不是素数               

CODE：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <cstdlib>
#include <algorithm>

using namespace std;
#define memset(a,v) memset(a,v,sizeof(v))

typedef long long LL;
int cnt=0;

LL ans[1000000+10];
LL flag=0;
LL a,b;

int prime(LL x)
{
    if(x==1)
        return 0;
    for(LL i=2; i<=sqrt(x); i++)
        if(x%i==0)
            return 0;
    return 1;
}

void Creat()
{
    for(int i=1; i<=9; i++)
        for(int j=0; j<=9; j++)
            ans[flag++] = i*100+j*10+i;

    for(int i=1; i<=9; i+=2)
        for(int j=0; j<=9; j++)
            for(int k=0; k<=9; k++)
                ans[flag++]=i*10000+j*1000+k*100+j*10+i;

    for(int i=1; i<=9; i++)
        for(int j=0; j<=9 ; j++)
            for(int k=0; k<=9; k++)
                for(int l=0; l<=9; l++)
                    ans[flag++]=i*1000000+j*100000+k*10000+l*1000+k*100+j*10+i;

    for(int i=0; i<flag&&ans[i]<=b; i++)
    {
        if(prime(ans[i])&&ans[i]>=a) // 这里一定要判断是不是大于 a，注意边界范围
        {
            printf("%lld\n",ans[i]);
        }
    }

}


int main()
{
    memset(ans,0);

    scanf("%lld %lld",&a,&b);
    if(a<=11)
    {
        for(int i=5; i<=11; i++)
        {
            if(i>=a&&i<=b)
            {
                if(prime(i))
                {
                    printf("%d\n",i);
                }
            }
        }
    }
    if(b>11)
    {
        Creat();
    }
}