bzoj 2618【半平面交模板】
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=505;
int d,b,n,m;
struct dian
{
double x,y;
dian(double X=0,double Y=0)
{
x=X,y=Y;
}
dian operator + (const dian &a) const
{
return dian(x+a.x,y+a.y);
}
dian operator - (const dian &a) const
{
return dian(x-a.x,y-a.y);
}
dian operator * (const double &a) const
{
return dian(x*a,y*a);
}
dian operator / (const double &a) const
{
return dian(x/a,y/a);
}
}a[N],p[N];
struct bian
{
dian s,v;//s表示向量的起点,v表示向量的方向和长度(从(0,0)射出)
bian(dian S=dian(),dian V=dian())
{
s=S,v=V;
}
}l[N],q[N];
double cj(dian a,dian b)//叉积
{
return a.x*b.y-a.y*b.x;
}
double mj(dian a,dian b,dian c)//求有向面积
{
return cj(b-a,c-a)/2.0;
}
dian jd(bian x,bian y)//求交点
{
return x.s+x.v*(cj(x.s-y.s,y.v)/cj(y.v,x.v));
}
bool px(bian x,bian y)//判断平行
{
return cj(y.v,x.v)==0;
}
bool bn(bian x,bian y)//x在y的逆时针方向(平行先左后右
{
double ar=cj(x.v,y.v);
return (ar>0)||((ar==0)&&cj(x.v,y.s-x.s)>0);
}
bool dn(dian x,bian y)//点在线的逆时针方向
{
return cj(y.v,x-y.s)<=0;
}
bool cmp(const bian &x,const bian &y)//极角排序
{
if(x.v.y==0&&y.v.y==0)//同与x轴平行
return x.v.x<y.v.x;
if((x.v.y<=0)==(y.v.y<=0))//同在x轴上或下(包括x轴)
return bn(x,y);
return x.v.y<y.v.y;//一上一下下在前
}
int main()
{
scanf("%d",&d);
for(int i=1;i<=d;i++)
{
scanf("%d",&b);
for(int j=1;j<=b;j++)
{
int x,y;
scanf("%d%d",&x,&y);
p[n+j].x=x,p[n+j].y=y;//cout<<p[n+j].x<<" "<<p[n+j].y<<endl;
if(j!=1)
l[++m]=bian(p[n+j-1],p[n+j]-p[n+j-1]);
}
n+=b;
l[++m]=bian(p[n],p[n-b+1]-p[n]);
}
//半平面交
sort(l+1,l+m+1,cmp);
int top=0;
for(int i=1;i<=m;i++)
{
if(i==1||!px(l[i],l[i-1]))//去掉平行边
top++;
l[top]=l[i];
}
m=top;
int ll=1,rr=2;
q[1]=l[1],q[2]=l[2];
for(int i=3;i<=m;i++)//每次新加入向量,就会删掉在向量右边的交点(线上的也要删),维护的凸包首尾都是要删除的,最后还要模拟插入队头,把队尾中多余的半平面去掉
{
while(ll<rr&&dn(jd(q[rr],q[rr-1]),l[i]))
rr--;
while(ll<rr&&dn(jd(q[ll],q[ll+1]),l[i]))
ll++;
s[++rr]=l[i];
}
while(ll<rr&&dn(jd(q[rr],q[rr-1]),q[ll]))
rr--;//cout<<rr<<endl;
if(rr-ll<=1)
{
puts("0.000");
return 0;
}
top=0;
q[ll-1]=q[rr];
for(int i=ll;i<=rr;i++)
a[++top]=jd(q[i],q[i-1]);//求出相邻两边的交点,转化为凸包的记录方法
double ans=0.0;
for(int i=3;i<=top;i++)
ans+=mj(a[1],a[i-1],a[i]);
printf("%.3f\n",ans);
return 0;
}
/*
2
6
-2 0
-1 -2
1 -2
2 0
1 2
-1 2
4
0 -3
1 -1
2 2
-1 0
*/
bzoj 2618【半平面交模板】的更多相关文章
- bzoj 2618 半平面交模板+学习笔记
题目大意 给你n个凸多边形,求多边形的交的面积 分析 题意\(=\)给你一堆边,让你求半平面交的面积 做法 半平面交模板 1.定义半平面为向量的左侧 2.将所有向量的起点放到一个中心,以中心参照进行逆 ...
- POJ 3525 /// 半平面交 模板
题目大意: 给定n,接下来n行逆时针给定小岛的n个顶点 输出岛内离海最远的点与海的距离 半平面交模板题 将整个小岛视为由许多半平面围成 那么以相同的比例缩小这些半平面 一直到缩小到一个点时 那个点就是 ...
- 半平面交模板(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 ...
- 半平面交模板(BZOJ1007)
#include<cstdio> #include<algorithm> #define LDB long double using namespace std; ]; str ...
- 【BZOJ 2618】 2618: [Cqoi2006]凸多边形 (半平面交)
2618: [Cqoi2006]凸多边形 Description 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n=2时,两个凸多边形如下图: 则相交部分的面积为5.233. Input 第一 ...
- 2018.07.04 BZOJ 2618 Cqoi2006凸多边形(半平面交)
2618: [Cqoi2006]凸多边形 Time Limit: 5 Sec Memory Limit: 128 MB Description 逆时针给出n个凸多边形的顶点坐标,求它们交的面积.例如n ...
- BZOJ - 2618 凸多边形 (半平面交)
题意:求n个凸多边形的交面积. 半平面交模板题. #include<bits/stdc++.h> using namespace std; typedef long long ll; ty ...
随机推荐
- SOJ 4482 忽悠大神【最小生成树】
题目链接: http://acm.scu.edu.cn/soj/problem.action?id=4482 题意: 给定边权和点权,从一个点出发并回到该点,减少尽量多的边,每路过点和边都要把权重加到 ...
- P1396 营救 洛谷
https://www.luogu.org/problem/show?pid=1396 题目描述 “咚咚咚……”“查水表!”原来是查水表来了,现在哪里找这么热心上门的查表员啊!小明感动的热泪盈眶,开起 ...
- Java多线程分析案例
1. 多线程的创建方式 (1).继承 Thread类:但Thread本质上也是实现了Runnable 接口的一个实例,它代表一个线程的实例,并且,启动线程的唯一方法就是通过 Thread 类的 sta ...
- MongoDB学习day07--mongoose入门,数据库增删改查,默认参数,模块化
一.mongoose介绍 Mongoose 是在 node.js 异步环境下对 mongodb 进行便捷操作的对象模型工具. Mongoose 是 NodeJS 的驱动, 不能作为其他语言的驱动. M ...
- Atom替换换行符
直接[Ctrl]+[F],然后选择正则,输入\n
- Java中Cookie常用操作类(Spring中操作Cookie)
说明:Cookie下用Key取值没有快速的方法,只能便利循环去取. 技巧:置0则cookie会立即删除,设置-1,负值则会在关闭浏览器后删除.切记一定要增加路径:setPath("/&quo ...
- html5 拖拽元素
利用html5实现元素的拖拽,让拖动元素在指定的容器中拖动. 注意点:1.被拖元素要设置拖动属性.draggable="true" 2.容器元素要设置阻止默认事件处理 实现效果图如 ...
- 使用python分析解压zip、jar包等
python内置zipfile import zipfile, os zipFile = zipfile.ZipFile(os.path.join(os.getcwd(), 'txt.zip')) f ...
- 【Nginx】定时器事件
转自:烟雨江南 Nginx事件管理主要是网络事件和定时器事件.下面介绍定时器事件管理,即超时管理. 为什么进行超时管理? Nginx有必要对可能发生超时的事件 进行统一管理,并在事件超时时作出相应的处 ...
- 使用HDP快速搭建Hadoop开发环境 | Debugo
本文简单记录了一下使用VMware workstation 10.CentOS和HDP 2.0.6(Hadoop 2.2)发行版构建Hadoop开发测试环境的全部流程.这个过程中我遇到了不少问题,也耽 ...