2

我已经编写了格雷厄姆算法,但它仍然给了我凸包的错误点。我需要帮助。认为我的符号功能中有一个错误,但不知道它是什么。

#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
}
4

1 回答 1

0

我的问题,为什么int f(pair<int, int>, pair<int, int>)takepair<int, int>而不是pair<double, double>

另外,为什么不将它命名为信息丰富的东西compare_blah

最后,为什么它不返回bool而不是返回int?两者都可以,但如果你返回一个bool. 让阅读它的人清楚你的程序应该是你的主要目标。让它做它应该做的事情是次要目标。毕竟,它做它应该做的只是一个暂时的状态。最终有人会希望它做其他事情。

事情可能是你的pair<int, int>问题。您在该函数中进行了几次隐式类型转换,intdouble在左右丢失信息。我怀疑那是你的意图。

如果您愿意根据自己的pair喜好使用 typedef typedef pair<double, double> point2d_t,然后point2d_t在任何地方使用,您可以保护自己免受此类错误的影响,并使您的程序在交易中更清晰。

我对格雷厄姆的算法不够熟悉,无法评估您对absinside of的使用f,尽管对此发表评论的人很可能是正确的。

于 2012-11-04T20:47:39.750 回答