洛谷P1217

题目

原题地址
题目描述
因为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;//(处理回文数...)
         }
     }
 }

题解

这题咋被分在递归里了。。。
我直接生成了5~100,000,000以内所有的回文数（我觉得可能会比生成素数快一点），然后再判断是不是在ab之间和是不是素数，输出。其实就是个自动打表的方法，本来用for(int count = 2; count < sqrt(num); count++) if(num%count == 0) return false;会TLE，然后我去网上找了个求素数的快捷办法然后就被我水过去了2333333333
话说我看到洛谷题解里好多直接在程序里打了个静态表的，我还能说什么呢。。。真是猛士。。。。。果然打表大发好233333333

代码

#include 
#include 
using namespace std;
int palindromes[100000];

bool isPrime2(int n){
    if((n == 2) || (n == 3))
        return true;
    if((n == 1) || ((n & 1) == 0) || ((n % 3) ==0))
        return false;
    int k = (int)sqrt(n);
    int step = 4;
    for(int i = 5; i <= k; i += step)
    {
        if(n % i == 0)
            return false;
        step = 6 - step;
    }
    return true;
}
void generate(){
    int d1 = 0, d2 = 0, d3 = 0, d4 = 0, count = 0;
    for (d1 = 5; d1 <= 9; d1+=2){
        palindromes[count] = d1;
        count++;
    }
    for (d1 = 1; d1 <= 9; d1+=2){
        palindromes[count] = d1*10 + d1;
        count++;
    }
    for (d1 = 1; d1 <= 9; d1+=2) {
        for (d2 = 0; d2 <= 9; d2++) {
            palindromes[count] = 100*d1 + 10*d2 + d1;
            count++;
        }
    }
    for (d1 = 1; d1 <= 9; d1+=2) {
        for (d2 = 0; d2 <= 9; d2++) {
            palindromes[count] = 1000*d1 + 100*d2 +10*d2 + d1;
            count++;
        }
    }
    for (d1 = 1; d1 <= 9; d1+=2) {
        for (d2 = 0; d2 <= 9; d2++) {
            for (d3 = 0; d3 <= 9; d3++) {
                palindromes[count] = 10000*d1 + 1000*d2 +100*d3 + 10*d2 + d1;
                count++;
            }
        }
    }
    for (d1 = 1; d1 <= 9; d1+=2) {
        for (d2 = 0; d2 <= 9; d2++) {
            for (d3 = 0; d3 <= 9; d3++) {
                palindromes[count] = 100000*d1 + 10000*d2 +1000*d3 + 100*d3 +10*d2 + d1;
                count++;
            }
        }
    }
    for (d1 = 1; d1 <= 9; d1+=2) {
        for (d2 = 0; d2 <= 9; d2++) {
            for (d3 = 0; d3 <= 9; d3++) {
                for (d4 = 0; d4 <= 9; d4++) {
                    palindromes[count] = 1000000*d1 + 100000*d2 + 10000*d3 + 1000*d4 + 100*d3 + 10*d2 + d1;
                    count++;
                }
            }
        }
    }
    for (d1 = 1; d1 <= 9; d1+=2) {
        for (d2 = 0; d2 <= 9; d2++) {
            for (d3 = 0; d3 <= 9; d3++) {
                for (d4 = 0; d4 <= 9; d4++) {
                    palindromes[count] = 10000000*d1 + 1000000*d2 + 100000*d3 + 10000*d4 + 1000*d4+ 100*d3 + 10*d2 + d1;
                    count++;
                }
            }
        }
    }

}
int main(int argc, const char * argv[]) {
    int a = -1, b = -1;
    cin >> a >> b;
    generate();
    for(int num: palindromes)
        if((num != 0) && (num>=a) && (num<= b) && isPrime2(num))
            cout << num << endl;
    return 0;
}