为什么 pow(x,1/p) 和 pow(x,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));给出了正确的答案。我猜这与完成1/p评估时获得0.33333有关,但我不明白为什么打印值10但在条件检查中不匹配。

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

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

对于float p1/p可以与1.0/p不同,因为第一个是使用float(或更宽)的数学来完成的,第二个使用double(或更宽)的,因为1.0是一个double

这反过来又调用floatdoublepow()。因此,可能存在不同的res结果。

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

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

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

以足够的精度打印以查看最终输出的差异。 建议至少使用 9 位有效数字float res


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

不过,最终输出可能仍然不是确切的预期整数(来自数学分析),pow()

最新更新