ACM中的浮点数精度处理

    在ACM中，精度问题很是常见。其中计算几何头疼的地方通常在于代码量大和精度问题，代码量问题只要平时注意积累模板通常就不成问题了。精度问题则很差说，有时候一个精度问题就可能成为一道题的瓶颈，让你debug半天都找不到错误出在哪。

1.浮点数为啥会有精度问题：

浮点数(以C/C++为准)，通常用的较多的是float, double。

	占字节数	数值范围	十进制精度位数
float	4	-3.4e-38～3.4e38	6~7
double	8	-1.7e-308～1.7e308	14~15

若是内存不是很紧张或者精度要求不是很低，通常选用double。14位的精度(是有效数字位，不是小数点后的位数)一般够用了。注意，问题来了，数据精度位数达到了14位，但有些浮点运算的结果精度并达不到这么高，可能准确的结果只有10~12位左右。那低几位呢？天然就是不可预料的数字了。这给咱们带来这样的问题：即便是理论上相同的值，因为是通过不一样的运算过程获得的，他们在低几位有可能(通常来讲都是)是不一样的。这种现象看似没太大的影响，却会一种运算产生致命的影响: ==。恩，就是判断相等。注意，C/C++中浮点数的==须要彻底同样才能返回true。来看下面这个例子:

#include<stdio.h>
#include<math.h>

int main()

{
      double a = asin( sqrt( 2.0 ) / 2 ) * 4.0;
      double b = acos( -1.0 );
      printf( "      a = %.20lf\n", a );
      printf( "      b = %.20lf\n", b );
      printf( " a -  b = %.20lf\n", a - b );
      printf( " a == b = %d\n", a == b );
      return 0;
}

输出：

      a = 3.14159265358979360000

      b = 3.14159265358979310000

      a - b = 0.00000000000000044409

      a == b = 0

咱们解决的办法是引进eps，来辅助判断浮点数的相等。

2. eps

       eps缩写自epsilon，表示一个小量，但这个小量又要确保远大于浮点运算结果的不肯定量。eps最多见的取值是1e-8左右。引入eps后，咱们判断两浮点数a、b相等的方式以下:

定义三出口函数以下: int sgn(double a){return a < -eps ? -1 : a < eps ? 0 : 1;}

则各类判断大小的运算都应作以下修正:

传统意义	修正写法1	修正写法2
a == b	sgn(a - b) == 0	fabs(a – b) < eps
a != b	sgn(a - b) != 0	fabs(a – b) > eps
a < b	sgn(a - b) < 0	a – b < -eps
a <= b	sgn(a - b) <= 0	a – b < eps
a > b	sgn(a - b) > 0	a – b > eps
a >= b	sgn(a - b) >= 0	a – b > -eps

这样，咱们才能把相差很是近的浮点数判为相等;同时把确实相差较大(差值大于eps)的数判为不相等。

PS: 养成好习惯，尽可能不要再对浮点数作==判断。例如，个人修正写法2里就没有出现==。

 

3. eps带来的函数越界

若是sqrt(a), asin(a), acos(a) 中的a是你本身算出来并传进来的，那就得当心了。

若是a原本应该是0的，因为浮点偏差，可能实际是一个绝对值很小的负数(好比1e-12),这样sqrt(a)应得0的，直接因a不在定义域而出错。

相似地，若是a原本应该是±1,则asin(a)、acos(a)也有可能出错。

所以，对于此种函数，必需事先对a进行校订。

4. 输出陷阱I

这一节和下一节同样，都是由于题目要求输出浮点数，致使的问题。并且都和四舍五入有关。

说到四舍五入，就再扯一下相关内容,据我所知有三种常见的方法:

1. printf(“%.3lf”, a);  //保留a的三位小数，按照第四位四舍五入

2. (int)a;  //将a靠进0取整

3. ceil(a); floor(a);   //顾名思义，向上取证、向下取整。须要注意的是，这两个函数都返回double，而非int

其中第一种很常见于输出(nonsense…)。

如今考虑一种状况,题目要求输出保留两位小数。有个case的正确答案的精确值是0.005,按理应该输出0.01,但你的结果多是0.005000000001(恭喜)，也有多是0.004999999999(悲剧),若是按照printf(“%.2lf”, a)输出，那你的遭遇将和括号里的字相同。

解决办法是，若是a为正，则输出a+eps, 不然输出a-eps

典型案例: POJ2826

5. 输出陷阱II

ICPC题目输出有个不成文的规定(有时也成文)，不要输出: -0.000

那咱们首先要弄清，何时按printf(“%.3lf\n”, a)输出会出现这个结果。

直接给出结果好了：a∈(-0.000499999……, -0.000……1)

因此，若是你发现a落在这个范围内，请直接输出0.000。更保险的作法是用sprintf直接判断输出结果是否是-0.000再予处理。

典型案例:UVA746

 

6. 范围越界

这个严格来讲不属于精度范畴了，不过凑数仍是能够的。请注意，虽然double能够表示的数的范围很大，却不是不穷大，上面说过最大是1e308。因此有些时候你得当心了，好比作连乘的时候，必要的时候要换成对数的和。

典型案例:HDU3558

7. 关于set<T>

有时候咱们可能会有这种需求，对浮点数进行 插入、查询是否插入过 的操做。手写hash表是一个方法(hash函数同样要当心设计)，但set不是更方便吗。但set好像是按==来判重的呀？貌似行不通呢。经观察,set不是经过==来判断相等的，是经过<来进行的，具体说来，只要a<b 和 b<a 都不成立，就认为a和b相等，能够发现，

若是将小于定义成:      bool operator < (const Dat dat)const{return val < dat.val - eps;}就能够解决问题了。 (基本类型不能重载运算符，因此封装了下)

8. 输入值波动过大

这种状况不常见，不过能够帮助你更熟悉eps。假如一道题输入说，给一个浮点数a, 1e-20 < a < 1e20。那你还敢用1e-8作eps么？合理的作法是把eps按照输入规模缩放到合适大小。

典型案例: HUSTOJ 1361

9. 一些建议

容易产生较大浮点偏差的函数有asin、 acos。欢迎尽可能使用atan2。

另外，若是数据明确说明是整数，并且范围不大的话，使用int或者long long代替double都是极佳选择，这样就不存在浮点偏差了