在本教程中,我们将讨论一个程序来查找给定点集的凸包。
凸包是最小的多边形凸图,其中包含图内边界上的所有给定点。
在此程序中,我们将使用蛮力将给定的点划分为较小的段,然后最终将随后的段合并以构造凸包。
示例#include
using namespace std;
//存储多边形的中心点
pair mid;
//的象限
//特定点
int quad(pair p){
if (p.first >= 0 && p.second >= 0)
return 1;
if (p.first <= 0 && p.second >= 0)
return 2;
if (p.first <= 0 && p.second <= 0)
return 3;
return 4;
}
//如果线接触多边形
int calc_line(pair a, pair b,
pair c){
int res = (b.second-a.second)*(c.first-b.first) -
(c.second-b.second)*(b.first-a.first);
if (res == 0)
return 0;
if (res > 0)
return 1;
return -1;
}
//比较功能进行排序
bool compare(pair p1, pair q1){
pair p = make_pair(p1.first - mid.first,
p1.second - mid.second);
pair q = make_pair(q1.first - mid.first,
q1.second - mid.second);
int one = quad(p);
int two = quad(q);
if (one != two)
return (one
return (p.second*q.first
}
//查找两个多边形的上切线
vector> merger(vector > a,
vector > b){
int n1 = a.size(), n2 = b.size();
int ia = 0, ib = 0;
for (int i=1; i
if (a[i].first > a[ia].first)
ia = i;
//计算b的最左点
for (int i=1; i
if (b[i].first
ib=i;
int inda = ia, indb = ib;
bool done = 0;
while (!done){
done = 1;
while (calc_line(b[indb], a[inda], a[(inda+1)%n1]) >=0)
inda = (inda + 1) % n1;
while (calc_line(a[inda], b[indb], b[(n2+indb-1)%n2]) <=0){
indb = (n2+indb-1)%n2;
done = 0;
}
}
int uppera = inda, upperb = indb;
inda = ia, indb=ib;
done = 0;
int g = 0;
//计算下切线
while (!done){
done = 1;
while (calc_line(a[inda], b[indb], b[(indb+1)%n2])>=0)
indb=(indb+1)%n2;
while (calc_line(b[indb], a[inda], a[(n1+inda-1)%n1])<=0){
inda=(n1+inda-1)%n1;
done=0;
}
}
int lowera = inda, lowerb = indb;
vector> ret;
//合并两个多边形以获得凸包
int ind = uppera;
ret.push_back(a[uppera]);
while (ind != lowera){
ind = (ind+1)%n1;
ret.push_back(a[ind]);
}
ind = lowerb;
ret.push_back(b[lowerb]);
while (ind != upperb){
ind = (ind+1)%n2;
ret.push_back(b[ind]);
}
return ret;
}
//蛮力算法找到凸包
vector> bruteHull(vector> a){
set >s;
for (int i=0; i
for (int j=i+1; j
int x1 = a[i].first, x2 = a[j].first;
int y1 = a[i].second, y2 = a[j].second;
int a1 = y1-y2;
int b1 = x2-x1;
int c1 = x1*y2-y1*x2;
int pos = 0, neg = 0;
for (int k=0; k
if (a1*a[k].first+b1*a[k].second+c1 <= 0)
neg++;
if (a1*a[k].first+b1*a[k].second+c1 >= 0)
pos++;
}
if (pos == a.size() || neg == a.size()){
s.insert(a[i]);
s.insert(a[j]);
}
}
}
vector>ret;
for (auto e:s)
ret.push_back(e);
//沿逆时针方向移动
mid = {0, 0};
int n = ret.size();
for (int i=0; i
mid.first += ret[i].first;
mid.second += ret[i].second;
ret[i].first *= n;
ret[i].second *= n;
}
sort(ret.begin(), ret.end(), compare);
for (int i=0; i
ret[i] = make_pair(ret[i].first/n, ret[i].second/n);
return ret;
}
//返回凸包的值
vector> divide(vector> a){
if (a.size() <= 5)
return bruteHull(a);
//左包含左半点
//右边包含右边的半点
vector>left, right;
for (int i=0; i
left.push_back(a[i]);
for (int i=a.size()/2; i
right.push_back(a[i]);
vector>left_hull = divide(left);
vector>right_hull = divide(right);
//合并凸包
return merger(left_hull, right_hull);
}
int main(){
vector > a;
a.push_back(make_pair(0, 0));
a.push_back(make_pair(1, -4));
a.push_back(make_pair(-1, -5));
a.push_back(make_pair(-5, -3));
a.push_back(make_pair(-3, -1));
a.push_back(make_pair(-1, -3));
a.push_back(make_pair(-2, -2));
a.push_back(make_pair(-1, -1));
a.push_back(make_pair(-2, -1));
a.push_back(make_pair(-1, 1));
int n = a.size();
sort(a.begin(), a.end());
vector >ans = divide(a);
cout <
for (auto e:ans)
cout <
return 0;
}
输出结果Convex Hull:
-5 -3
-1 -5
1 -4
0 0
-1 1