精确表达浮点数 -- 编程之美

以下描述来自<<编程之美>>：

在计算机中，使用float或者double来存储小数是不能得到精确值的。如果你希望得到精确计算结果，最好是用分数形式来表示小数。有限小数或者无限循环小数都可以转化为分数。比如：
0.9 = 9/10
0.333（3）= 1/3（括号中的数字表示是循环节）
当然一个小数可以用好几种分数形式来表示。如：
0.333（3）= 1/3 = 3/9
       给定一个有限小数或者无限循环小数，你能否以分母最小的分数形式来返回这个小数呢？如果输入为循环小数，循环节用括号标记出来。下面是一些可能的输入数据，如0.3、0.30、0.3（000）、0.3333（3333）、……

拿到这样一个问题，我们往往会从最简单的情况入手，因为所有的小数都可以分解成一个整数和一个纯小数之和，不妨只考虑大于0，小于1的纯小数，且暂时不考虑分子和分母的约分，先设法将其表示为分数形式，然后再进行约分。题目中输入的小数，要么为有限小数X=0.a1a2…an，要么为无限循环小数X=0.a1a2…an（b1b2…bm），X表示式中的字母a1a2…an，b1b2…bm都是0~9的数字，括号部分（b1b2…bm）表示循环节，我们需要处理的就是以上两种情况。

对于有限小数X=0.a1a2…an来说，这个问题比较简单，X就等于（a1a2…an）/10n。
对于无限循环小数X=0.a1a2…an（b1b2…bm）来说，其复杂部分在于小数点后同时有非循环部分和循环部分，我们可以做如下的转换：
X= 0.a1a2…an（b1b2…bm）

10n* X= a1a2…an.（b1b2…bm）

10n* X= a1a2…an+0.（b1b2…bm）

X =（a1a2…an+0.（b1b2…bm））/10n

    对于整数部分a1a2…an，不需要做额外处理，只需要把小数部分转化为分数形式再加上这个整数即可。对于后面的无限循环部分，可以采用如下方式进行处理：
令Y=0. b1b2…bm，那么
10m *Y=b1b2…bm.（b1b2…bm）

10m *Y=b1b2…bm+0.（b1b2…bm）

10m *Y-Y=b1b2…bm

Y= b1b2…bm/（10m-1）

将Y代入前面的X的等式可得：

X=（a1a2…an+Y）/10n
=（a1a2…an+ b1b2…bm/（10m-1））/10n
=（（a1a2…an）*（10m-1）+（b1b2…bm））/（（10m-1）*10n）  -----重点
       至此，便可以得到任意一个有限小数或无限循环小数的分数表示，但是此时分母未必是最简的，接下来的任务就是让分母最小，即对分子和分母进行约分，这个相对比较简单。对于任意一个分数A/B，可以简化为（A/Gcd（A,B））/（B/Gcd（A,B）），其中Gcd函数为求A和B的最大公约数，这就涉及本书中的算法（2.7节“最大公约数问题”），其中有很巧妙的解法，请读者阅读具体的章节，这里就不再赘述。
      综上所述，先求得小数的分数表示方式，再对其分子分母进行约分，便能够得到分母最小的分数表现形式。

以下代码未处理大数的情况，所以会溢出！只是为了展现原理

#include <iostream>
#include <cmath>
#include <string>
#include <sstream>

using namespace std;

struct Fraction
{
	int Denominator;
	int numerator;
};

int gcd(int x,int y);
int simple_gcd(int x,int y);
Fraction deal_num(string str);

int main()
{
	cout<<"start"<<endl;
	string str;
	cin>>str;
	Fraction frac;
	frac = deal_num(str);
	cout<<frac.Denominator<<"/"<<frac.numerator<<endl;
	cout<<"end"<<endl;
	return 1;
}
//未考虑数值大溢出问题！！
Fraction deal_num(string str)
{
	int pos_dot;
	int pos_Bracket1;
	int pos_Bracket2;

	pos_dot = str.find('.');
	pos_Bracket1 = str.find('(',pos_dot+1);
	pos_Bracket2 = str.find(')',pos_Bracket1+1);
	
	int len_a,len_b;
	string str_a(str.substr(pos_dot+1,pos_Bracket1-pos_dot-1));
	string str_b(str.substr(pos_Bracket1+1,pos_Bracket2-pos_Bracket1-1));
	len_a = str_a.size();
	len_b = str_b.size();

	unsigned int a;
	unsigned int b;
	stringstream str_stream(str_a);
	str_stream>>a;
	str_stream.clear();
	str_stream.str("");
	str_stream<<str_b; //出问题
	str_stream>>b;
	str_stream.clear();
	
	int temp_a,temp_b;
	temp_a = (int)pow((float)10,len_a); //容易溢出
	temp_b = (int)pow((float)10,len_b);
	Fraction frac;
	frac.Denominator = a*(temp_b - 1)+b;
	frac.numerator = temp_a*(temp_b - 1);//容易溢出
	
	int gcd_val = gcd(frac.Denominator,frac.numerator);
	frac.Denominator /= gcd_val;
	frac.numerator /= gcd_val;

	return frac;
}


int gcd(int x,int y)
{
	if (x < y)
	{
	  return gcd(y,x);//用小的 --或者直接互换
	}
	if(y == 0)return x;

	if ((x&0x1) == 0)//偶数
	{
		if((y&0x1) == 0)//偶数
		{
			return 2*gcd(x>>1,y>>1);
		}
		else//奇数
		{
			return gcd(x>>1,y);
		}
	}
	else//奇数
	{
		if((y&0x1) == 0)
		{
			return gcd(x,y>>1);
		}
		else//奇数
		{
			return gcd(y,x-y);
		}
	}


}

int simple_gcd(int x,int y)
{
   return !y?x:simple_gcd(y,x%y);
}