【USACO 1.5.2】回文质数

【题目描述】

因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。

写一个程序来找出范围[a,b](5 <= a < b <= 100,000,000)( 一亿)间的所有回文质数;

【格式】

INPUT FORMAT:

(file pprime.in)

第 1 行: 二个整数 a 和 b .

OUTPUT FORMAT:

(file pprime.out)

输出一个回文质数的列表，一行一个。

【分析】

晕，交了3遍才过。（不要鄙视我了==）

先判断回文数后判断质数。

 1 #include <iostream>
 2 #include <cstdio>
 3 #include <cmath>
 4 #include <cstring>
 5 #include <algorithm>
 6 #include <queue>
 7 const int Max=10000000;
 8 using namespace std;
 9 long long a,b,ans[Max],point=0;
10 long long mi[15]={1,10,100,1000,10000,100000,1000000,10000000,100000000,1000000000,10000000000ll};
11 void work(long long num,long long len);//len代表的数字的长度 
12 bool prime(long long num);
13 int main()
14 {
15     //文件操作 
16     freopen("pprime.in","r",stdin);
17     freopen("pprime.out","w",stdout);
18     scanf("%lld%lld",&a,&b); 
19     for (long long i=0;i<=9;i++) 
20     {
21         work(i,1);
22         if (i!=0) work(i*10+i,2);
23     }
24     sort(ans,ans+point);
25     for (long long i=0;i<point;i++) printf("%lld
",ans[i]);
26     return 0;
27 }
28 bool prime(long long num)
29 {
30      for (long long i=2;i<=(long long)sqrt((double)num)+1;i++) if (num%i==0) return 0;
31      return 1;
32 }
33 void work(long long num,long long len)
34 {
35 
36      if (num>b) return;
37      if (num>=a && prime(num)) ans[point++]=num;
38      for (long long i=1;i<=9;i++) 
39      for (long long j=1;j<=5;j++)
40      {
41          if ((len+j*2-1)<=10)
42          work(num*mi[j]+i+i*mi[len+j*2-1],len+j*2);
43      }
44 }