USACO-Section1.5 Prime Palindromes【暴力枚举】

题目描述：

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

INPUT FORMAT:

(file pprime.in)
第 1 行: 二个整数 a 和 b .

OUTPUT FORMAT:

(file pprime.out)
输出一个回文质数的列表，一行一个。

---

SAMPLE INPUT

5 500

---

SAMPLE OUTPUT

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

---

解题思路：

借鉴题目给的思想，枚举回文数再计算（如果用筛法空间占用太大）。偶数回文数能被11整除（注意11符合），下面是代码。

#include<stdio.h>#include<string.h>#include<math.h>#include<stdlib.h>int ans[100010],count=0;int fun(int a){//判断素数    int i;    for(i=2;i<=sqrt(a);i++){        if(a%i==0)return 0;    }    return 1;}int main(){    FILE *fin  = fopen ("pprime.in", "r");    FILE *fout = fopen ("pprime.out", "w");    int a,b;    fscanf(fin,"%d%d",&a,&b);    int i,a1,a2,a3,a4,flag=0,flag1=0;//flag判断是否大于b，flag1判断11是否添加过了    for(a1=0;a1<=9;a1++){//4重循环~        if(flag==1)break;        for(a2=0;a2<=9;a2++){            if(flag==1)break;            for(a3=0;a3<=9;a3++){                if(flag==1)break;                for(a4=0;a4<=9;a4++){                    int temp;                    if(a1){temp=a1*1000000+a2*100000+a3*10000+a4*1000+a3*100+a2*10+a1;}                    else if(a2){temp=a2*10000+a3*1000+a4*100+a3*10+a2;}                    else if(a3){                    if(b>=11&&flag1==0){//将b添加到ans数组，注意限制条件，只添加一次                    if(a<=11)                    ans[count++]=11;                    flag1=1;                    }                    temp=a3*100+a4*10+a3;}                    else temp=a4;                     if(temp>b){//如果大于b就退出循环                    flag=1;                    break;                    }                    else if(fun(temp)&&temp>=a)//如果大于a且是素数就添加                           ans[count++]=temp;                }            }        }    }    for(i=0;i<count;i++){        fprintf(fout,"%d\n",ans[i]);    }    exit(0);}