bzoj1837: [CROATIAN2009]cavli 凸包1
Description
Input
Output
预处理删点顺序,逆序加点,用set和链表维护凸包上的点的极角序(以最初凸包内一点为原点),每次更新凸包时计算面积的变化量。
#include<bits/stdc++.h>
typedef long long i64;
int _(){
int x=,c=getchar();
while(c<)c=getchar();
while(c>)x=x*+c-,c=getchar();
return x;
}
int n;
double xs=,ys=;
struct _pos{
int id;
double a;
bool operator<(const _pos&w)const{return a<w.a;}
};
std::set<_pos>ss;
struct pos{
int x,y,id;
_pos get(){return(_pos){id,std::atan2(y-ys,x-xs)};}
}ps[],ps2[],ps0[];
bool cx(const pos&a,const pos&b){return a.x<b.x;}
bool cy(const pos&a,const pos&b){return a.y<b.y;}
char s[];
bool del[];
int ds[],dp=,nx[],pv[],cs[],ap=;
i64 ans=,as[];
bool chk(int a,int b,int c){
return i64(ps0[a].x-ps0[b].x)*(ps0[c].y-ps0[b].y)-i64(ps0[c].x-ps0[b].x)*(ps0[a].y-ps0[b].y)>=;
}
void cal(int a,int b,int c){
ans+=i64(ps0[c].x-ps0[b].x)*(ps0[a].y-ps0[b].y)-i64(ps0[a].x-ps0[b].x)*(ps0[c].y-ps0[b].y);
}
void dels(int w){
cal(nx[w],w,pv[w]);
nx[pv[w]]=nx[w];
pv[nx[w]]=pv[w];
ss.erase(ss.find(ps0[w].get()));
}
void inss(int l,int a,int r,_pos w){
cal(l,a,r);
nx[pv[a]=l]=pv[nx[a]=r]=a;
ss.insert(w);
}
int main(){
n=_();
for(int i=;i<n;++i){
ps[i].x=_();
ps[i].y=_();
ps[i].id=i;
}
memcpy(ps2,ps,sizeof(pos)*n);
memcpy(ps0,ps,sizeof(pos)*n);
std::sort(ps,ps+n,cx);
std::sort(ps2,ps2+n,cy);
int l1=,r1=n-,l2=,r2=n-;
scanf("%s",s);
for(int i=;s[i];++i){
if(s[i]=='L'){
while(del[ps[l1].id])++l1;
del[ds[dp++]=ps[l1++].id]=;
}else
if(s[i]=='R'){
while(del[ps[r1].id])--r1;
del[ds[dp++]=ps[r1--].id]=;
}else
if(s[i]=='D'){
while(del[ps2[l2].id])++l2;
del[ds[dp++]=ps2[l2++].id]=;
}else{
while(del[ps2[r2].id])--r2;
del[ds[dp++]=ps2[r2--].id]=;
}
}
int cp=;
for(int i=l1;i<=r1;++i)if(!del[ps[i].id]){
xs+=ps[i].x;
ys+=ps[i].y;
cs[cp++]=ps[i].id;
}
xs/=;ys/=;
ss.insert(ps0[cs[]].get());
ss.insert(ps0[cs[]].get());
pv[cs[]]=nx[cs[]]=cs[];
pv[cs[]]=nx[cs[]]=cs[];
while(dp){
int a=ds[--dp];
_pos w=ps0[a].get();
std::set<_pos>::iterator it=ss.lower_bound(w);
if(it==ss.end())it=ss.begin();
int r=it->id,l=pv[r];
if(!chk(l,a,r)){
while(chk(a,r,nx[r]))dels(r),r=nx[r];
while(chk(pv[l],l,a))dels(l),l=pv[l];
inss(l,a,r,w);
}
as[ap++]=ans;
}
while(ap){
ans=as[--ap];
printf("%lld.%lld\n",ans>>,*(ans&));
}
return ;
}
bzoj1837: [CROATIAN2009]cavli 凸包1的更多相关文章
- bzoj AC倒序
Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...
- [poj1113][Wall] (水平序+graham算法 求凸包)
Description Once upon a time there was a greedy King who ordered his chief Architect to build a wall ...
- ZOJ 3871 Convex Hull(计算几何、凸包)
题意:给n个点,|x[i]|,|y[i]| <= 1e9.求在所有情况下的子集下(子集点数>=3),凸包的面积和. 这题主要有几个方面,一个是凸包的面积,可以直接用线段的有向面积和求得,这 ...
- UVALive 2453 Wall (凸包)
题意:给你一个多边形的城堡(多个点),使用最短周长的城墙将这个城堡围起来并保证城墙的每个点到城堡上的每个点的距离都不小于l 题解:因为两点间的直线一定比折线短,所以这样做 先使用所有点求得一个凸包,接 ...
- UVA 11168 Airport(凸包+直线方程)
题意:给你n[1,10000]个点,求出一条直线,让所有的点都在都在直线的一侧并且到直线的距离总和最小,输出最小平均值(最小值除以点数) 题解:根据题意可以知道任意角度画一条直线(所有点都在一边),然 ...
- 关于2016.12.12——T1的反思:凸包的意义与应用
2016.12.12 T1 给n个圆,保证圆圆相离,求将圆围起来的最小周长.n<=100 就像上图.考场上,我就想用切线的角度来做凸包.以圆心x,y排序,像点凸包一样,不过用两圆之间的下切线角度 ...
- poj1228--稳定凸包
题目大意:给你一个凸包上的某些点(可能在凸包内),询问是否能确定这个凸包. 思路:先求出题目给出的点的凸包,看看在凸包的每条边内(不包括端点)有没有点,若有,则这条边是确定的,若没有,则这条边不确定, ...
- POJ 2225 / ZOJ 1438 / UVA 1438 Asteroids --三维凸包,求多面体重心
题意: 两个凸多面体,可以任意摆放,最多贴着,问他们重心的最短距离. 解法: 由于给出的是凸多面体,先构出两个三维凸包,再求其重心,求重心仿照求三角形重心的方式,然后再求两个多面体的重心到每个多面体的 ...
- HDU 4946 Area of Mushroom(构造凸包)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4946 题目大意:在一个平面上有n个点p1,p2,p3,p4....pn,每个点可以以v的速度在平面上移 ...
随机推荐
- CentOS 6.x 一键安装PPTP VPN脚本
环境 CentOS 6.x 32位/64位XEN/KVM/OpenVZ 步骤 依次运行下列命令 #wget http://www.hi-vps.com/shell/vpn_centos6.sh #ch ...
- CSS布局:两列,左边宽度自适应,右边宽度固定200px
<div id="box1"> <div id="left1"></div> <div id="right1 ...
- jquery异步加载json格式的数据
1.直接使用$.getJSON()方法是加载不了与静态界面同级别的本地的json后缀的文件. 2.解决办法:将json后缀的文件改为js后缀,这样就相当于加载了一个js文件. 解决办法:用$.getS ...
- CodeForces #367 div2 C
题目链接: Hard problem 题意:每个字符串可以选择反转或者不反转,给出反转每个字符串的代价,问使最少的代价使得这个字符串序列成字典序. dp[i][j] = x : 第一维是第i个字符串, ...
- iOS开发网络篇—数据安全
iOS开发网络篇—数据安全 一.简单说明 1.说明 在开发应用的时候,数据的安全性至关重要,而仅仅用POST请求提交用户的隐私数据,还是不能完全解决安全问题. 如:可以利用软件(比如Charles)设 ...
- INSTALL_FAILED_INSUFFICIENT_STORAGE
现象:运行程序,进行安装时,ANDROID模拟器启动失败,在Eclipse的控制台里log显示如下错误信息 Installation error: INSTALL_FAILED_INSUFFI ...
- oracle 11g ORA-12541: TNS: 无监听程序 (DBD ERROR: OCIServerAttach)
From :http://www.cnblogs.com/wangyt223/archive/2012/12/11/2812931.html em无法浏览,同时监听起不来.同时他的监听服务还是正常的, ...
- linux 2.6.21版本的内核合法的MAC地址
当执行ifconfig eth0 hw ether 11:22:33:44:55:66时,当前内核显示修改成功,但是ping时只无限发送ARP包,PC机也已经给板子回ARP包,但没有任何ICMP包的信 ...
- matlab 非平稳变化时域分析
对于非平稳信号,由于傅立叶变换核心函数-正弦函数具有无限性,因此选用短时窗来分析局域信号: 需要注意的时,选取完滑动的时间窗一般是中心对称而且为奇数,这时被分析的时间点正好是滑动窗的中点. 因此,时域 ...
- Python开发入门与实战7-Django Form
7. Django Form 7.1. Form表单 Django带有一个form库,称为django.forms,这个库可以处理上一章提到的包括HTML表单的自动生成以及数据验证. 我们在inven ...