循环 计算圆周率 新手求帮忙

[cwl353253492]

循环-02. 计算圆周率(15)

根据下面关系式，求圆周率的值，直到最后一项的值小于给定阈值。


输入格式：

输入在一行中给出小于1的阈值。

输出格式：

在一行中输出满足阈值条件的近似圆周率，输出到小数点后6位。

输入样例：
0.01
输出样例：
3.132157
#include <stdio.h>
int main()
{ int a,i,c,k;
float x,p,t;
scanf("%f",&x);
a=1;
c=1;
p=1;
i=0;
k=1;
t=0;
while(t>=x)
{
//p是总和。依次加上后来的t新项
i++; //i是上面的递增数
a=a*i; //分子依次兼程

k=k+2; //k是下面的递增数
c=c*k;
t=(1.0*a)/c; //后面的一项来与阈值比较
p=p+t; //p是总和。依次加上后来的t新项

}


printf("%.6f",2*p);



return 0;
}

[「已注销」]

#include<stdio.h>
int main()
{
	double a[50];
	double f(int, double *);
	double x,pi;
	int n = 0;
	pi = 0.0;
	scanf("%f",&x);
	do
	{
		n++;
		a[n] = f(n,a);
		pi = pi +a[n];
	}while(a[n] > x);
	printf("%f", 2 * pi);
}

double f(int n, double * a)
{
	if(n == 1) return 1.0;
	return a[n -1] * ( n - 1) / (2 * n - 1);
}

谁来调试这个烂摊子

[百曉生]

楼主想问什么

[一个广电工程师]

楼主要问什么？程序编译不过?运行结果不对？

[emailtome]

Bigint ? Python for the task. Perl for the task. OpenSSL for the task. GNU MP (multi-precision) library

[赵4老师]

仅供参考

#include <stdio.h>
long a=10000;
long b;
long c=2800;
long d;
long e;
long f[2801];
long g;
int main() {
    for(;b-c;) f[b++]=a/5;
//  while (1) {
//      if (0==b-c) break;
//      f[b]=a/5;
//      b++;
//  }

    //f[0 - 2800] = 10000/5
    for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a)
          for(b=c;d+=f[b]*a,f[b]=d%--g,d/=g--,--b;d*=b);
//  while (1) {
//      d=0;
//      g=c*2;
//      if (0==g) break;
//      c-=14;
//      printf("%.4d",e+d/a);
//      e=d%a;
//      b=c;
//      while (1) {
//          d+=f[b]*a;
//          f[b]=d%--g;
//          d/=g--;
//          --b;
//          if (0==b) break;
//          d*=b;
//      }
//  }

    return 0;
}