3

我应用以下方法将 Double/Float 形式转换为有理形式。

对于前。

a=0.125
Find t=1000 
Find k=GCD(t,a*t);
Print (a*t/k) "/" t/k

它为 0.125 提供了正确的 o/p，但对于 0.12 则没有。问题在于找到t。请帮我解决这个问题。这是代码，

#include<iostream>
#include<stdio.h>
#include<stdlib.h>
using namespace std;
int gcd(int a,int b)
{
  if(b==0)
  return a;
  else
  return gcd(b,a%b);
}
void DoubleToRational(float a)
{
  int t=1;
  while((float)t*a!=(int)(t*a))
  {
   //cout<<t*a<<" "<<(float)t*a<<" "<<(int)(t*a)<<endl; // For Checking
   t=t*10;
  }
  int k=gcd(t*a,t);
  cout<<(int)t*a/k<<" / "<<t/k;
}
int main()
{
  DoubleToRational(0.125); //0.021 seconds Works Perfectly
 // DoubleToRational(0.12); //Didnt Work ????
  return 0;
}

另外，我不希望这种方法在精度很高时有效。在这种情况下，请建议修改/其他方法。

4

1 回答 1

4

该数字0.125可以用二进制浮点数精确表示。但数字0.12不能。所以你正在寻找最接近的可表示数字的有理值0.12。绝对不是3/25。事实上，最接近的双精度值是0.12：

0.11999 99999 99999 99555 91079 01499 37383 83054 73327 63671 875

要在这项任务中取得任何合理进展，您可能需要使用十进制数据类型而不是二进制类型（如double.

如果以上内容对您没有任何实际意义，那么我建议您好好阅读：每个计算机科学家都应该知道的关于浮点运算的知识。

于 2012-12-30T20:07:07.560 回答