平面几何常用算法.pptVIP

* A E D C B 只需要把多边形的每一条边看成一条直线。这n条直线所组成的半平面交就是答案。（如图中绿色边框围起来的多边形中任意一点都满足条件） 增量法求半平面交： 按照逆时针顺序考虑多边形的所有顶点，保留在直线左侧（如果是顺时针给出当然是右侧咯）和直线上的点，删除直线右边的点。如果有向直线和多边形相交时产生了新的点，这些点应加在半平面交上。具体来说考虑完一个顶点p[i]，在考虑p[i+1]之前要先判断p[i]p[i+1]是否和别的直线在多边形内部相交，如果是，则吧该点加入到半平面交上（如左图）只要最终加入的点的个数不为0就存在。。 * 半平面交模板 #include iostream #include cstdio #include algorithm #include cmath using namespace std; const double eps = 1e-8; const int maxn = 105; int dq[maxn], top, bot, pn, order[maxn], ln; //dp数组最终保存半平面交的交点 struct Point { double x, y; } p[maxn];//保存输入的点 struct Line { Point a, b;double angle; } l[maxn]; //没存每两个相邻点组成的直线 int dblcmp(double k) { if (fabs(k) eps) return 0; return k 0 ? 1 : -1; } double multi(Point p0, Point p1, Point p2) //向量的叉积运算 { return (p1.x-p0.x)*(p2.y-p0.y)-(p1.y-p0.y)*(p2.x-p0.x); } bool cmp(int u, int v) { int d = dblcmp(l[u].angle-l[v].angle); if (!d) return dblcmp(multi(l[u].a, l[v].a, l[v].b)) 0; return d 0; } void getIntersect(Line l1, Line l2, Point p) //得到L1 L2的交点p { double dot1,dot2; dot1 = multi(l2.a, l1.b, l1.a); dot2 = multi(l1.b, l2.b, l1.a); p.x = (l2.a.x * dot2 + l2.b.x * dot1) / (dot2 + dot1); p.y = (l2.a.y * dot2 + l2.b.y * dot1) / (dot2 + dot1); } bool judge(Line l0, Line l1, Line l2) { Point p; getIntersect(l1, l2, p); return dblcmp(multi(p, l0.a, l0.b)) 0; } * void addLine(double x1, double y1, double x2, double y2) //每相邻的两个点组成一条直线 { l[ln].a.x = x1;l[ln].a.y = y1;l[ln].b.x = x2;l[ln].b.y = y2;l[ln].angle = atan2(y2-y1, x2-x1); order[ln] = ln; //order是各条直线的标号 ln++; } void halfPlaneIntersection() //增量法求半平面交 { int i, j; sort(order, order+ln, cmp); for (i = 1, j = 0; i ln; i++) if (dblcmp(l[order[i]].angle-l[order[j]].angle) 0) order[++j] = order[i]; ln = j + 1;dq[0] = order[0];dq[1] = order[1];bot = 0;top = 1; for (i = 2; i ln; i++) { while (bot top judge(l[order[i]], l[dq[top-1]], l[dq[top]])) top--; while (bot top judge(