计算PI值

计算PI值

1. 一个PI值计算的算法。
2. 分析该算法。

3. 实现该算法。

   算法：PI/2=1/3 + 1/3 * 2/5 + 1/3 * 2/5 * 3/7 + …….+ 1/3* ……* n/(n+1)

   分析：如果直接写成程序的花很难，既要算乘法又要算除法还要算加法还要保证精度。
   所以只有化简公式，网上有大量的关于这个公式的化简和写成算法的教程，所以我的分析不是要分析如何化简和如何把简化的公式写成算法，我的分析主要是解决，网上博客没有解决的问题，一个就是精度控制，网上大多的程序只是给出了精确到800位的算法，也就是当n=2800的时候。另一个就是算法的时间复杂度和空间复杂的分析，随着n的增加该算法的耗时如何增长，占用的内存空间有时如何增长。

算法实现如下：
1 #define M 2800
2 #include <stdio.h>
3 int main(){
4 long yushu[M+1], k = 0, i = 0, j = 0, sum = 0, PI = 0;
5
6 for (k = 0; k <= M; k++)
7 yushu[k] = 1000 * 2;
8
9
10 for (i = 1; i <= 800 / 4; i++){
11
12 for (j = M, sum = 0; j > 0; j--){
13 sum += yushu[j] * 10000;
14 yushu[j] = sum % (2 * j + 1);
15 sum /= (2 * j + 1);
16 sum *= j;
17 }
18
19
20 sum += yushu[j] * 10000;
21 yushu[j] = sum % (2 * j + 1);
22 printf("%.4ld", PI + sum / 10000);