【随笔】Π （pi）的表达——C语言

方法一、直接在库函数<math.h>的定义中使用

方法如下

#define _USE_MATH_DEFINES
#include<stdio.h>
#include<math.h>
int main()
{
	printf("%2.16lf\n", M_PI);
	return 0;
}

输出结果

3.1415926535897931

对比定义中的数值

#define M_E        2.71828182845904523536   // e
#define M_LOG2E    1.44269504088896340736   // log2(e)
#define M_LOG10E   0.434294481903251827651  // log10(e)
#define M_LN2      0.693147180559945309417  // ln(2)
#define M_LN10     2.30258509299404568402   // ln(10)
#define M_PI       3.14159265358979323846   // pi
#define M_PI_2     1.57079632679489661923   // pi/2
#define M_PI_4     0.785398163397448309616  // pi/4
#define M_1_PI     0.318309886183790671538  // 1/pi
#define M_2_PI     0.636619772367581343076  // 2/pi
#define M_2_SQRTPI 1.12837916709551257390   // 2/sqrt(pi)
#define M_SQRT2    1.41421356237309504880   // sqrt(2)
#define M_SQRT1_2  0.707106781186547524401  // 1/sqrt(2)

可以看到，从第16位开始是错误的，根据方法二对比，怀疑是double的极限，当使用long double型时也发生错误，所以姑且只用double就能满足大部分要求，如果需要更高位数，建议手搓（还不会）。

方法二、使用三角函数计算

代码如下

#include<stdio.h>
#include<math.h>
int main()
{
	double a = acos(-1);
	double b = atan(1) * 4;
	printf("%2.16lf\n", a);	
	printf("%2.16lf\n", b);
	return 0;
}

输出结果

3.1415926535897931
3.1415926535897931

可以看到，与方法一相同，仍然是在第16位出现错误，但相比方法一，省了一个定义（

方法三、背圆周率

3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679
显然是不可能放进double的，可以考虑放数组，（如果需要的话）

总结

直接使用库函数<math.h>中的三角函数，用acos(-1)代替pi是最为方便的，直接使用定义也不是不行，但感觉那一串定义会更加麻烦。