USACO代码解析PrimePalindromes(pprime).docVIP

描述 因为151既是一个质数又是一个回文数(从左到右和从右到左是看一样的)，所以 151 是回文质数。 写一个程序来找出范围[a,b](5 = a b = 100,000,000)( 一亿)间的所有回文质数; [编辑]格式 PROGRAM NAME: pprime 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 暴力： #includecstdio #includecmath int huiwen(int x) { int n,a=x; int w[50]={0}; n=0; while (a!=0) { w[n]=a%10; a=a/10; n++; } for(int i=0;i=n/2;i++) if(w[i]!=w[n-i-1]) return 0; return 1; } int zhishu(int a) { If( a%2==0) return 0; for(int i=3;i*i=a;i+=2) //不用sqrt！！！！！！！！！ if(a%i==0) return 0; return 1; } int main() { freopen(zhishubiao.txt,w,stdout); int sum=0; for(int i=5;ii++) if(zhishu(i) huiwen(i)) { sum++; printf(%d,,i); } printf(%d,sum); return 0; } Prime Palindromes (pprime) 我的： #include iostream #include fstream using namespace std; ifstream fin(pprime.in); ofstream fout(pprime.out); int a,b; int h[10]; int hui(int n) { int p=0; while(n!=0) { h[p++]=n%10; n/=10; } for(int i=0;i=p/2;i++) { if(h[i]!=h[p-1-i]) return 0; } return 1; } int prime(int n) { if(n==2) return 1; // even only 2 is prime if(n%2==0) return 0; // 2 is the only even divider for(int i=3;i*i=n;i+=2) //from then on its odd { if(n%i==0) return 0; } return 1; } int main() { finab; if(a%2==0) a++; for(int n=a;n=b;n+=2) { if(n%3==0) continue; if(!hui(n)) continue; if(!prime(n)) continue; foutnendl; } return 0; } 就是5-100000000过不了。。。艹