谁能用C/C++编一个计算π小数点后无限位的程序（意思就是边计算，边显示。）具体见补充

再在显示的最后一行显示时间和小数点后尾数。

小弟有的财富都给你。

举报该文章

相关建议 2013-05-09

#include <stdio.h>
#define MAXLEN 1000
#define N MAXLEN/4 + 1
// MAXLEN 为位数，N是array长度
void add(int*, int*, int*);
void sub(int*, int*, int*);
void div(int*, int, int*);
int main(void) {
int s[N + 3] = {0};
int w[N + 3] = {0};
int v[N + 3] = {0};
int q[N + 3] = {0};
int n = (int)(MAXLEN/1.39793 + 1);
int k,count = 1;
w[0] = 16*5;
v[0] = 4*239;
for(k = 1; k <= n; k++) {
// 套用公式
div(w, 25, w);
div(v, 239, v);
div(v, 239, v);
sub(w, v, q);
div(q, 2*k-1, q);
if(k%2) // 奇数项
add(s, q, s);
else // 偶数项
sub(s, q, s);
}
printf("%d.", s[0]);
for(k = 1; k < N; k++) {
if((count % 10) == 0) printf("\n");
printf("%04d", s[k]);
count++;
}
printf("\n");
return 0;
}
void add(int *a, int *b, int *c) {
int i, carry = 0;
for(i = N + 1; i >= 0; i--) {
c[i] = a[i] + b[i] + carry;
if(c[i] < 10000)
carry = 0;
else { // 进位
c[i] = c[i] - 10000;
carry = 1;
}
}
}
void sub(int *a, int *b, int *c) {
int i, borrow = 0;
for(i = N+1; i >= 0; i--) {
c[i] = a[i] - b[i] - borrow;
if(c[i] >= 0)
borrow = 0;
else { // 借位
c[i] = c[i] + 10000;
borrow = 1;
}
}
}
void div(int *a, int b, int *c) { // b 为除数
int i, tmp, remain = 0;
for(i = 0; i <= N+1; i++) {
tmp = a[i] + remain;
c[i] = tmp / b;
remain = (tmp % b) * 10000;
}
}

追问

如何增加计算位数？？
最大多少位？？
最后直接退出加getchar()不顶用呀

追答

MAXLEN就是给定的位数，下面是在VC++ 6,0环境下的运行结果。
3.141592653589793238462643383279502884
1971693993751058209749445923078164062862
0899862803482534211706798214808651328230
6647093844609550582231725359408128481117
4502841027019385211055596446229489549303
8196442881097566593344612847564823378678
3165271201909145648566923460348610454326
6482133936072602491412737245870066063155
8817488152092096282925409171536436789259
0360011330530548820466521384146951941511
6094330572703657595919530921861173819326
1179310511854807446237996274956735188575
2724891227938183011949129833673362440656
6430860213949463952247371907021798609437
0277053921717629317675238467481846766940
5132000568127145263560827785771342757789
6091736371787214684409012249534301465495
8537105079227968925892354201995611212902
1960864034418159813629774771309960518707
2113499999983729780499510597317328160963
1859502445945534690830264252230825334468
5035261931188171010003137838752886587533
2083814206171776691473035982534904287554
6873115956286388235378759375195778185778
0532171226806613001927876611195909216420
1989
Press any key to continue

追问

我能运行
但把MAXLEN加大后显示“int s[N + 3] = {0};”有错
我不是要结果

追答

不能说说是什么样的错吗？我这儿把MAXLEN增加到35000仍可以正常运行，更大的没测试。

温馨提示：内容为网友见解，仅供参考

当前网址：https://11.t2y.org/zz/2s72mff42.html

其他看法

第1个回答 2013-05-09

计算机的计算精度是有限的，所有编出的程序也是只能算到有限位，现在能算出的也就小数点后几百万位。追问

也可以呀

第2个回答 2013-05-09

PI能计算么

第3个回答 2013-05-09

如果你要很多位的话可以用matlab来求。

求一个用C++编过计算器的,就是那种加减乘除三角函数可以写一排算的
='

C语言\/C++用投针法计算圆周率
代码如下：include <iostream> #include <fstream>#include <math.h>using namespace std;double randx,randnum;unsigned long randi=1;const unsigned long randa=16807;unsigned long randm=pow(2,31)-1;const double pi=3.1415926;\/\/ random number creationvoid random_number(){ if ( randi =...

C语言程序设计,用级数计算圆周率π
include<stdio.h>int main(){ double pi=2.0; double i=1.0; double a= i; double b= 2*i+1; double t=2*a\/b; while(t>=1e-100) { pi = pi + t; i=i+1; a=a*i; b=b*(2*i+1); t=2*a\/b; } printf("%g\\n",pi)...

π是如何算出来的?
用公式编程序计算，程序算到的无穷级数项数越多，得到的精确位数越多。计算的速度和计算出的位数随着公式的改进会得到提高。比较著名的公式有：pi\/4=1-1\/3+1\/5-1\/7+1\/9……pi^2\/6=1+1\/4+1\/9+1\/16……更先进的迭代法之类的可以看看 http:\/\/elephant.linux.net.cn\/articles\/pi.php ...

用函数重载方法编写求面积程序,要求能计算圆,矩形三角形面积_百度知 ...
double dbWide){ return dbLong*dbWide;}\/\/ 计算三角形面积，a、b、c分别为三角形的三条边double getArea(double a, double b, double c){ double p = (a + b + c) \/ 2.0; return sqrt(p*(p - a)*(p - b)*(p - c));}int main(){ double a, b, r, m,...

c++编写程序精确计算100!的值
void Compute(int n){ \/\/一个int值n乘以一个数组 \/\/数组的第0位表示大数的最低位 \/\/思路： \/\/先用最低位乘以n \/\/最低位的进位肯定是0，所以result[j] += aboveTen;没有问题 \/\/十位数以上的值都是应该加到下一位的，暂时存储在变量aboveTen里面 aboveTen = result[j] \/10; \/\/然后将...

...π的值,精确到小数点后6位,结果老是不对。C++程序哪里错了?_百度知 ...
include<cmath> include<string> using namespace std;int main(){ double sum=0.0,item=1.0;for(int n=1;abs(item)>1e-6;++n){ item *=(-1.0)*(2.0*n-3.0)\/(2.0*n-1.0);sum+=item;} cout<<"PI="<<setiosflags(ios::fixed)<<setprecision(6)<<sum*4<<endl;double ...

C语言编程题
\/\/输入15个学生的信息void insertTest(STU *stu);\/\/测试,用随机数代替输入，我懒得输入那么多成绩，你要手输，调用insertSTU替换该函数void printfSTU(STU *stu);\/\/按输入顺序打印15名学生7模块成绩及每个模块前三名的成绩int main(){ STU stu[stuNum];\/\/定义15个学生 \/\/ int i,...

我的C++程序怎么没法运行?(AGM公式算π)
将声明部分放到main外面也就是声明成全局变量。因为一个很大的数组声明在main里容易分配失败

请推导一个圆周率兀的计算公式,谢啦!
fft.cpp FFT 的 C++ 源程序利用 FFT，要计算 n1 位和 n2 位的两个多位数乘法，可以这样进行：开辟两个长度为 n(n>=n1+n2，取 2m 最佳) 的复数数组，将两个多位数从低位到高位分别填入，高位补 0。对两个数组分别进行正向傅立叶变换。将得到的两个变换后的数组的对应项相乘，然后进行反向...

相似回答

大家正在搜