Question

我在C ++（CodeBlocks）中做了一些事情，但我发现了一个奇怪的问题。我将我的代码发送给了我的朋友（他在DevC ++中测试过）并且它有效。我尝试了这两个代码：

#include <iostream>
#include <math.h>
using namespace std;

int main() //this function works
{
    if (pow(3, 2) + pow(4, 2) == pow(5, 2)) {
        cout << "Works" << endl;
    } else { cout << "Nope" << endl; } 
}

但是，然后我改变了这样的主要功能（它没有用）：

int main() //this function doesn't work
{
    int t1 = 3, t2 = 4, t3 = 5;
    if (pow(t1, 2) + pow(t2, 2) == pow(t3, 2)) {
        cout << "Works" << endl;
    } else { cout << "Doesn't work" << endl; }
}

任何人都知道这是什么问题？

Answer 1

除非你告诉我们你的“怪异”错误是什么，我假设你的第二个代码段打印出来：

不起作用

问题是，您要比较浮点数，因为pow()会返回floating points，see the definition of pow() function。

由于Rounding Errors，浮点数学不准确。像9.0这样的简单值无法使用二进制浮点数精确表示，浮点数的有限精度意味着操作顺序的微小变化可能会改变结果。不同的编译器和CPU架构以不同的精度存储临时结果，因此结果将根据您的环境的详细信息而有所不同。例如：

    float a = 9.0 + 16.0
    float b = 25.0
    if(a == b) // can be false!
    if(a >= b) // can also be false!

甚至

    if( Math.abs(a-b) < 0.00001) // wrong - don't do this

这是一个糟糕的方法，因为选择固定的epsilon（0.00001）因为它“看起来很小”实际上可能太大了，因为被比较的数字也非常小。

我个人使用以下方法，

public static boolean nearlyEqual(float a, float b, float epsilon) {
        final float absA = Math.abs(a);
        final float absB = Math.abs(b);
        final float diff = Math.abs(a - b);

        if (a == b) { // shortcut, handles infinities
            return true;
        } else if (a == 0 || b == 0 || diff < Float.MIN_NORMAL) {
            // a or b is zero or both are extremely close to it
            // relative error is less meaningful here
            return diff < (epsilon * Float.MIN_NORMAL);
        } else { // use relative error
            return diff / Math.min((absA + absB), Float.MAX_VALUE) < epsilon;
        }
    }

不要忘记阅读What Every Computer Scientist Should Know About Floating-Point Arithmetic！

参考： This是我的回答的参考。

修改：由于OP的问题与C++有关，所以这里是nearlyEqual()的编辑版本：

#include <iostream>     // std::cout
#include <cmath>        // std::abs
#include <algorithm>    // std::min
using namespace std;

#define MIN_NORMAL 1.17549435E-38f
#define MAX_VALUE 3.4028235E38f

bool nearlyEqual(float a, float b, float epsilon) {
    float absA = std::abs(a);
    float absB = std::abs(b);
    float diff = std::abs(a - b);

    if (a == b) {
        return true;
    } else if (a == 0 || b == 0 || diff < MIN_NORMAL) {
        return diff < (epsilon * MIN_NORMAL);
    } else {
        return diff / std::min(absA + absB, MAX_VALUE) < epsilon;
    }
}

int main(void) {
    float t1 = 3.0, t2 = 4.0, t3 = 5.0, epsilon = 0.0000000001; // don't use int here!
    if (nearlyEqual((pow(t1, 2) + pow(t2, 2)), pow(t3, 2), epsilon)) {
        cout << "Works" << endl;
    } else {
        cout << "Doesn't work" << endl;
    }
    return 0;
}

输出结果为：

作品

编译器： Cygwin C ++编译器。

Cygwin版本： 1.7.25

Answer 2

变量类型。

double pow (double base     , double exponent);
     float pow (float base      , float exponent);
long double pow (long double base, long double exponent);
     double pow (Type1 base      , Type2 exponent);        // additional overloads

http://www.cplusplus.com/reference/cmath/pow/

试

> int main() {
>     double d1 = 3.0, d2 = 4.0, d3 = 5.0;
>     if(pow(d1,2) + pow(d2,2) == pow(d3,2)) {
>         cout << "Works" << endl;
>     }
>     else {
>         cout << "Nope" << endl;
>     } }

以上在VS 2013中也不起作用。

我通过使用函数

让它工作

bool CheckPyth(double a, double b, double c) {
  double aa = pow(a,2);
  double bb = pow(b,2);
  double cc = pow(c,2);

  if(aa + bb == cc) {
    return(true);
  }
  else {
    return(false);
  }
}

这似乎表明，正如CoryKramer指出的那样，它是一个浮点舍入错误。关于d后缀对双重文字无效，他也是正确的，所以我也改变了代码来纠正它。

Answer 3

它也适用于ideone.com！ http://ideone.com/3XcY4G

    int main() 
{
    int t1 = 3, t2 = 4, t3 = 5;
    if (pow(t1, 2) + pow(t2, 2) == pow(t3, 2)) {
        cout << "Works" << endl;
    } else { cout << "Doesn't work" << endl; }
}

Answer 4

问题是math.h不起作用，所以我必须使用cmath。第二，不是int，而是float或double。

代码：

#include <iostream>
#include <cmath>

using namespace std;

int main()
{

    float t1 = 3, t2 = 4, t3 = 5;

    if (pow(t1, 2) + pow(t2, 2) == pow(t3, 2)) {
        cout << "PT" << endl;
    }
    else {
        cout << pow(t1, 2) + pow(t2, 2) << endl;
        cout << pow(t3, 2) << endl;
    }

}

CodeBlocks C ++ Bug

4 个答案: