首先解决问题:什么是半平面? 顾名思义,半平面就是指平面的一半,我们知道,一条直线可以将平面分为两个部分,那么这两个部分就叫做两个半平面。
然后,半平面怎么表示呢? 二维坐标系下,直线可以表示为ax + by + c = 0,那么两个半平面则可以表示为ax + by + c >= 0 和ax + by + c < 0,这就是半平面的表示方法。
还有,半平面的交是神马玩意? 其实就是一个方程组,让你画出满足若干个式子的坐标系上的区域(类似于线性规划的可行域),方程组就是由类似于上面的这些不等式组成的。
另外,半平面交可以干什么? 半平面交虽然说是半平面的问题,但它其实就是关于直线的问题。一个一个的半平面其实就是一个一个有方向的直线而已。
半平面交的一个重要应用就是求多边形的核 。 多边形的核又是神马玩意? 它是平面简单多边形的核是该多边形内部的一个点集,该点集中任意一点与多边形边界上一点的连线都处于这个多边形内部。就是一个在一个房子里面放一个摄像 头,能将所有的地方监视到的放摄像头的地点的集合即为多边形的核。经常会遇到让你判定一个多边形是否有核的问题。(以上是转载的http://blog.csdn.net/xuechelingxiao/article/details/40859973。 侵删)
关于nlogn半平面交算法的论文。。
https://wenku.baidu.com/view/c750720bf78a6529647d53ae.html
这个题其实就相当于模板了。。不上题面了。。
#include <iostream>
#include <stdio.h>
#include <string.h>
#include <queue>
#include <cmath>
#include <algorithm>
using namespace std;
const int maxn = 1550;
const double esp = 1e-8;
int n, pn, dq[maxn], top, bot;
struct Point {
double x, y;
}p[maxn];
struct Line {
Point a, b;
double angle;
}l[maxn];
//精度
int dcmp(double k) {
if (fabs(k) < esp)
return 0;
return k > 0 ? 1 : -1;
}
//叉积
double multi(Point p1, Point p2, Point p3) {
return (p2.x - p1.x) * (p3.y - p1.y) - (p3.x - p1.x) * (p2.y - p1.y);
}
bool cmp(const Line& l1, const Line& l2) {
int d = dcmp(l1.angle - l2.angle);
if (!d)
return dcmp(multi(l1.a, l2.a, l2.b)) < 0;
return d < 0;
}
void add_line(Line& l, double x1, double y1, double x2, double y2) {
l.a.x = x1;
l.a.y = y1;
l.b.x = x2;
l.b.y = y2;
l.angle = atan2(y2 - y1, x2 - x1);
}
//get 交点
void get_intersect(const Line& l1, const Line& l2, Point& p) {
double a1 = l1.b.y - l1.a.y;
double b1 = l1.a.x - l1.b.x;
double c1 = (l1.b.x - l1.a.x) * l1.a.y - (l1.b.y - l1.a.y) * l1.a.x;
double a2 = l2.b.y - l2.a.y;
double b2 = l2.a.x - l2.b.x;
double c2 = (l2.b.x - l2.a.x) * l2.a.y - (l2.b.y - l2.a.y) * l2.a.x;
p.x = (c2 * b1 - c1 * b2) / (a1 * b2 - a2 * b1);
p.y = (c1 * a2 - c2 * a1) / (a1 * b2 - a2 * b1);
}
bool judge(const Line& l1, const Line& l2, const Line& l3) {
Point p;
get_intersect(l2, l3, p);
return dcmp(multi(p, l1.a, l1.b)) > 0;
// > 0 p在向量l1.a->b 左边
}
void half_plan_intersect() {
sort(l, l + n, cmp);
int i, j;
for (i = 0, j = 0; i < n; ++i) {
if (dcmp(l[i].angle - l[j].angle) > 0)
l[++j] = l[i];
//去掉极角相同的
}
n = j + 1;
dq[0] = 0;
dq[1] = 1;
top = 1;
bot = 0;
for (i = 2; i < n; ++i) {
while (top > bot && judge(l[i], l[dq[top]], l[dq[top - 1]]))
top--;
while (top > bot && judge(l[i], l[dq[bot]], l[dq[bot + 1]]))
bot++;
dq[++top] = i;
}
while (top > bot && judge(l[bot], l[dq[top]], l[dq[top - 1]]))
top--;
while (top > bot && judge(l[top], l[dq[bot]], l[dq[bot + 1]]))
bot++;
dq[++top] = dq[bot];
for (pn = 0, i = bot; i < top; i++, pn++)
get_intersect(l[dq[i + 1]], l[dq[i]], p[pn]);
}
double get_area() {
if (pn < 3)
return 0;
double area = 0;
for (int i = 1; i < pn - 1; ++i) {
area += multi(p[0], p[i], p[i + 1]);
}
return fabs(area) / 2;
}
int main() {
//freopen("in.txt", "r", stdin);
int t, i;
cin >> t;
while (t--) {
cin >> n;
for (i = 0; i < n; ++i) {
cin >> p[i].x >> p[i].y;
}
for (i = 0; i < n - 1; ++i) {
add_line(l[i], p[i].x, p[i].y, p[i + 1].x, p[i + 1].y);
}
add_line(l[i], p[i].x, p[i].y, p[0].x, p[0].y);
half_plan_intersect();
printf("%.2lf\n", get_area());
}
}