BZOJ2618凸多边形
来源:互联网 发布:正规网络彩票关注网站 编辑:程序博客网 时间:2024/05/18 17:05
obviously这是一道半平面交模板题目
有的dalao最开始会对其合理性产生疑惑,BUT它是逆时针给点,所以你可以用半平面交搞一搞
于是就水过去了
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<ctime>
#include<cmath>
#include<algorithm>
#include<cctype>
#include<iomanip>
#define eps 1e-8
using namespace std;
inline int read(){
int i=0,f=1;
char ch;
for(ch=getchar();ch>'9'||ch<'0';ch=getchar())
if(ch=='-') f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())
i=(i<<3)+(i<<1)+(ch^48);
return i*f;
}
int n,cnt,tot;
double ans;
struct P{
double x;
double y;
}p[3000],a[3000];
P vec(P a,P b){
P s;
s.x=a.x-b.x;
s.y=a.y-b.y;
return s;
}
double xplus(P a,P b){
return a.x*b.y-a.y*b.x;
}
struct L{
P a,b;
double slope;
friend bool operator <(const L &a,const L &b){
if(a.slope!=b.slope) return a.slope<b.slope;
return xplus(vec(a.b,a.a),vec(b.b,a.a))>0;//想一想,为什么
}
}l[3000],que[3000];
P inter(L a,L b){//依然叉乘求交点
double k1,k2,t;
k1=xplus(vec(b.b,a.a),vec(a.b,a.a));
k2=xplus(vec(a.b,a.a),vec(b.a,a.a));
t=k1/(k1+k2);
P ret;
ret.x=b.b.x+(b.a.x-b.b.x)*t;
ret.y=b.b.y+(b.a.y-b.b.y)*t;
return ret;
}
bool check(L a,L b,L t){
P p=inter(a,b);
return xplus(vec(t.b,t.a),vec(p,t.a))<0;
}
void hpi(){
sort(l+1,l+cnt+1);
int ll=1,rr=0;tot=0;
for(int i=1;i<=cnt;++i){
if(l[i].slope!=l[i-1].slope) ++tot;
l[tot]=l[i];
}
cnt=tot;tot=0;
que[++rr]=l[1];
que[++rr]=l[2];
for(int i=3;i<=cnt;++i){
while(ll<rr&&check(que[rr-1],que[rr],l[i])) --rr;//弹一弹
while(ll<rr&&check(que[ll+1],que[ll],l[i])) ++ll;//弹一弹
que[++rr]=l[i];
}
while(ll<rr&&check(que[rr-1],que[rr],que[ll])) --rr;//弹一弹
while(ll<rr&&check(que[ll+1],que[ll],que[rr])) ++ll//弹一弹;
que[rr+1]=que[ll];
for(int i=ll;i<=rr;++i)
a[++tot]=inter(que[i],que[i+1]);
return;
}
void getans(){
if(tot<3) return;
a[++tot]=a[1];
for(int i=1;i<=tot;++i)
ans+=xplus(a[i],a[i+1]);//当然是叉积爆算求面积
ans=fabs(ans)/2;
}
int main(){
n=read();
ans=0.0000;
for(int i=1;i<=n;++i){
int k=read();
for(int j=1;j<=k;++j)
p[j].x=read(),p[j].y=read();
p[k+1]=p[1];
for(int j=1;j<=k;++j)
l[++cnt].a=p[j],l[cnt].b=p[j+1];
}
for(int i=1;i<=cnt;++i)
l[i].slope=atan2(l[i].b.y-l[i].a.y,l[i].b.x-l[i].a.x);
hpi();
getans();
printf("%.3lf",ans);
return 0;
}
- BZOJ2618凸多边形
- 【CQOI2006】【BZOJ2618】凸多边形
- bzoj2618凸多边形面积交
- bzoj2618【CQOI2006】凸多边形
- bzoj2618: [Cqoi2006]凸多边形
- [BZOJ2618]CQOI2006凸多边形|扫描线
- bzoj2618 凸多边形 半平面交
- BZOJ2618 凸多边形 半平面交
- [BZOJ2618][Cqoi2006]凸多边形(半平面交)
- 半平面交 模板 bzoj2618 【Cqoi2006】凸多边形
- bzoj2618 [Cqoi2006]凸多边形(半平面+S)
- 【BZOJ2618】【Cqoi2006】凸多边形 半平面交 、算法的深度细节剖析。
- 凸多边形
- BZOJ2618 (半平面交)
- 划分凸多边形
- 凸多边形交
- 【CQOI2006】凸多边形
- hiho1360 凸多边形
- PTA习题 5-4 是否同一棵二叉搜索树(二叉搜索树的构建以及相同树的判断)
- 网易游戏2017实习生招聘-小结
- Java的SAX解析代码二
- hadoop学习总结
- 如何把自己的项目上传到GitHub上
- BZOJ2618凸多边形
- 读书笔记之--Java虚拟机精讲第六章内存分配与垃圾回收
- 一起学python(四)
- Delphi XE2 and C++ Builder XE2 Update4 HotFix1
- JUnit4单元测试报错问题 :method initializationerror not found
- 通用Adapter的实现套路(ListView篇)
- 【JSOI2007】动态最值 Splay
- poj3067 Japan 树状数组 逆序数
- 浅谈css中浮动和清除浮动带来的影响