hdu 1431 素数回文 （模拟）

链接 ： http://acm.hdu.edu.cn/showproblem.php?pid=1431

解题思路：采用构造的思路来构造回文数字，怎么构造呢？首先将回文数分为两类，奇数长度的回文数字和偶数长度的回文数字。

1）对于偶数长度的回文数字，可以由1-n每一个数字将其对称翻转构造

2）对于奇数长度的回文数字，可以在跳过最后数字一位的情况下翻转构造

还有一个事实就是由同一个数字构造出来的回文数字，偶数长度的要大于所有的奇数长度。

所以要生成第N个回文数，只需要找到他是由谁构造而来即可。

AC代码：

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
//每一个数字都可以构造出两个回文数
ll num[20];//每一数位中所包含的数字个数*2
ll pre[20];//回文数字个数前缀和，按照数位统计
ll ans[20005];
int bit[50];

bool is_prime(ll x)
{
    for(int i=2; i*i<=x; i++)
        if(x%i == 0)
            return false;
    return true;
}

int solve(ll num)//将数字拆分为每一个数位上面的数字
{
    int k = 0,len;
    ll nn = num;
    while(nn)
    {
        nn /= 10;
        ++k;
    }
    len = k;
    while(num)
    {
        bit[k--] = num%10;
        num /= 10;
    }
    return len;
}
int main()
{
    memset(bit,-1,sizeof(bit));
    num[1] = 18;
    pre[1] = 18;
    for(int i = 2; i<=18; i++)
    {
        num[i] = num[i-1]*10;
        pre[i] = pre[i-1] + num[i];
    }
    for(int n = 1; n<= 20005; n++)
    {
        int pos = lower_bound(pre,pre+18,n) - pre;//找到对应第几位可以构造出这么多回文数
        pos--;
        if( (n - pre[pos]) <= num[pos+1]/2)//构造奇数位回文数
        {
            ll number = n - pre[pos]/2;
            int len = solve(number);
            for(int i = 2*len-1; i>len; i--)
                bit[i] = bit[2*len-i];
        }
        else//构造偶数位回文数
        {
            ll number = n - pre[pos]/2  - num[pos+1]/2;
            int len = solve(number);
            for(int i = 2*len; i>len; i--)
                bit[i] = bit[2*len+1-i];
        }
        for(int i = 1; i<=50; i++)//生成回文数
        {
            if(bit[i]!=-1)
                ans[n] = ans[n]*10 + bit[i];
            else
                break;
        }
    }
    int a,b;
    while(cin>>a>>b)
    {
        int l = lower_bound(ans,ans+20000,a) - ans;
        int r = upper_bound(ans,ans+20000,b) - ans;
        for(int i = l; i<=r; i++)
        {
            if(is_prime(ans[i]))
                cout << ans[i] << endl;
        }
        cout << endl;
    }
    return 0;
}