c++ - 点和三角形

Question

任务 - 点在哪里，在三角形内部或外部。我检查了小三角形和大三角形的 Heron 公式。然后比较它们的总和。很可能是功能“区域”中的错误。

 #include <iostream>
    #include <math.h>

    using namespace std;

    struct point
    {
        int x;
        int y;
    };

    double length(point first, point second)
    {
        double katet1=second.x-first.x;
        double katet2=second.y-first.y;
        return(sqrt(katet1*katet1+katet2*katet2));
    }

    double area(point a, point b, point c)
    {
        double ab = length(a,b);
        double bc = length(b,c);
        double ca = length(c,a);
        double p =(ab+bc+ca)/2;
        double s = sqrt(p*(p-ab)*(p-bc)*(p-ca));
        return(s);
    }

    int main()
    {
        int num;
        cin>>num;
        int win=0;
        for(int i=0; i<num; i++)
        {
            point d,a,b,c;
            cin>>d.x>>d.y>>a.x>>a.y>>b.x>>b.y>>c.x>>c.y;
            double s1,s2,s3,s4;
            s1=area(d,a,b);
            s2=area(d,b,c);
            s3=area(d,a,c);
            s4=area(a,b,c);
            if((s1+s2+s3)==s4)
                        win++;
        }
        cout<<win;
        cin.get();
    }

未通过所有测试。例如 test 1 1 1 0 0 2 2 0 3 必须返回 1 但返回 0。有什么问题？对不起我的英语不好。

Answer 1

您正在尝试将 double 与 operator== 进行比较。
计算可能会有所不同，您需要使用 epsilon 值作为误差范围：

const double Epsilon = 0.0001;
if (((s1+s2+s3) >= s4 - Epsilon) && ((s1+s2+s3) <= s4 + Epsilon))
{
}

Answer 2

您不应该期望在计算机世界(s1+s2+s3) == s4中工作正常。double/float尝试：

abs(s1 + s2 + s3 - s4) < Epsilon

其中 Epsilon 是您所需的精度，例如将其设置为0.00001、0.01、1.0，甚至比10.0基于您的应用程序的精度更高。

--

“每个计算机科学家都应该知道的关于浮点运算的知识”