USACO-Section1.5 Prime Palindromes

2017-07-19

题目大意：

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

样例输入：

5 500

样例输出：

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

题解：

这道题我一开始使用找回文数，再判断是不是质数，然后提交的时候，发现最后一组测试数据死也过不了。。
最后一组测试数据是5 100000000。好像过不了是很正常的，毕竟算法的时间复杂度过高。后来看了大佬的博客，发现不是一定要找出回文数的，可以生成回文数，这样10000000（一亿）的数据只需要四位数就能生成八位的回文数，计算量大大减小。

代码：

#include<fstream>
#include<cmath>
#include<cstring>
#include<vector>
#include<algorithm>
using namespace std;
vector<int> v;
int flag = false;
ofstream cout("pprime.out");
ifstream cin("pprime.in");

bool is_Prime(int _num){
    if(_num == 2 || _num == 3){
        return true;
    }
    if(_num % 6 != 1 && _num % 6 != 5){
        return false;
    }
    int temp = sqrt(_num);
    for(int i = 5;i <= temp;i += 6){
        if(_num % i == 0 || _num % (i + 2) == 0){
            return false;
        }
    }
    return true;
}

void create_Palindromes(int m ,int n){
    int tmp1 , tmp2;
    for(int i = 1;i <= n;i++){
         tmp1 = i , tmp2 = i;
         //奇数位数
         tmp2 /= 10;
         while(tmp2){
             tmp1 *= 10;
             tmp1 += tmp2 % 10;
             tmp2 /= 10;
         }
         if(tmp1 > n) break;
         if(tmp1 >= m && is_Prime(tmp1)){
             v.push_back(tmp1);
         }
         //偶数位数
         if(!flag){
            tmp1 = i , tmp2 = i;
            while(tmp2){
                tmp1 *= 10;
                tmp1 += tmp2 % 10;
                tmp2 /= 10;
            }
            if(tmp1 > n) flag = true;
            if(tmp1 >= m && is_Prime(tmp1)){
                v.push_back(tmp1);
            }
         }

    }
}

int main(){
    int m , n;
    cin >> m >> n;
    create_Palindromes(m,n);
    sort(v.begin(),v.end());
    for(vector<int>::iterator it = v.begin();it != v.end();it++){
        cout << *it << endl;
    }
    return 0;
}