Double / Float TO Rational(X / Y形式)

时间:2012-12-30 20:05:08

标签: c++

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

对于Ex。

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;
}

当精度很高时,我不希望这种方法有效地工作。 在这种情况下,请建议修改/其他方法。

1 个答案:

答案 0 :(得分:4)

数字0.125可以用二进制浮点精确表示。但数字0.12不能。因此,您要找到最接近的可表示数字的合理值0.12。绝对不是3/25。事实上,closest double precision value0.12是:

0.11999 99999 99999 99555 91079 01499 37383 83054 73327 63671 875

要使用此任务取得任何合理的进展,您可能需要使用十进制数据类型而不是像double这样的二进制类型来表示您的值。

如果上述内容对您没有任何意义,那么我建议您好好阅读:What Every Computer Scientist Should Know About Floating-Point Arithmetic