二维平面内判断点是否在一个简单多边形内部，在程序设计中我们一般采用射线法，或者内角和法。
如果这个简单多边形是一个凸多边形，可以在logN的时间复杂度内判断点是否在N个顶点的凸多边形中。
如图 判断点P是否在凸多边形内 设凸多边形顶点保存在convex[0..n-1]中
首先必须满足 向量convex[0]-P X convex[0]-convex[1]<0 convex[0]-P X convex[0]-convex[6]>0
X代表叉乘 如果允许点在多边形边上 <0 >0可以改写为<=0 >=0
这样确定点在角106范围内后 二分枚举2--6每个点x 利用叉乘可以判断向量convex[0]-P在convex[0]-convex[x]逆时针还是顺时针 如果是在逆时针 继续向x--6寻找顶点 否则在2--x-1寻找顶点 直到寻找到最接近P且在P顺时针方向的边convex[0]--convex[x]
此时可以确定点在角x-1 o x范围内
然后向量convex[x-1]-P X convex[x-1]-convex[x] 判断P是否在线段x-1--x的左侧 在左侧点在多边形内 否则不在
因为任意凸多边形都可以按照其中一个顶点三角剖分 所以有了logN的二分高效算法
/**************************************
Author : lxglbk
Date : 2012/08/17
Function : To determine whether a point is inside of a convex polygon in 2D plane
O(logN)
*********/
double cross(cpoint p0, cpoint p1, cpoint p2)
{
return (p1.x - p0.x) * (p2.y - p0.y) - (p2.x - p0.x) * (p1.y - p0.y);
}
bool Compl_inside_convex(const cpoint & p,cpoint *con,int n)
{
if(n<3) return false;
if(cross(con[0],p,con[1])>-eps) return false;
if(cross(con[0],p,con[n-1])>1;
if(cross(con[0],p,con[mid])>-eps)
{
line=mid;
j=mid-1;
}
else i=mid+1;
}
return cross(con[line-1],p,con[line])<-eps;
}