为什么powx,1 / p和powx,1.0 / p不相等,即使打印它们的值会得到相同的结果

问题描述

问题是:

您将获得2个数字(N,M);任务是找到N√M(M的第N个根)。

输入:

输入的第一行包含一个整数T,它表示测试用例的数量。然后是T测试用例。每个测试用例包含两个以空格分隔的整数N和M。

输出

对于每个测试用例,在新行中打印一个表示M的第N个根的整数(如果该根是整数,则打印-1)。

现在我对这个问题的解决方案是:

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

int main() {
    int t;
    float x,p;
    cin>>t;
    for(int i=0;i<t;i++)
    {
        cin>>p>>x;
        if(p==0)
        {
            cout<<"1"<<endl;

        }
        
        else
        {
            float res=pow(x,(1/p));
            cout<<"res="<<res<<endl;
            if(res==int(res))
                cout<<res<<endl;
            else
                cout<<"-1"<<endl;
        }    
    }
    return 0;
}

这在测试用例中引起了问题:

1
3
1000

尽管我打印res时得到的结果是10,但在条件检查if(res==int(res))期间却是假的。

我还注意到从float res=pow(x,(1/p));float res=pow(x,(1.0/p));的变化给出了正确的答案。我猜测它与完成0.33333的评估时获取1/p有关,但是我不明白为什么打印的值是10但在条件检查中不匹配。>

解决方法

为什么pow(x,1 / p)和pow(x,1.0 / p)不相等,即使打印它们的值得到相同的结果(?)

计算因精度而异。计算在数学上并不精确:pow()未收到精确的 1/3。

对于float p1/p1.0/p可能不同,因为第一个使用float(或更宽泛的数学)完成,第二个使用double(或较宽)为1.0是double

依次调用floatdouble pow()。因此,可能会有不同的res结果。

在OP的情况下:pow(1000.0f,(1/3.0f))执行float精度计算,类似于pow(1000,near_one_third)而不是cbrt(1000.0)-结果不完全是10.0fpow(1000.0f,(1.0/3.0f))执行了类似的double精度计算,当四舍五入为float时,精确计算为10.0f

为什么打印值是10但在条件检查中不匹配。

如果res的计算值大于或小于10.0f,那么==是不正确的。

以足够的精度打印,以查看最终输出中的差异。推荐float res的至少9位有效数字。


至少,我建议始终使用相同的浮点类型和数学运算。 (将1.0f与float一起使用。)进一步建议使用double

最终输出仍可能不是pow()所精确的期望整数(来自数学分析)。