本 站 搜 索     推 荐 文 章        More... 汉语是伟大的，英语是愚蠢的 　　汉语是伟大的，英语是愚蠢的 　..     文 章 阅 读 排 行 复杂成常态，现代网站技术就是.. 《复杂成常态，现代网站技术就是个笑.. 百种中华典籍将走向世界 　　由于中国文字的特殊性，世界各国.. 利润四定律 　　所有的企业都应该为利润而努力，.. 中式英语风靡全球 语言贡献彰.. 　　这些年中国贡献的英语词句： .. 圆周率简介及历史 　　【圆周率简介】 　　圆周率是.. 汉字简化，得不偿失 　　几十年来流行的观点认为：汉字的.. 中国电话区号的编制规律 　　中国电话区号的分布是有规律的。.. 一个苹果的沉痛代价 　　河南夏邑县第一高级中学分校高三.. 圆周率的推算 　　梁羽生兄在《数学与逻辑》一文中.. 圆周率的探索者(图) Archimedes (BC287-BC212) 祖冲之 (..

文 章 信 息 圆周率，计算机C++编译器中的运算程序 〖评论(0)〗〖留言〗〖收藏〗 〖文章分类：话题 / 科技〗〖阅读选项〗 微机 Windows XP 中 Dev-cpp 中的运算程序（30000位）（C++） #include <cstdlib.h> #include <iostream.h> #include <fstream.h> #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[]) { 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)) break; int *t; t=p; p=pi; pi=t; if (i%50==0) cout << i << " "; } cout << endl; mult(p,2,pi); ofstream fout("pi.txt"); fout << *pi << "."; for (int i=1;i<=N;i++) { fout << *(pi+i); if (i%10==0) fout << " "; if (i%80==0) fout << endl; } return EXIT_SUCCESS; } 注： ①运行时会有数据弹出，那是无关紧要的，只为了加快了感觉速度； ②最后的txt文本里有30015位，其中最后15位可能是错的。 ③程序中的<，>是大写的请改成小写。 文章作者：未知  更新日期：2008-09-10 〖文章浏览：今天：0，本周：0，本月：0，总计：0〗〖发送文章〗〖打印文章〗 上下篇章 ·圆周率小数点后21500位 2008-09-10 ·圆周率的计算历史 2008-09-10 ·圆周率的推算 2008-09-10 相关文章 ·圆周率简介及历史 2008-09-10 ·圆周率的计算及方法 2008-09-10 ·圆周率小数点后21500位 2008-09-10 ·圆周率的计算历史 2008-09-10 ·圆周率的推算 2008-09-10