POJ - 3608 Bridge Across Islands（凸包+旋转卡壳模板 求两凸包点之间的最近距离）

链接：https://cn.vjudge.net/problem/POJ-3608

题意：求两凸包点之间的最近距离，注意这个意思是说，如果求点到线段的距离，那么要保证垂足落到线段上。

思路：枚举一个凸包的边，去找另一凸包中的最远点，和求在一个凸包中的对踵点差不多。

PS:本想自己整理个模板，发现还是kuangbin大大的模板好用。

#include 
#include 
#include 
#include 
#define ll long long
using namespace std;
const int N = 1e4+10;
const double eps = 1e-8;
int sgn(double x)
{
	if(fabs(x)0) return 1;
	else return 0;
}
//线段L上，距离P最近的点 
Point NearestPointToLineSeg(Point P,Line L)
{
    Point result;
    double t = ((P-L.s)*(L.e-L.s))/((L.e-L.s)*(L.e-L.s));
    if(t >= 0 && t <= 1)
    {
        result.x = L.s.x + (L.e.x - L.s.x)*t;
        result.y = L.s.y + (L.e.y - L.s.y)*t;
    }
    else
    {
        if(dis(P,L.s) < dis(P,L.e))
            result = L.s;
        else result = L.e;
    }
    return result;
}
//点p0到线段p1p2的距离
double pointtoseg(Point p0,Point p1,Point p2)
{
    return dis(p0,NearestPointToLineSeg(p0,Line(p1,p2)));
}
//平行线段p0p1和p2p3的距离
double dispallseg(Point p0,Point p1,Point p2,Point p3)
{
    double ans1 = min(pointtoseg(p0,p2,p3),pointtoseg(p1,p2,p3));
    double ans2 = min(pointtoseg(p2,p0,p1),pointtoseg(p3,p0,p1));
    return min(ans1,ans2);
}
//得到向量a1a2和b1b2的位置关系
double Get_angle(Point a1,Point a2,Point b1,Point b2)
{
    Point t = b1 - ( b2 - a1 );
    return (a2-a1)^(t-a1);
} 
//求凸包时比较 
bool check(Point a,Point b,Point c)
{
	x=(b-a)^(c-a);
	return sgn(x)<=0;
}
//将点按逆时针排序 
void sort_point(Point *p,int n)
{
		pos=0;		
		for(int i=0;i1&&check(st[top-1],st[top],p[i])) 
			top--;			
		st[++top]=p[i];
	}	
	n=top+1;
	for(int i=0;i0)
			qq=i;
	for(int i=0;i