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 ...
随机推荐
- hdu——3861 The King’s Problem
The King’s Problem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Other ...
- 学习日常笔记<day15>mysql基础
1.数据库入门 1.1数据库软件 数据库:俗称数据的仓库,方便管理数据的软件(或程序) 1.2市面上数据库软件 Oracle,甲骨文公司的产品. 当前最流行应用最广泛的数据库软件.和java语言兼容非 ...
- jfree-生成xy图
需要导入的包: import org.jfree.chart.*; import org.jfree.chart.plot.PlotOrientation; import org.jfree.data ...
- java IO与NIO
场景:IO适用于大而少,NIO适用于小而多 转载:https://www.cnblogs.com/kzfy/p/5063467.html 传统的socket IO中,需要为每个连接创建一个线程,当并发 ...
- MongoDB学习day08--mongoose预定义修饰符和getter、setter修饰符
一.mongoose预定义修饰符 lowercase. uppercase . trim var UserSchema=mongoose.Schema({ name:{ type:String, tr ...
- foobar2000使用cue文件播放时出现Unable to open item for playback (Object not found):的问题解决
如下错误: 一般是找不到APE文件导致的.解决方法如下: 1.打开APE文件,对一下路径修改即可.
- ssh 执行多条命令包含awk的用法
格式:ssh user@ip command 单条命令:ssh user@ip command1 多条命令:ssh user@ip "command1;command2" 不加双引 ...
- js中的自定义异常处理函数
1. Can I suppress JavaScript error messages? 2. Can I set up my own JavaScript error handler? 3. Can ...
- Linux 将一般的用户加入sudo组is_not_in_the_sudoers_file._This_incident_will_be_reported解决方法
在一般用户下执行sudo命令提示xxx is not in the sudoers file. This incident will be reported.解决方法: $where ...
- vue 获取当前时间 格式YYYY-MM-DD
函数封装: /** * 获取当前时间 * 格式YYYY-MM-DD */ Vue.prototype.getNowFormatDate = function() { var date = new Da ...