浮子的精度损失

时间:2017-04-30 18:17:10

标签: java floating-point floating-accuracy

我正在尝试计算浮点数的总和。所有小数字输出正确,当我使用非常大的数字作为输入时,输出总是由几个整数关闭。例如, H = 5764801 W = 1679616,在纸面上,计算为335923 30275911 。但是在我的程序中,会打印335923 30275908 。这是代码:

public void printOutput(int H, int W) // The inputs
{
    if(H == 1 && W == 1)
    {
        System.out.println(0 + " " + 1);
        return;
    }

    List<Integer> pfw = primeFactors(W);

    int y = 1;

    while(H != (int) (Math.pow(Math.pow(W,  1f/y) + 1f, y))) y++; 

    final float N = findWholeNumber(pfw);

    float height = 0;
    for(int x = 1; x <= y + 1; x++)
    {
        height += (float) (W * Math.pow((N + 1f) / N, x-1f) + 1e-8); //Here is the summation
    }
    float cats = 1;
    for(int x = 2; x <= y + 1; x++)
        cats += (float) (Math.pow(N, x-1));

    int notWorking = (int) (cats - W);
    System.out.println(notWorking + " " + (int)height); //Outputs printing
}

private int findWholeNumber(List<Integer> factors)
{
    List<Integer> common = new ArrayList<Integer>();
    for(int i = 0; i < factors.size(); i++)
    {
        if(common.contains(factors.get(i))) continue;
        common.add(factors.get(i));
    }
    int num = common.get(0);
    for(int i = 1; i < common.size(); i++)
        num *= common.get(i);
    return num;
}

private List<Integer> primeFactors(int num)
{
    List<Integer> pf = new ArrayList<Integer>();

    if(num == 1)
    {
        pf.add(1);
        return pf;
    }

    for(int j = 2; j <= num; j++)
        while(num % j == 0) // is prime
        {
            pf.add(j);
            num /= j;
        }

    return pf;
}

}

2 个答案:

答案 0 :(得分:2)

浮点数的精度有限,因为尾数的宽度有限。

你可以尝试double为你的情况精确度更高(因为它的尾数更宽),但它也是有限的。

更多信息:https://en.wikipedia.org/wiki/IEEE_floating_point#IEEE_754-2008What is the maximum number in the mantissa part of a Java float?

如果您需要无限制的精度,请尝试BigDecimal。有效数字的数量仅受记忆量的限制。

如果您只需要整数值,BigInteger是一个选项。

答案 1 :(得分:1)

研究每个计算机科学家应该知道的关于浮点运算的知识,David Goldberg,1991。 https://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html