半平面交总结and模板
博客原文地址:http://blog.csdn.net/xuechelingxiao/article/details/40859973
这两天刷了POJ上几道半平面交,对半平面交有了初步的体会,感觉半平面交还是个挺有用的知识点。
半平面交主要是看的ZZY的国家队论文,他提出的是一种O(n×log(n))的排序增量法。
附论文地址: 算法合集之《半平面交的新算法及其有用价值》。
题目大意:
World finals 要開始了,比赛场地是一个多边形。裁判跟教练坐在多边形的边上,问场地里是不是存在一个区域使得这个区域内的每一个点都能看到多边形边上的计分板。
解题思路:
非常单纯的推断多边形的核是不是存在。
题目大意:
一个360度摄像头。给你一个多边形。问多边形中是否存在一个区域。将摄像头放进去能够监控整个多边形区域。
解题思路:
相同的半平面交推断多边形是否存在核。
POJ 3130 How I Mathematician Wonder What You Are!
题目大意:
给出星型的定义:一个多边形F是星形的条件是有一点C∈F,随意点P∈F。线段CP都在多边形F内。
给出一个多边形。推断他是不是一个星型。
解题思路:
一样的求多边形的核是否存在。
题目大意:
给出一个画廊(不一定是凸多边形),找到一块区域。使得这块区域内的每一个点都能看到整个画廊。求出这块区域的面积。
解题思路:
半平面交求多边形的核,须要求核的面积。
id=2451">POJ 2451 Uyuw's Concert
题目大意:
在一个10000*10000的正方形内举办演唱会,中间会有一些警戒线,舞台仅仅能布置在全部警戒线的左側,给定一些线段,问最后舞台的面积是多少。
解题思路:
这是zzy当初为他半平面交排序增量法而出的一个题。可是后来POJ改过数据,所以O(n*n)的算法也是能够过去的。
事实上还是一个半平面交,求交出来的面积的问题。
半平面交的题就做了差点儿相同这些,主要还是看出来题目类型是求半平面交就能够了,剩下的就好说了。
因为题目都差点儿相同,代码没有太大的变化,所以就仅仅贴一个代码了。
POJ 2451 的AC代码,能够当做模板来用了,当时还由于数组开小了而WA了10多遍。。。。。。为什么不给我返回个RE。。!
#include <stdio.h>
#include <math.h>
#include <algorithm>
using namespace std;
const double eps = 1e-8;
const int maxn = 20010; int tail, head, it;
struct Point {
double x, y;
} P[maxn];
struct Line {
Point a, b;
double angle;
void getAngle() {angle = atan2(b.y-a.y, b.x-a.x);}
} L[maxn], deq[maxn]; int dcmp(double x) {
return x < -eps ? -1 : x > eps;
}
double xmult(Point a, Point b, Point c) {
return (a.x-c.x)*(b.y-c.y)-(a.y-c.y)*(b.x-c.x);
}
bool cmp(Line u, Line v) {
int d = dcmp(u.angle-v.angle);
if(d) return d > 0;
return dcmp(xmult(u.a, v.a, v.b)) > 0;
///Clockwise:大于0取向量左半部分为半平面,小于0。取右半部分
}
Point intersection(Line u, Line v)
{
Point ret = u.a;
double t = ((u.a.x-v.a.x)*(v.a.y-v.b.y)
-(u.a.y-v.a.y)*(v.a.x-v.b.x))
/ ((u.a.x-u.b.x)*(v.a.y-v.b.y)
-(u.a.y-u.b.y)*(v.a.x-v.b.x));
ret.x += (u.b.x-u.a.x)*t, ret.y += (u.b.y-u.a.y)*t;
return ret;
}
bool judge(Line l1, Line l2, Line l3) {
Point p = intersection(l2, l3);
return dcmp(xmult(p, l1.a, l1.b)) < 0;
///Clockwise:大于小于符号与上面cmp()中凝视处相反
}
void HPI(Line L[], int n){
sort(L, L+n, cmp);
int tmp = 1;
for(int i = 1; i < n; ++i) {
if(dcmp(L[i].angle-L[tmp-1].angle) != 0) {
L[tmp++] = L[i];
}
}
n = tmp;
deq[0] = L[0], deq[1] = L[1];
head = 0, tail = 1;
for(int i = 2; i < n; ++i) {
while(head < tail && judge(L[i], deq[tail-1], deq[tail]))
tail--;
while(head < tail && judge(L[i], deq[head+1], deq[head]))
head++;
deq[++tail] = L[i];
}
while(head < tail && judge(deq[head], deq[tail-1], deq[tail]))
tail--;
while(head < tail && judge(deq[tail], deq[head+1], deq[head]))
head++;
if(head == tail) return ;
it = 0;
for(int i = head; i < tail; ++i) {
P[it++] = intersection(deq[i], deq[i+1]);
}
if(tail > head+1) {
P[it++] = intersection(deq[head], deq[tail]);
}
}
double getArea(Point p[], int n) {
double area = 0;
for(int i = 1; i < n-1; ++i) {
area += xmult(P[0], P[i], P[i+1]);
}
return fabs(area)/2.0;
} int main()
{
int n;
while(~scanf("%d", &n)) {
n += 4;
L[0] = (Line){(Point){0, 10000}, (Point){0, 0}};
L[1] = (Line){(Point){10000, 10000}, (Point){0, 10000}};
L[2] = (Line){(Point){10000, 0}, (Point){10000, 10000}};
L[3] = (Line){(Point){0, 0}, (Point){10000, 0}};
L[0].getAngle(), L[1].getAngle(), L[2].getAngle(), L[3].getAngle();
for(int i = 4; i < n; ++i) {
scanf("%lf%lf%lf%lf", &L[i].a.x, &L[i].a.y, &L[i].b.x, &L[i].b.y);
L[i].getAngle();
}
HPI(L, n);
printf("%.1f\n", getArea(P, it));
} return 0;
}
半平面交总结and模板的更多相关文章
- 半平面交模板(O(n*n)&& O(n*log(n))
摘自http://blog.csdn.net/accry/article/details/6070621 首先解决问题:什么是半平面? 顾名思义,半平面就是指平面的一半,我们知道,一条直线可以将平面分 ...
- POJ 半平面交 模板题 三枚
给出三个半平面交的裸题. 不会的上百度上谷(gu)歌(gou)一下. 毕竟学长的语文是体育老师教的.(卡格玩笑,别当真.) 这种东西明白就好,代码可以当模板. //poj1474 Video Surv ...
- 再来一道测半平面交模板题 Poj1279 Art Gallery
地址:http://poj.org/problem?id=1279 题目: Art Gallery Time Limit: 1000MS Memory Limit: 10000K Total Su ...
- bzoj 2618 半平面交模板+学习笔记
题目大意 给你n个凸多边形,求多边形的交的面积 分析 题意\(=\)给你一堆边,让你求半平面交的面积 做法 半平面交模板 1.定义半平面为向量的左侧 2.将所有向量的起点放到一个中心,以中心参照进行逆 ...
- [模板] 计算几何2: 自适应Simpson/凸包/半平面交/旋转卡壳/闵可夫斯基和
一些基本的定义在这里: [模板] 计算几何1(基础): 点/向量/线/圆/多边形/其他运算 自适应Simpson Simpson's Rule: \[ \int ^b_a f(x)dx\approx ...
- POJ 3525 /// 半平面交 模板
题目大意: 给定n,接下来n行逆时针给定小岛的n个顶点 输出岛内离海最远的点与海的距离 半平面交模板题 将整个小岛视为由许多半平面围成 那么以相同的比例缩小这些半平面 一直到缩小到一个点时 那个点就是 ...
- 三道半平面交测模板题 Poj1474 Poj 3335 Poj 3130
求半平面交的算法是zzy大神的排序增量法. ///Poj 1474 #include <cmath> #include <algorithm> #include <cst ...
- bzoj 2618【半平面交模板】
#include<iostream> #include<cstdio> #include<algorithm> #include<cmath> usin ...
- 【模板】凸包向内推进求不严格的半平面交——poj3384
想不明白这题写严格的半平面交为什么会错 /* 凸包所有边向内推进r */ #include<iostream> #include<cstring> #include<cs ...
随机推荐
- Google Map API V2密钥申请
之前用的都是v1,用的是MapView,好吧,仅仅能认命了.废话不再多说,開始android 的Google Maps Android API v2吧 之前參考了http://www.cnblogs. ...
- C# 6.0 (C# vNext) 的新功能:Expression Bodied Functions and Properties
Expression Bodied Function 它可以用在: methods user-defined operators type conversions read-only properti ...
- 采用PopupWin控制消息推送功能
最近的项目需要,急需实现消息推送功能.接连试了很多办法,让我们用JavaScript为了实现啊,其效果是不咋好,最后,我发现了一个PopupWin这个,看着眼前的成绩不错,开始使用. 1.准备工作.先 ...
- 【Android开发经验】使用反射,得到的类的字段、方法、并实现了简单的调用
本文后推出Android的ICO框架做准备,所以,假设你想要一个最近的一项研究Android的ICO学生框架.你可以稍微看一下. 首先,简介一下Java里面的反射. JAVA反射机制是在执行状态中,对 ...
- Andorid Async-HttpClient阅览
Async-httpclient它是一个用于Android应用程序开发http访问开源框架.开源GitHub在,因为今天GitHub下载没有看到.我在这里提供下载地址.是之前从GitHub上下载的,版 ...
- c++日历改进版
#include<iostream> # include<fstream> #include<time.h> #include<string> #inc ...
- POJ 1915 Knight Moves(BFS+STL)
Knight Moves Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 20913 Accepted: 9702 ...
- JavaFX横幅类游戏开发 教训 游戏贴图
上一节课,我们即将完成战旗Demo有了一个大概的了解.教训这,我们将学习绘制游戏地图. 由于JavaFX 2.2中添加了Canvas相关的功能,我们就能够使用Canvas来实现游戏绘制了. 游戏地图绘 ...
- 【转】linux建立软链接
实例:ln -s /home/gamestat /gamestat linux下的软链接类似于windows下的快捷方式 ln -s a b 中的 a 就是源文件,b是链接文件名,其作用是当进入 ...
- std::list.pop_back() 弹空了列表导致的崩溃
core文件输出: (gdb) bt # # ) at xxxxx/sql/signal_handler.cc: # <signal handler called> # # # # # 0 ...