利用公式π/4≈1-1/3+1/5-1/7+1/9-…求解
我的代码如下：
#include <stdio.h>
#include <math.h>
int main(void)
{
int sign=1;//定义sign为符号，sum为总和，deno为分母，term为一项值
double sum=0.0,deno,term;
for(deno=1.0,term=1.0;fabs(term)>=1e-6;deno=deno+2)
{ term = sign/deno;
sum = sum+term;
sign = -sign;
}
printf("%10.8lf\n",4*sum);
return 0;
}
为啥输出结果是3.14159465啊？
补充一点，同样是利用上述公式，这样得到的值是3.14159065
#include <stdio.h>
#include <math.h>
int main(void)
{
int sign=1;//定义sign为符号，sum为总和，deno为分母，term为一项值
double sum=0.0,deno=1,term=1;
for(;fabs(term)>=1e-6;)
{ sum = sum+term;
deno = deno+2;
sign = -sign;
term = sign/deno;
}
printf("%10.8lf\n",4*sum);
return 0;
}
补充：
你比较一下第一个和第二个，第一个里面deno第一项中是1，第二个里面deno第一项是3，这明显不一样啊
浮点型小数点后面有效数字为6~7位
双精度浮点型小数位后有效位数15~16位。
明白了吧
用g++编译器
//圆周率的计算方法!(能精确到800多位)
#include <stdio.h>
const int N=2800,A=10000;
int c=N,e=0,f[N+10];
int main()
{
for (int& i:f) i=A/5;
while (c)
{
int g=c*2,d=0;
for (int i=c;i!=0;i--)
{
d+=f[i]*A;
f[i]=d%(g-1),d/=g-1;
g-=2;
if (i-1) d*=i-1;
}
int digit=e+d/A;
if (digit==3141)
{
printf("Pi=%.3f",digit*1.0/1000);
e=d%A;
continue;
}
else if (digit<10)printf("000");
else if (digit<100) printf("00");
else if (digit<1000) printf("0");
printf("%d",digit);
e=d%A;
c-=14;
}
getchar();
}
#include <stdio.h>
const int N=2800,A=10000;
int c=N,e=0,f[N+10];
int main()
{
for (int& i:f) i=A/5;
while (c)
{
int g=c*2,d=0;
for (int i=c;i!=0;i--)
{
d+=f[i]*A;
f[i]=d%(g-1),d/=g-1;
g-=2;
if (i-1) d*=i-1;
}
int digit=e+d/A;
if (digit==3141)
{
printf("Pi=%.3f",digit*1.0/1000);
e=d%A;
continue;
}
else if (digit<10)printf("000");
else if (digit<100) printf("00");
else if (digit<1000) printf("0");
printf("%d",digit);
e=d%A;
c-=14;
}
getchar();
}本回答被网友采纳
