section 1.5 pprime

-先生成回文数
-生成回文数的迭代函数比较复杂，因为是新手，这里处理应该没那么好
-质数判定，先列出sqrt(b)以内的质数表（100,000,000最多只需要列10,000以内），用待测数N依次除以这些质数，求余
-生成质数表的函数很通用
-没有必要用Miller-Rabin算法

/*
ID: penglin3
PROG: pprime
LANG: C++11
*/
#include <iostream>
#include <fstream>
#include <vector>
#include <math.h>
using namespace std;//fstream 需要std合法范围


vector<int> viPrime{};//质数表
vector<long> answer{};//答案表

//建质数表的函数,数n:3,5,7...maxnum数组最大下标,然后写入vbPrime
void listPrime(vector<int> &viPrime, const int& maxnum);

//判断是否为质数的函数
bool judgePrime(const vector<int>&viPrime, const long& num);

//求num的位数
int numbit(long num);//这里必须用拷贝参数

//生成回文数pal的第n位并填入，如果生成完成，将pal带入到f中去。N为总位数,b为生成的最大数
void Palindromes(const int& n, long pal, const int& N, const long& b, void f(long pal));

//生成回文数之后的处理函数
void select(long pal);

int main() {
    ifstream fin("pprime.in");
    ofstream fout("pprime.out");

    //数据输入
    long a, b;
    fin >> a >> b;//5 <= a, b <= 100,000,000

    //质数表,大于2的偶数值虽然是false，但是之后操作不会遍历到这些数，
    //可以当做无关项

    listPrime(viPrime, sqrt(b));

    int bita=numbit(a),bitb = numbit(b);
    long pal = 0;

    for(int i=bita; i<=bitb; ++i)
        Palindromes(1, pal, i, b, select);

    //输出
    for (auto &itan : answer)
        if(itan >= a) fout << itan << endl;

    fin.close();
    fout.close();
    return 0;
}

//建质数表的函数,数n:3,5,7...maxnum数组最大下标,然后写入vbPrime
void listPrime(vector<int> &viPrime, const int& maxnum) {
    bool bPrime[10001] = {};
    bPrime[0] = true;
    bPrime[1] = true;

    int n = 3;
    while (n * n < maxnum) {
        for (int i = n * n; i <= maxnum; i += n) {
                bPrime[i] = true;
            }
        do//这里必须用do-while
            n += 2;
        while (n * n < maxnum && bPrime[n]);
    }

    for (int i = 3; i <= maxnum; ++i)
        if (i % 2 && !bPrime[i])viPrime.push_back(i);

    return;
}

/*---------------------------------------------------------*/

//判断是否为素数,直接从3开始，因为本题传入的是大于2的奇数
bool judgePrime(const vector<int>&viPrime, const long& num) {
    for (int i = 0; i != viPrime.size(); ++i) {
        if (viPrime[i] > sqrt(num)) break;
        else if (!(num % viPrime[i]))return false;
    }
    return true;
}

/*---------------------------------------------------------*/

//求num的位数
int numbit(long num) {
    int i = 0;
    while (num) {
        i++;
        num /= 10;
    }
    return i;
}

/*---------------------------------------------------------*/

//生成回文数pal的第n位并填入，如果生成完成，将pal带入到f中去。N为总位数,b为生成的最大数
void Palindromes(const int& n,long pal, const int& N, const long& b, void f(long pal)) {
    if (pal > b) return;
    //回文数第一位
    if (n == 1) {
        long pal0 =(N == 1)? 1: pow(10, N - n) + 1;
        for (int i = 1; i <= 9; i += 2) {
            Palindromes(n + 1, pal + i * pal0, N, b, f);
        }
    }
    //奇数回文数中间项
    else if (N % 2 && n == N / 2 + 1) {
        long pal0 = pow(10, n - 1);
        for (int i = 0; i <= 9; ++i) {
            Palindromes(n + 1, pal + i * pal0, N, b, f);
        }
    }
    //回文数补全
    else if (n <= N / 2) {
        long pal0 = pow(10, n - 1) + pow(10, N - n);
        for (int i = 0; i <= 9; ++i) {
            Palindromes(n + 1, pal + i * pal0, N, b, f);
        }
    }
    else {
        f(pal);
    }
    return; 
}

/*--------------------------------------------------------*/
//生成回文数之后的处理函数
void select(long pal) {
    if (judgePrime(viPrime, pal)) answer.push_back(pal);
    return;
}