用C语言求圆周率，要求用到for循环。

【C++编译器中的运算程序】
微机WindowsXP中Dev-cpp中的运算程序（30000位）（C++）
#include <cstdlib>
#include <iostream>
#include <fstream>
#define N 30015
using namespace std;
void mult (int *a,int b,int *s)
{
for (int i=N,c=0;i>=0;i--)
{
int y=(*(a+i))*b+c;
c=y/10;
*(s+i)=y%10;
}
}
void divi (int *a,int b,int *s)
{
for (int i=0,c=0;i<=N;i++)
{
int y=(*(a+i))+c*10;
c=y%b;
*(s+i)=y/b;
}
}
void incr(int *a,int *b,int *s)
{
for (int i=N,c=0;i>=0;i--)
{
int y=(*(a+i))+(*(b+i))+c;
c=y/10;
*(s+i)=y%10;
}
}
bool eqs(int *a,int *b)
{
int i=0;
while (((*(a+i))==(*(b+i)))&&(i<=N)) i++;
return i>N;
}
int main(int argc, char *argv[])
{
cout << "正在计算 . . . (0%)";
int lpi[N+1],lls[N+1],lsl[N+1],lp[N+1];
int *pi=lpi,*ls=lls,*sl=lsl,*p=lp;
for (int i=0;i<=N;i++)*(pi+i)=*(ls+i)=*(sl+i)=*(p+i)=0;
memset(pi,0,sizeof(pi));
memset(ls,0,sizeof(ls));
memset(sl,0,sizeof(sl));
memset(p,0,sizeof(p));
*pi=*ls=*sl=1;
for (int i=1;true;i++)
{
mult(ls,i,sl);
divi(sl,2*i+1,ls);
incr(pi,ls,p);
if (eqs(pi,p))
{
cout << "\b\b\b\b100%)\n";
break;
}
int *t;
t=p;
p=pi;
pi=t;
//if (i%1000==0) cout << i << " ";
if(i%1000 == 0)
{
/*cout << i/1000 << "% ";
if(i%5000 == 0)
cout << endl;*/
if(i/1000 < 11)
{
cout << "\b\b\b";
} else {
cout << "\b\b\b\b";
}
cout << i/1000 << "%)";
}
}
cout << endl;
cout << "计算完成\n正在保存 . . .\n";
mult(p,2,pi);
ofstream fout("pi.txt");
fout << *pi << ".";
for (int i=1;i <= N - 15;i++)
{
fout << *(pi+i);
if (i%10==0) fout << " ";
if (i%80==0) fout << endl;
}
cout << "保存完成\n";
cout<< "按回车键退出";
cin.peek();
return EXIT_SUCCESS;
}
注：①运行时会有数据弹出，这无关紧要，只为了加快了感觉速度；
注：程序中有语法错误。请高人改正。
运行环境 CodeBlocks C++
#include <iostream>
using namespace std;
long long a=1000000, b, c=2800000, d, e, f[2801000], g;
int main()
{
for( ;b-c; ) f[b++] =a/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 ) ;
return 0;
}

参考资料：http://baike.baidu.com/view/3287.html?wtp=tt

温馨提示：答案为网友推荐，仅供参考

当前网址：https://verywind.cn/ee/vxjx2223x.html