c++ - 计算pi的值-我的代码有什么问题

Question

我正在做另一个 C++ 练习。我必须从无限级数中计算 pi 的值:

pi=4 - 4/3 + 4/5 – 4/7 + 4/9 -4/11+ 。 . .

程序必须在该系列的前 1,000 项中的每一项之后打印 pi 的近似值。这是我的代码:

#include <iostream>
using namespace std;

int main()
{
    double pi=0.0;
    int counter=1;

    for (int i=1;;i+=2)//infinite loop, should "break" when pi=3.14159
    {
        double a=4.0;
        double b=0.0;

        b=a/static_cast<double>(i);

        if(counter%2==0)
            pi-=b;
        else
            pi+=b;

        if(i%1000==0)//should print pi value after 1000 terms,but it doesn't  
            cout<<pi<<endl;

        if(pi==3.14159)//this if statement doesn't work as well
            break;

        counter++;
    }

    return 0;
}

它编译没有错误和警告，但执行后只出现空的控制台窗口。如果我删除行“if(i%1000==0)”，我可以看到它确实运行并打印每个 pi 值，但它不会停止，这意味着第二个 if 语句也不起作用。我不知道还能做什么。我假设这可能是一个简单的逻辑错误。

Answer 1

好吧，i % 1000 永远不会 = 0，因为您的计数器从 i = 1 开始运行，然后以 2 为增量运行。因此，i 始终是奇数，永远不会是 1000 的倍数。

它永不终止的原因是算法没有精确地收敛到 3.14157 - 它将是一个更高的精度，无论是低于还是高于近似值。你想说“当在 3.14157 的给定增量内时”，所以写

if (fabs(pi - 3.14157) < 0.001)
  break

或类似的东西，无论你想在停止前“接近”什么。