C语言经典算法 - 超长PI计算的代码
下边内容段是关于C语言经典算法 - 超长PI计算的内容，应该能对大伙有一些用。
#include <stdio.h>
#define L 1000
#define N L/4+1
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)(L / 1.39793 + 1);
int k;
for (k = 1; k <= n; k++)
{
div(w, 25, w);
div(v, 239, v);
div(v, 239, v);
sub(w, v, q);
if (k % 2)
add(s, q, s);
else
sub(s, q, s);
}
printf("%d.", s[0]);
for (k = 1; k < N; k++)
printf("%04d", s[k]);
printf("n");
return 0;
}
{
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;
}
}
}
{
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;
}
}
}
{
int i, tmp, remain = 0;
for (i = 0; i <= N + 1; i++)
{
tmp = a[i] + remain;
c[i] = tmp / b;
}
}
