Question

我试图调试＆＃34;凸包＆＃34;贾维斯的算法。＆＃34;凸壳＆＃34;问题是，给定平面中n个点的集合P，找到形成包含所有其他点的凸多边形的顶点的子集CH（P）。递归写这个函数，但永远坚持循环并返回分段错误

    int main()
{
    vector<Point> hull(20);
    int n,x,y;
    cin >> n;
    vector<Point> ps;
    Point p;
//  Point p1,p2,q1,q2;

    while(cin >> x >> y)
    {   
        p.x = x;
        p.y = y;
        ps.push_back(p);
    }
    int base = find_leftmost_point(ps, n);
    hull.push_back(ps[base]);
    vector<Point> po = convexHull(ps, base, hull);
    cout << perimeter(po) << endl;
    return 0;
}

vector<Point> convexHull(vector<Point> points, int base, vector<Point> &hull)
{

    int p, q;

    p = base;   
    q = (p+1) % points.size();
    if (points.size() <= 3) 
    {

        return hull;
    }
    if(q == base)
    {
        return hull;
    }
    else
    {   
                for (int i = 0; i < points.size(); i++)
        {   
            if (orientation(points[p], points[i], points[q]) == 2)
            {
                q = i;
            }

        }
        cout<<points[q].x<<points[q].y<<endl;
        hull.push_back(points[q]);
        return convexHull(points, q, hull);
    }
}

double perimeter(vector<Point> P)
{
    double r = 0;
    for(int i = 1;i < P.size(); i++)
        r += sqrt(pow(P[i].x - P[i-1].x, 2) + pow(P[i].y - P[i-1].y, 2));
    return r;
}
int orientation(Point p, Point q, Point r)
{
    int val = (q.y - p.y) * (r.x - q.x) - (q.x - p.x) * (r.y - q.y);

    if (val == 0) 
        return 0;
    return (val > 0) ? 1 : 2;
}

int find_leftmost_point(vector<Point> points, int n)
{
    int l = 0;
    for (int i = 1; i < n; i++)
        if (points[i].x < points[l].x)
            l = i;
    return l;

}

Answer 1

你当然可以返回矢量。这本身并不会导致段错误。可能导致此类错误的原因是：

hull.push_back(points[p]);因为无法保证向量中有p+1个点
orientation(points[p], points[i], points[q])因为无法保证您的媒介中有p+1或n或q点
你最终会进行无限递归。应使用您使用的输入数据进行分析。但是对于if n＆gt; base，你最终会递归调用函数而没有可以阻止它减少的整数。

编辑＆amp;解决方案：

根据您提供的其他信息，确保n==points.size()和base<n。从那里可以清楚地看出，p，i和q总是小于n。这消除了两个首先可能的错误。

但是用一个小样本运行你的代码，表明你无休止地骑自行车：一旦你将最后一个点添加到船体，你又开始添加第一个。缺少什么，所以要确保你添加的点不在船体中。

您可以通过在for循环之后添加以下代码来执行此操作：

    auto srch=find(hull.begin(), hull.end(), points[q]); 
    if (srch!=hull.end()) {
        cout << "ALREADY IN"<<endl; 
        return hull;
    }

这里是online demo。

c ++凸包在递归方法中

1 个答案: