凸包问题（含叉积讲解）

预备知识

叉积

1.概念的引入

在平面中我们为了度量一条直线的倾斜状态，为引入倾斜角这个概念。而通过在直角坐标系中建立tan α = k，我们实现了将几何关系和代数关系的衔接，这其实也是用计算机解决几何问题的一个核心，计算机做的是数值运算，因此你需要做的就是把几何关系用代数关系表达出来。而在空间中，为了表示一个平面相对空间直角坐标系的倾斜程度，我们利用一个垂直该平面的法向量来度量(因为这转化成了描述直线倾斜程度的问题)。

2.定义

3.应用

（1）求解三角形（平行四边形）面积

相对海伦公式更精准（减少了开根号的误差）
基于这个结论还可以推出n边形面积计算公式

严谨的来说，用叉乘算面积也是能够算出负值的

（2）点定位

（3）极角排序

极角，指的就是以x轴正半轴为始边，建立极坐标，逆时针转过的角，这个角的范围是[0,2π]。

实现方法：

1.用叉积计算极角（精度高，时间慢）

struct point
{
   
    double x,y;
    point(double x=0, double y=0):x(x), y(y){
   }
    point operator - (const point &t)const
    {
   
        return point(x-t.x, y-t.y);
    }//a - b
    double operator *(const point &t)const
    {
   
        return x*t.x + y*t.y;
    }//a * b
    double operator ^(const point &t)const
    {
   
        return x*t.y - y*t.x;
    }//a X b
};
double compare(point a,point b,point c)//计算极角 ab × ac 
{
   
    return (b-a)^(c-a);
}
bool cmp(point a,point b)
{
   
	double f=compare(p[pos],a,b);
	if(f==0) return a.x-p[pos].x<b.x-p[pos].x;
	else if(f>0) return true;
	else return false;
}

如果取的点不是边角的点，那么需要先按照象限排序。

int Quadrant(point a)//象限排序，注意包含四个坐标轴
{
   
    if(a.x>0&&a.y>=0)  return 1;
    if(a.x<=0&&a.y>0)  return 2;
    if(a.x<0&&a.y<=0)  return 3;
    if(a.x>=0&&a.y<0)  return 4;
}


bool cmp2(point a,point b)//先象限后极角
{
   
    if(Quadrant(a)==Quadrant(b))//返回值就是象限
        return cmp(a,b);
    else Quadrant(a)<Quadrant(b);
}

2.atan2函数（时间快，精度较差）

atan2(y,x)函数返回的是原点至点(x,y)的方位角，即与 x 轴的夹角。也可以理解为复数 x+yi 的辐角。返回值的单位为弧度，这里的这个极角的范围是（−π,π]， 一二象限为正，三四象限为负（结果为正表示从 X 轴逆时针旋转的角度，结果为负表示从 X 轴顺时针旋转的角度）。所以我们从小到大排完序后，实际上是 第三象限<第四象限<第一象限<第二象限

struct point
{
   
	double x,y;
	double angle;
	bool operator <(const point &t)
	{
   
		return angle<t.angle;
	}
}p[N];
bool cmp(point a,point b)
{
   
	if(a.angle==b.