V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
spencerqiu
V2EX  ›  C

算法竞赛中判断浮点数是否为整数和浮点数与浮点数是否相等的安全做法?

  •  
  •   spencerqiu · 2015-07-24 23:07:42 +08:00 · 1023 次点击
    这是一个创建于 3191 天前的主题,其中的信息可能已经有所发展或是发生改变。
    在不调用其他各类函数的情况下。
    11 条回复    2015-07-25 17:22:08 +08:00
    razrlele
        1
    razrlele  
       2015-07-24 23:15:31 +08:00
    浮点数的相等可以判断 fabs(a - b) < 1e12

    判断浮点数是否是整数的话就判断 fabs(a - (int)a) < 1e12
    lzhtony
        2
    lzhtony  
       2015-07-24 23:16:24 +08:00
    令eps=10e-5
    然后fabs(a-b) < eps即认为a b相等
    当然eps可以根据需要取值
    lzhtony
        3
    lzhtony  
       2015-07-24 23:18:24 +08:00
    @razrlele 是1e-12吧..话说12位小数已经不精确了吧?
    razrlele
        4
    razrlele  
       2015-07-24 23:25:00 +08:00
    @lzhtony 嗯嗯。。。。是-12.。。。不过好像一般-6左右就差不多够了。。。
    lzhtony
        5
    lzhtony  
       2015-07-24 23:32:27 +08:00
    @razrlele 嗯,一般都没那么小的.
    zhicheng
        6
    zhicheng  
       2015-07-24 23:34:29 +08:00 via Android
    winoros
        7
    winoros  
       2015-07-24 23:42:22 +08:00
    const double eps = blabla; //自己根据情况设一个
    int sgn(double x) {
    return x < -eps ? -1 : x > eps;
    }

    x-y的大小根据sgn(x - y)来看
    husinhu
        8
    husinhu  
       2015-07-25 10:37:36 +08:00 via iPhone
    你可以去看看softfloat.c的实现
    zodiac1111
        9
    zodiac1111  
       2015-07-25 11:52:13 +08:00
    NaN呢?正负无穷呢?

    这个网站可能有点用

    isNaN

    a==a

    http://babbage.cs.qc.cuny.edu/IEEE-754/
    http://www.binaryconvert.com/
    zodiac1111
        10
    zodiac1111  
       2015-07-25 11:56:49 +08:00
    试下这几个case

    float f1=0.0/0.0
    float f2=1.0/0.0
    float f3=-1.0/0.0

    ? f1==f1
    ? f1+1.0==f1
    ? f2+f2>f2
    JamesRuan
        11
    JamesRuan  
       2015-07-25 17:22:08 +08:00
    正确的做法是:
    #include<float.h>

    float f1,f2;
    ...
    ...

    if (((f1 - FLT_EPSILON) < f2) &&
    ((f1 + FLT_EPSILON) > f2)) // equal

    或者

    #include<math.h>
    if(fabsf(f1-f2) < FLT_EPSILON)

    同理,比较double,用DBL_EPSILON,long double用LDBL_EPSILON

    这些都是标准库的做法。
    关于   ·   帮助文档   ·   博客   ·   API   ·   FAQ   ·   我们的愿景   ·   实用小工具   ·   5297 人在线   最高记录 6543   ·     Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 · 29ms · UTC 06:56 · PVG 14:56 · LAX 23:56 · JFK 02:56
    Developed with CodeLauncher
    ♥ Do have faith in what you're doing.