USACO Prime Palindromes

两个double之间取余，要用到math.h里的函数

取整与取余
double modf (double, double*); 将参数的整数部分通过指针回传, 返回小数部分
double fmod (double, double); 返回两参数相除的余数

两个double之间判断是否相等，不能用==，目前最优的方案是做差。 
if   (fabs(result1   -   result2)   <   FLT_EPSILON) 
        printf( "equals\n "); 
详见 45楼 http://topic.csdn.net/t/20061104/21/5133629.html#

判断一个int数是否为回文

int main()
{
	int m,n,i,j,k,c=0;

	cin>>m>>n;
	for(i=m;i<=n;i++)
	{
		for(j=i,k=0;j;j/=10) //这个循环就是求反序数的部分
			k=k*10+j%10;
		if(k==i)
			c++;
	}
	cout<<c<<endl;
	return 0;
}

取一个double的整数部分——
1、直接赋值给整数变量。如：  
 inti=2.5;或i=(int)2.5;  
 这种方法采用的是四舍五入，不适用于你的问题。  
 2、使用floor函数。floor(x)返回的是x的整数部分。如：  
 floor(2.5)=2  
 floor(-2.5)=-2  
 3、使用ceil函数。ceil(x)返回的是不大于x的最小整数。如：  
 ceil(2.5)=2  
 ceil(-2.5)=-3  
 floor和ceil对于正数没有区别，但是对于负数结果不同。

如果你只是想简单的不用科学计算法来表示一个数，那么在输出的时候用 

cout<<std::fixed<<yournumber; 

这样就可以了，注意要包含<iomanip>头文件

我的第一版本代码（cao，今天才知道c++里int类型可以表达到2^32,原来只能表示到3万多的那个是short int），从a到b，一个个先判断是否回文数，然后是否质数。test10超时……被无情的鄙视了

/*
ID: wangxin12
PROG: pprime
LANG: C++
*/

#include <iostream>
#include <fstream>
#include <vector>
#include <math.h>
#include <string>

using namespace std;

int A, B;
vector<int> ans;

bool isPalindrome(int n);
bool isPrime(int n);

int main() {
	ifstream fin("pprime.in");
	fin>>A>>B;
	cout<<A<<" "<<B<<endl;
	fin.close();

	for(int i = A; i <= B; i += 1) {
		if(isPalindrome(i)) {
			if(isPrime(i)) ans.push_back(i);
		}
	}


	ofstream fout("pprime.out");
	int j = 0;
	while(j < ans.size()) {
		fout<<std::fixed<<ans[j]<<endl;
		j += 1;
	}

	fout.close();

    return 0;
}

bool isPalindrome(int index) {
	bool flag = false;
	
	int m = 0, n = index;
	for(m = 0; n > 0; n = n / 10) {
		m = m * 10 + n % 10;
	}
	if( m == index ) flag = true;

	return flag;
}

//bool isPalindrome(double index) {
//	bool flag = false;
//	
//	double m = 0, n = index;
//	for(m = 0; fabs(n - 0) > 0.0000001; n = floor(n / 10)) {
//		m = m * 10 + fmod(n, 10);
//	}
//	if( fabs(m - index) < 0.000001 ) flag = true;
//
//	return flag;
//}

bool isPrime(int n) {
	bool flag = true;

	int start = 2;
	int end = sqrt((double)n);
	for(int i = start; i <= end; i++) {
		int mod = n % i;
		if( mod == 0) {
			flag = false;
			break;
		}
	}

	return flag;
}

//bool isPrime(double n) {
//	bool flag = true;
//
//	double start = 2;
//	double end = floor(sqrt(n));
//	for(double i = start; i <= end; i++) {
//		double mod = fmod(n ,i);
//		if( fabs(mod - 0) < 0.0000001) {
//			flag = false;
//			break;
//		}
//	}
//
//	return flag;
//}

优化方向
http://blog.csdn.net/sjtuyunlei/article/details/6629357
http://acm.swust.edu.cn/algwiki/doku.php?id=usaco:1:5:prime_palindromes

/*
ID: wangxin12
PROG: pprime
LANG: C++
*/

#include <iostream>
#include <fstream>
#include <vector>
#include <math.h>
#include <string>

using namespace std;

vector<int> palindromes;
int A, B;

void genOdd();
void genEven();
bool isPrime(int n);
void Sort(vector<int> &vec);
void Merge(vector<int> &vec, vector<int> &v1, vector<int> &v2);

int main() {
	ifstream fin("pprime.in");
	fin>>A>>B;
	cout<<A<<" "<<B<<endl;
	fin.close();

	//collecting all the odd digit palindromes
	genOdd();

	//collecting all the odd digit palindromes, actually ,we can only add 11,cause others are not prime
	//getEven();
	palindromes.push_back(11);

	Sort(palindromes);

	ofstream fout("pprime.out");
	for(int k = 0; k < palindromes.size(); k++) {
		int temp = palindromes[k];
		if(temp >= A && temp <=B) {
			if(isPrime(temp)) fout<<temp<<endl;
		}
	}
	fout.close();

    return 0;
}

void genOdd() {
	for(int i = 5; i <= 9999; i++) {
		int n = i, j = i / 10;
		while(j) {
			n = n * 10 + j % 10;
			j /= 10;
		}
		palindromes.push_back(n);
	}
}

void genEven() {
	for(int i = 5; i <= 9999; i++) {
		int n = i, j = i;
		while(j) {
			n = n * 10 + j % 10;
			j /= 10;
		}
		palindromes.push_back(n);
	}
}

bool isPrime(int n) {
	if(n % 2 == 0 && n > 2) return false;
	if(n > 5 && n % 5 == 0) return false;

	bool flag = true;

	int start = 2;
	int end = sqrt((double)n);
	for(int i = start; i <= end; i++) {
		int mod = n % i;
		if( mod == 0) {
			flag = false;
			break;
		}
	}
	return flag;
}

void Sort(vector<int> &vec) {
	int n = vec.size();
	if(n <= 1) return;   // already sorted

	vector<int> v1, v2;

	int i = 0;
	while(i < n) {
		if(i < n/2) {
			v1.push_back(vec[i++]);
		} else {
			v2.push_back(vec[i++]);
		}
	}

	Sort(v1);
	Sort(v2);
	vec.clear();
	Merge(vec, v1, v2);
}

void Merge(vector<int> &vec, vector<int> &v1, vector<int> &v2) {
	int n1 = v1.size();
	int n2 = v2.size();

	int p1 = 0;
	int p2 = 0;

	while(p1 < n1 && p2 < n2) {
		if(v1[p1] <= v2[p2]) {
			vec.push_back(v1[p1++]);
		} else {
			vec.push_back(v2[p2++]);
		}
	}

	while(p1 < n1) vec.push_back(v1[p1++]);
	while(p2 < n2) vec.push_back(v2[p2++]);
}