编程计算圆周率的探索

1、引言

偶尔在刷抖音时，突然想到能否利用编程去计算圆周率。

随即立马去查资料，然后看到了利用单位圆+内接正六边形2N等分方式，利用电脑的快速计算能力去分隔单位圆。图形如下所示：

如上图所示：AB为正六边形的一条边，OQ垂直于AB，和外接圆的交点 为Q，与AB交于C。

Len_n初值为单位圆的半径长，也等于内接正六边形的边长。

Len_2n的初值为线段AQ的长度。

单位圆的半径默认为1。

2、原理和步骤

2.1、第一次分隔

从内接正六边形开始，此时没有开始分隔，即PI=周长/直径=Len_n*6/2。

单位圆的半径为1，此时圆周率PI=3。

下面开始第一次分隔，也就是上述图中的AQ=AB=Len_2n，计算Len_2n初值如下：

利用勾股定理：

其中：，，

注意：以上公式都是最简单的勾股定理应用，这里需要强调的是，当第二次对AQ进行同样分割后，上面的推导过程和涉及到的变量是否通用？显然，第二次分割后的总长度不变，再加上圆的半径不变，上述计算过程通用。小伙伴们可以自行推导试试。

最后结果如下：

即

2.2、代码如下

#include "stdio.h"
#include "math.h"
/*
 * 初值 Len_n 为圆 的 正六边形的 边长， 
 * 初值 圆的半径为 1， 则正六边形的边长为 1， 周长为 6， Pi = 6 / 2 = 3
 * Len_2n = sqrt( 2 - sqrt( 4 - Len_n*Len_n )  )
 *
*/
int main(int argc, char* argv[])
{
	double len_n = 1.0;
	double len_2n = 0.0;
	unsigned char r = 1;
	double C_len = 6*len_n*1.0; //正多边形周长
	unsigned long N = 6; //以正六边形为基准

	double m_PI = 0.0f;
	
	unsigned long num = 800; //分割倍数
	unsigned long i = 0;

	for( i = 1; i< num; i=i*2 )
	{
		len_2n = sqrt( 2 - sqrt( 4 - pow( len_n, 2)));
		C_len = len_2n*2*N*i*1.0;
		m_PI = C_len*1.0/2;

		len_n = len_2n;

		printf("%d 边形， PI = %08f, i=%d\n", 2*N*i, m_PI, i);
	}
	
	printf("Hello World!\n");
	return 0;
}

运算结果为：

12 边形， PI = 3.105829, i=1
24 边形， PI = 3.132629, i=2
48 边形， PI = 3.139350, i=4
96 边形， PI = 3.141032, i=8
192 边形， PI = 3.141452, i=16
384 边形， PI = 3.141558, i=32
768 边形， PI = 3.141584, i=64
1536 边形， PI = 3.141590, i=128
3072 边形， PI = 3.141592, i=256
6144 边形， PI = 3.141593, i=512
Hello World!
Press any key to continue

目前的C环境下，小数位数有限，推荐大家在Python环境下测试。

这个时候的你，有没有感受到那个历史上最混乱的时候，竟然诞生了如此了不起的人物，吾辈还需努力吧。