编程之美 - 浮点数的精确表示

问题描述: 如何将浮点数(小数)转换为分数的形式,这样可以更精确的表达出浮点数的值。

问题分析: 小数可以分为两种,有限小数和无限循环小数。

有限小数的处理相对简单: 只要把它乘以10的N次方,转换为整数,再和10的N次方都除以它们的最大公约数即可。 例如:0.25  把它乘以 10的2次方 变成25, 25/100 分子分母除以最大公约数25变成1/4即可。

无限循环小数的处理要复杂一些: 书中使用了一个有趣的计算方法,它把无限循环小数分为一小数部分和一个循环部分 例如: X=0.a1a2…am(b1b2…bn) 这里a是小数部分,b是循环部分。 1)  X先乘以10的m次方,把小数部分变为整数   X*10^m ==> a1a2…am + 0
.b1b2…bn(b1b2…bn)      X = (a1a2…am + 0
.b1b2…bn(b1b2…bn))/(10^m)

2)  假设 Y = 0.b1b2…bn(b1b2…bn),那么Y可以转换为 Y*10^n = b1b2…bn + 0.b1b2…bn(b1b2…bn)       Y*10^n = b1b2…bn +
 Y   ==>   Y =  b1b2…bn/(10^n – 1)

3)  把Y的部分代入X的等式中  X = (a1a2…am +
 b1b2…bn/(10^n – 1))/(10^m)      X = (a1a2…am * (10^n – 1) + b1b2…bn) / ((10^n – 1)*10^m)

例如 0.3(3) = (3*9 + 3) / (9*10) = 30 / 90 = 1/3



实例程序: 程序中的缺陷,如果循环小数位数太多,会造成乘积超过int的范围,造成计算错误。

#include <iostream>

using namespace std;

typedef struct _data
{
    unsigned long _a;
    unsigned long _b;   // the loop part, if = 0 means no
} data_t;

unsigned long gcb(unsigned long a, unsigned long b)
{
    unsigned long temp = 0, r = 0;

    // make sure a < b
    if (a > b)
    {
        temp = a;
        a = b;
        b = temp;
    }

    r = b % a;
    while(r != 0)
    {
        b = a;
        a = r;
        r = b%a;
    }

    return a;
}

void calc(data_t data)
{
    unsigned long temp = 0, gcb_val = 0;
    unsigned long power_a = 1, power_b = 1;

    temp = data._a;
    while (temp != 0)
    {
        temp = temp/10;
        power_a *= 10;
    }

    temp = data._b;
    while (temp != 0)
    {
        temp = temp/10;
        power_b *= 10;
    }

    cout << "a->" << data._a << ":" << power_a << "   b->" <<data._b << ":" << power_b << endl;

    // no loop part
    if (data._b == 0)
    {
        gcb_val = gcb(data._a, power_a);
        cout << data._a/gcb_val << "/" << power_a/gcb_val << endl;
    }
    else
    {
        temp = (data._a * (power_b - 1)) + data._b;
        power_b = power_a*(power_b - 1);
        gcb_val = gcb(temp, power_b);
        cout << temp/gcb_val << "/" << power_b/gcb_val << endl;
    }
}


int main()
{
    data_t test;
    int val = 0;
    //test._a = 1; test._b = 0;
    //test._a = 23; test._b = 0;
    //test._a = 3; test._b = 3;
    test._a = 1; test._b = 1;
    //test._a = 285714; test._b = 285714;   // error: int is overflow

    calc(test);

    cin >> val;
    return 0;
}

    原文作者:zy__
    原文地址: https://blog.csdn.net/wangzhiyu1980/article/details/51155354
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞