凸包的定义就是,由平面上n个点组成的集合,,其凸包就是包含这些点的最小凸边多边形。
给定二维平面上的点集,凸包就是将最外层的点连接起来构成的凸边多边形,它能包含点集中所有的点。
举例,如果是两个点,则包含这两点的凸包就是他们构成的线段;如果是三点,且不共线,则凸包就是以这三个点构成的三角形,共线,则是三点构成的直线,如果三点以上,如下图,
蛮力法求解
凸包多变形有一个性质,就是如果两点的连线属于凸包的边,则其余的点全部位于同一侧,我们可以利用这个性质,蛮力法通过扫面所有的点,如果满足以上性质,则是该凸包多边形的边。
我们利用直线方程,定义两个点的坐标是(x1,y1),(x2,y2):
则:ax+by=c(其中:a=y2-y1,b=x2-x1,c=x1y2-y1x2)
这样的直线,把平面分成了两个部分,其中一个平面满足ax+by>c,另外一个满足ax+by<c。
步骤:
1、以x,y轴,用两个循环,一次来扫面每个点。
2、求出a,b,c。
3、一个循环扫描剩余的点,判断出在直线的哪一边,在这个循环完毕判断是否全部位于一个方向是的话则这条边是凸边性的边。
完整代码(c++):
#include<iostream>
#include<cmath>
#include<vector>
#include<map>
using namespace std;
vector<pair<float, float> > convexHull(vector<pair<float, float> > & point){
vector<pair<float, float> > pointPairVec;
if(point.size() <= 2) //假如点少于2
return point;
//两个循环扫描每个点
for(int i=0; i<point.size()-1; i++){
for(int j=i+1; j<point.size(); j++){
float a = point[j].second - point[i].second; //a=y2-y1,b=x2-x1,c=x1y2-y1x2
float b = point[i].first - point[j].first;
float c = point[i].first*point[j].second - point[i].second*point[j].first;
vector<pair<float, float> > signEqualVec;
signEqualVec.push_back(point[i]);
signEqualVec.push_back(point[j]);
int sign1 = 0, sign2 = 0; //记录直线两边的点
//扫描直线外的所有的点
for(int k=0; k<point.size(); k++){
if(k==i || k==j)
continue;
else{
float sign = a*point[k].first + b*point[k].second - c;
if(sign == 0)
signEqualVec.push_back(point[k]);
else if(sign > 0)
sign1++;
else if(sign < 0)
sign2++;
}
}
if(sign1 > 0 && sign2 > 0)
continue;
else{
for(int l=0;l<signEqualVec.size(); l++)
pointPairVec.push_back(signEqualVec[l]);
}
}
}
return pointPairVec;
}
int main(){
vector<pair<float, float> > pointsArray(8);
pointsArray[0] = pair<float, float>(9.0, -16.0);
pointsArray[1] = pair<float, float>(9.5, 13.0);
pointsArray[2] = pair<float, float>(-36.0, -10.0);
pointsArray[3] = pair<float, float>(-8.5, -5.5);
pointsArray[4] = pair<float, float>(27.0, 7.0);
pointsArray[5] = pair<float, float>(20.5, 19.0);
pointsArray[6] = pair<float, float>(-22.5, 19.0);
pointsArray[7] = pair<float, float>(-14.0, 8.0);
vector<pair<float, float> > convex = convexHull(pointsArray);
for(int i=0; i<convex.size(); i++){
cout<<'('<<convex[i].first<<", "<<convex[i].second<<')'<<'\t';
if(i%2 == 1)
cout<<endl;
}
return 0;
}