计算PI值
- 一个PI值计算的算法。
- 分析该算法。
实现该算法。
算法: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;