问题描述
问题是:
您将获得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 p
,1/p
与1.0/p
可能不同,因为第一个使用float
(或更宽泛的数学)完成,第二个使用double
(或较宽)为1.0是double
。
依次调用float
或double
pow()
。因此,可能会有不同的res
结果。
在OP的情况下:pow(1000.0f,(1/3.0f))
执行float
精度计算,类似于pow(1000,near_one_third)
而不是cbrt(1000.0)
-结果不完全是10.0f
。 pow(1000.0f,(1.0/3.0f))
执行了类似的double
精度计算,当四舍五入为float
时,精确计算为10.0f
。
为什么打印值是10但在条件检查中不匹配。
如果res
的计算值大于或小于10.0f
,那么==
是不正确的。
以足够的精度打印,以查看最终输出中的差异。推荐float res
的至少9位有效数字。
至少,我建议始终使用相同的浮点类型和数学运算。 (将1.0f与float
一起使用。)进一步建议使用double
。
最终输出仍可能不是pow()
所精确的期望整数(来自数学分析)。