如何判断点是否在凸包或任意多边形内部

判断点否在凸包内部

例题

https://vjudge.net/problem/UVALive-7281
直接上例题，稍稍梳理一下，如下图

• 思考一下，最简单的办法，按照$Andrew$求凸包的思路，从凸包的某个点沿着逆时针顺序行进，同时看这个点是不是在这条线段左侧，如果到某条线段发现点在这条线段右侧，那么说明这个点就在凸包外侧，判断点在线段的哪一侧使用叉积即可
• 这样的判断时间复杂度是$O(n)$的，对于大规模询问显然不行
• 我们有更好的办法，我们在凸包上进行二分，看下面的图
  
• 这个点的位置只有四种情况（边界一样），在外侧的两种情况可以通过叉积来判断，这是非法情况，内侧的情况可能是合法情况，只需要判断一下这个点是不是在上面的线内部即可，使用叉积判断，比较简单就不细说了，做一下上面的那道题，整理一下模板

#include <bits/stdc++.h>

using namespace std;
#define db double
const db eps = 1e-10;
const int MAXN = 1e4 + 100;
int sgn(db x){
    if(fabs(x) < eps) return 0;
    return x < 0 ? -1 : 1;
}
struct Point{
    db x, y;
    Point(){}
    Point(db x, db y): x(x), y(y){}
    Point operator + (const Point &B)const{
        return Point(x + B.x, y + B.y);
    }
    Point operator - (const Point &B)const{
        return Point(x - B.x, y - B.y);
    }
    bool operator < (const Point &B)const{
        return sgn(x - B.x) < 0 || (sgn(x - B.x) == 0 && sgn(y - B.y) < 0);
    }
    bool operator == (const Point &B)const{
        return sgn(x - B.x) == 0 && sgn(y - B.y) == 0;
    }
}s[MAXN], ch[MAXN];
typedef Point Vector;
db Cross(Vector A, Vector B){
    return A.x * B.y - A.y * B.x;
}
int Convex_hull(int n){
    int v = 0;
    for(int i=0;i<n;i++){
        while(v > 1 && sgn(Cross(ch[v - 1] - ch[v - 2], s[i] - ch[v - 2])) <= 0){
            v -= 1;
        }
        ch[v++] = s[i];
    }
    int j = v;
    for(int i=n-2;i>=0;i--){
        while(v > j && sgn(Cross(ch[v - 1] - ch[v - 2], s[i] - ch[v - 2])) <= 0){
            v -= 1;
        }
        ch[v++] = s[i];
    }
    if(n > 1) v -= 1;
    return v;
}
bool check(Point A, int n){
    int l = 0;
    int r = n - 1;
    while(r - l > 1){
        int mid = ((r - l) >> 1) + l;
        db a1 = Cross(ch[mid] - ch[0], A - ch[0]);
        db a2 = Cross(ch[mid + 1] - ch[0], A - ch[0]);
        if(sgn(a1) >= 0 && sgn(a2) <= 0){
            if(sgn(Cross(ch[mid + 1] - ch[mid], A - ch[mid])) >= 0) return true;
            return false;
        }else if(sgn(a1) < 0){
            r = mid;
        }else{
            l = mid;
        }
    }
    return false;
}
int main(){
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    int n, m;
    db x, y;
    while(cin >> n){
        for(int i=0;i<n;i++){
            cin >> s[i].x >> s[i].y;
        }sort(s, s + n);
        n = unique(s, s + n) - s;
        int len = Convex_hull(n);
        cin >> m;
        int ans = 0;
        for(int i=0;i<m;i++){
            cin >> x >> y;
            if(check(Point(x, y), len)) ans += 1;
        }
        cout << ans << '\n';
    }
    return 0;
}

判断点是否在凹多边形内部

• 我们可以用一个简单的方法来判断，如果一个点在多边形内部，它向外的射线一定与多边形边缘有奇数个交点；否则有偶数个，这很容易理解，因为如果点在外部，进入多边形就一定要出去，如果在内部只要越过一次就可以出去了
• 首先应该确保给定的点是按顺序来的，然后我们可以找一个无穷远的点，枚举所有多边形上的边，然后看交点总数是奇数还是偶数即可

练习题

#include <bits/stdc++.h>

using namespace std;
#define db double
const db eps = 1e-10;
const int MAXN = 1e4 + 100;
int sgn(db x){
  if(fabs(x) < eps) return 0;
  return x < 0 ? -1 : 1;
}
struct Point{
  db x, y;
  Point(){}
  Point(db x, db y): x(x), y(y){}
  Point operator + (const Point &B)const{
      return Point(x + B.x, y + B.y);
  }
  Point operator - (const Point &B)const{
      return Point(x - B.x, y - B.y);
  }
  bool operator < (const Point &B)const{
      return sgn(x - B.x) < 0 || (sgn(x - B.x) == 0 && sgn(y - B.y) < 0);
  }
  bool operator == (const Point &B)const{
      return sgn(x - B.x) == 0 && sgn(y - B.y) == 0;
  }
};
typedef Point Vector;
db Cross(Vector A, Vector B){
    return A.x * B.y - A.y * B.x;
}
const int offset = 1e9;
/**
 * @brief 计算一个点是否在多边形内部
 * 
 * @param x 目标的横坐标
 * @param y 目标节点的纵坐标
 * @param coords 多边形点集
 * @return int 返回2表示在多边形边缘, 如果返回1表示在内部, 返回0表示在外部
 */
int isPointInPolygon(double x, double y, vector<double>& coords) {
  Point p;
  Point target = {x, y};
  int len = coords.size();
  vector<Point> vs(len >> 1);
  for(int i=1;i<len;i+=2){
    vs[i >> 1] = {coords[i - 1], coords[i]};
  }
  len >>= 1;
  int number = 0;// 射线与多边形边的交点个数
  int i = 0;
  while(i < len){
    p.x = (double)rand() + offset;// 随机一个点, 越远越好
    p.y = (double)rand() + offset;
    for(;i<len;i++){
      if(sgn(Cross(vs[i] - target, vs[(i + 1) % len] - target)) == 0 &&
          sgn((vs[i].x - target.x) * (vs[(i + 1) % len].x - target.x)) <= 0 &&
          sgn((vs[i].y - target.y) * (vs[(i + 1) % len].y - target.y)) <= 0){
        return 2;// 在多边形边缘 
      }else if(sgn(Cross(p - target, vs[i] - target)) == 0){// 在射线上, 换个p点
        break;
      }else if(sgn(Cross(vs[i] - target, p - target) * Cross(vs[(i + 1) % len] - target, p - target)) < 0
            && sgn(Cross(target - vs[i], vs[(i + 1) % len] - vs[i]) * Cross(p - vs[i], vs[(i + 1) % len] - vs[i])) < 0){
        // 判断两线段是否相交
        number += 1;
      }
    }
  }
  return number & 1;// 如果交点为奇数, 那么说明在内部
}