【模板】凸包问题

模板题：洛谷P2742 圈奶牛 [USACO5.1]Fencing the Cows

Part 1 Concept

凸包概念：
（非数学用语）简单来讲就是，二维坐标系中给定若干点坐标，尽可能少地选择其中若干个点进行两两连线，使连线能包含所有给定点。【精准数学定义详见计算几何中对凸包的定义】
凸包性质：

1. 如果要将给定所有点全部围住，且要求最终多边形周长最短，则凸包为这一问题的最优方案（证明：两点间直线最短）（注：凸包并非面积最小的最优解；面积最小 很可能是凹多边形）可借助下图思考。（图源：https://blog.csdn.net/viafcccy/article/details/87483567）
   
2. 点集中最靠下（y值最小）的点一定是凸包上的点；同理，最靠右、最靠左、最靠上的点，均一定为凸包上的点。
3. 根据性质2，延申性质3： 考虑到，如果我们将点集与坐标系同步旋转任意角度，性质2始终成立； 由此推出，点集任意方向向外延申，均有“该方向上最外部点，一定为凸包上的点”。反之，若某一点在任意方向上，都不是最外部的点，则该点一定不是凸包上的点。
4. 根据性质2,3，推得：若将某一凸包上点，连向其余各个点，夹角最大的两条线所连接的另外两点，一定也是凸包上的点，且三点为凸包上相邻点。(如：p0连出p1和p12；或者p3向外连，可以找到p1,p10）
5. 根据性质2,3，推证下述解法的成立性。

其实根据性质4，我们就可以暴力解题了，对于每个已知凸包点，排序其和其余点的夹角，依次得出每个点。复杂度O(n²logn)。暴力复杂度过高，因此优化。

Part 2 凸包算法

学名：Graham扫描。
时间复杂度O(nlogn)
总得来说就是利用了一些单调栈的思路。
Step 1: 根据性质2，我们不妨找最低的点作为起点 P0。
Step 2: 计算各点与p0的连线与y=yp0的夹角。并按此夹角从小到达排序。如果夹角相等，则按与p0距离从小到大（可以预见的是，排序后，点集顺序从P0开始逆时针旋转）
Step 3: 设立栈。任意时刻，栈中元素表示 当前我们认为是凸包上的点编号。初始状态是，将p0和p1放入栈中。【依据：性质4】
Step 4: 按先前排序，考虑下一点p2。暂时放入栈内。接着再考虑下一点p3，这时p3与前面两点所组成的连线有一定位置关系。
如果，点p3在直线p1p2的外部（即向量p1p3在向量p1p2的右侧），则p2一定不是凸包点，将其退栈; p3可能是凸包点，将其入栈。【证明：这一步其实相当于在判断，在向量p1p3的垂直方向上，p2是否为这一方向上的最外部点（性质3）】
如果，点p3在直线p1p2的内部或在直线上（即向量p1p3在向量p1p2的左侧），则p2目前为之仍是凸包上点；p3也可能是凸包点，将其入栈。
Step 5: 按顺序对各个入栈点进行step4，更新栈内元素。
最终栈内所留点即为凸包各点。

Part 3 Code

注1：step4中，栈顶top，当前在循环到第i个点，则while( Line( top-2,top-1) 在 Line( top-2, i) ) top–;。
注2：叉积的正负可用于判断向量的左右位置情况。
a x b=a.xb.y-a.yb.y; 正，a在右侧，负，b在右侧
叉积判断两线之间的位置关系

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const double INF=(1<<30);
const int N=1e5+5;
struct ss{
	double x,y;
	void operator =(ss a){
		x=a.x;y=a.y;
	}
	ss operator -(ss a){
		return (ss){x-a.x,y-a.y};
	}
}p[N];//点集
inline double rate(ss a,ss b){//直线ab斜率 
	if(a.x==b.x) return INF;
	return (a.y-b.y) / (a.x-b.x);
}
inline double Xji(ss a,ss b){//向量ca与cb的叉积 
	return a.x*b.y-b.x*a.y;
}
inline double dis(ss a,ss b){
	return sqrt( (a.y-b.y)*(a.y-b.y)+(a.x-b.x)*(a.x-b.x) );
}
bool cmp1(ss a,ss b){if(a.y==b.y)return a.x<b.x;return a.y<b.y;}//找最低点中最左的点 
bool cmp2(ss a,ss b){
//叉积=0，共线，取离得近的 
	double ls=Xji(a-p[0],b-p[0]);
	if(!ls) return dis(a,p[0])<dis(b,p[0]);
	return ls<0;
}
int n;
int stack[N],top;
double ans;
int main(){
	cin>>n;
	for(int i=1;i<=n;i++)
		cin>>p[i].x>>p[i].y;
	sort(p+1,p+n+1,cmp1);p[0]=p[1];
	sort(p+1,p+n+1,cmp2);int pos=3;
	stack[++top]=1;stack[++top]=2;
	while(pos<=n){
		int a,b,c;double rate1,rate2;
		a=stack[top-1],b=stack[top],c=pos;
		while(Xji(p[c]-p[a],p[b]-p[a])<=0&&top>=2){//别忘了卡空栈 
			top--;
			a=stack[top-1],b=stack[top],c=pos;
		}
		stack[++top]=pos;pos++;
	}
	stack[0]=stack[top];
	for(int i=1;i<=top;i++){
		int a=stack[i-1],b=stack[i];
		ans+=dis(p[a],p[b]);
	}
	printf("%.2f\n",ans);
	return 0;	
}

因为忘了防止top<0，RE了半天我也是醉了…