Question

我编写了格雷厄姆的算法，但它仍然给出了凸包的错误点。我需要帮助。我认为我的签名功能有一个错误，但不知道它是什么。

#include <cstdio>
#include <algorithm>
#include <math.h>
#define pb push_back
#define mp make_pair
#include <vector>

using namespace std;

vector <pair<double, double> > st;
pair<double, double> p[1000];
double x, y;

int f(pair <double,double> a, pair<double, double> b)
{
    double x1 = x - a.first, x2 = x - b.first;
    double y1 = y - a.second, y2 = y - b.second;    
    return ((x1*y2-y1*x2) < 0);
}

void setlast(double &x1, double &y1, double &x2, double &y2)
{    
    x2 = st[st.size()-1].first;
    y2 = st[st.size()-1].second;
    x1 = st[st.size()-2].first;
    y1 = st[st.size()-2].second;
}

符号改进我使用双打

    double sign(double x1,double y1, double x2,double y2, double y3,double x3)
    {
        double xx1 = x2 - x1, xx2 = x3 - x1;
        double yy1 = y2 - y1, yy2 = y3 - y1;
        return (xx1*yy2-yy1*xx2);
    }

int main()
{    
    int n;
    x = 0x3f3f3f3f;
    y = 0x3f3f3f3f;
    scanf("%d", &n);
    for(int i = 0; i < n; i++)
    {
        scanf("%lf %lf", &p[i].first, &p[i].second);
        if(p[i].first <= x && p[i].second <= y)
            x = p[i].first,
            y = p[i].second;
    }
    sort(p, p + n, f);
    p[n].first = x;
    p[n].second = y;
    st.pb(mp(p[0].first, p[0].second));
    st.pb(mp(p[1].first, p[1].second));
    double x1, x2, x3, y1, y2, y3;

这里我遍历所有向量并尝试确定凸包的点

    for(int i = 2; i < n; i++)
    {
        x3 = p[i].first;
        y3 = p[i].second;
        setlast(x1,y1,x2,y2);
        while(1)
            if(sign(x1,y1,x2,y2,x3,y3) < 0)
            {
                st.pb(mp(x3, y3));
                break;
            }
            else
                st.pop_back(),
                setlast(x1, y1, x2, y2);
    }

这里打印凸包

for(int i = 0; i < st.size(); i++)
        printf("%lf %lf\n", st[i].first, st[i].second);
    return 0
}

Answer 1

我的问题是，为什么int f(pair<int, int>, pair<int, int>)采用pair<int, int>代替pair<double, double>？

另外，为什么不将它命名为compare_blah之类的信息？

最后，为什么不返回bool而不是int？当然要么是有效的，但如果你返回一个bool，那么这只是作为一个比较函数更清楚。并且让阅读它的人明白你的程序应该是你的主要目标。让它做它应该做的是次要目标。毕竟，它所做的事情只是一种短暂的事态。最终有人会希望它做其他事情。

pair<int, int>可能就是你的问题。您在int和double之间的该函数中进行了几次隐式类型转换，并且左右丢失了信息。我怀疑那是你的意图。

如果您为pair typedef pair<double, double> point2d_t使用typedef，然后在任何地方使用point2d_t，您可以保护自己免受此类错误的影响，并使您的计划更便宜。

我对格雷厄姆的算法不太熟悉，无法评估你在abs内对f的使用情况，尽管对此发表评论的人很可能是正确的。

影响格雷厄姆寻找凸壳的算法的未知错误

1 个答案: