c++ - 用C++计算Ramanujan常数

Question

可以看到Ramanujan's constant .
这是我的代码:

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

int main ()
{
long double s=sqrt(163);
long double P=M_PI;
long double R=exp(s*P);

cout.precision(150);

cout<<"Pi=  "<<P<<"\n"<<"sqrt(163)=  "<<s<<"\n"<<"R=  "<<R;
return(0);
}   

这是我的输出:

Pi=  3.141592653589793115997963468544185161590576171875
sqrt(163)=  12.767145334803704059822848648764193058013916015625
R=  262537412640768256  

我的程序有什么问题？
正确的输出是262537412640768743.99999999999925...

Answer 1

double 大约有 15 decimal digits of precision .

您的程序使用M_PI，它是一个double，not a long double .

因此，您的答案大约有 15 位精度。