Question

假设从a [0]到a [14]有15个数字满足：

a[0]=4
a[i]=pow(sqrt(a[i-1])*2-1, 2) for i>=1 and i<=14

我需要打印出从a [0]到a [14]的所有值，每个数字都在一行上，因此我编写了以下代码（C ++，Code :: Blocks）：

#include <iostream>
#include <cmath>
using namespace std;
#define maxN 15
int i[maxN]; int n=15;
int main()
{
    i[0]=2;
    cout<<4<<endl;
    for (int k=1; k<n; k++){
        i[k]=i[k-1]*2-1;
        cout<<pow(i[k],2)<<"\t"<<endl;
    }
    return 0;
}

结果是：

4
9
25
81
289
1089
4225
16641
66049
263169
1.05062e+006
4.1984e+006
1.67854e+007
6.71252e+007
2.68468e+008

最后五个数字不正确（由于整数溢出）。

在上面的“ for”循环中，我更改了行

cout<<pow(i[k],2)<<"\t"<<endl;

到

cout<<(long) pow(i[k],2)<<"\t"<<endl;

这次，结果是：

手动检查后，我意识到许多数字仍然不正确：24； 4224; 66048; 263168； 4198400； 16785408; 67125248； 268468224（它们全都比正确的数字低1）。我该怎么做才能解决此问题？

Answer 1

1.05062e + 006和1050625是相同的值，但打印方式不同。 pow(i[k],2) returns a floating point type，在此不会发生整数溢出，而是在将此类浮点转换为long时会四舍五入。

格式Xe+Y被称为科学记号，其值为X * 10 ^ Y。默认情况下，std::printf和std::ostream::operator<<打印浮点值，其值超出其科学计数法的特定范围。有关更多信息，请参见https://en.cppreference.com/w/cpp/io/manip/fixed。

Answer 2

该如何解决？

更改

cout<<pow(i[k],2)<<"\t"<<endl;

到

cout<<i[k]*i[k]<<"\t"<<endl;

pow是一个浮点函数，可能不精确。请注意，您的情况有些奇怪：如果pow使用64位IEEE-754双精度，则它将为您的输入打印准确的数字。因此，也许您使用32位float数字。

为什么我的代码输出错误的结果？整数溢出了吗？

2 个答案: