【BZOJ2300】[HAOI2011]防线修建

Description

近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了。可是A国上层现在还犹豫不决,到底该把哪些城市作为保护对象呢?又由于A国的经费有限,所以希望你能帮忙完成如下的一个任务:
1.给出你所有的A国城市坐标
2.A国上层经过讨论,考虑到经济问题,决定取消对i城市的保护,也就是说i城市不需要在防线内了
3.A国上层询问对于剩下要保护的城市,修建防线的总经费最少是多少
你需要对每次询问作出回答。注意单位1长度的防线花费为1。
A国的地形是这样的,形如下图,x轴是一条河流,相当于一条天然防线,不需要你再修建
A国总是有两个城市在河边,一个点是(0,0),一个点是(n,0),其余所有点的横坐标均大于0小于n,纵坐标均大于0。A国有一个不在(0,0)和(n,0)的首都。(0,0),(n,0)和首都这三个城市是一定需要保护的。

上图中,A,B,C,D,E点为A国城市,且目前都要保护,那么修建的防线就会是A-B-C-D,花费也就是线段AB的长度+线段BC的长度+线段CD的长度,如果,这个时候撤销B点的保护,那么防线变成下图

Input

第一行,三个整数n,x,y分别表示河边城市和首都是(0,0),(n,0),(x,y)。
第二行,一个整数m。
接下来m行,每行两个整数a,b表示A国的一个非首都非河边城市的坐标为(a,b)。
再接下来一个整数q,表示修改和询问总数。
接下来q行每行要么形如1 i,要么形如2,分别表示撤销第i个城市的保护和询问。

Output

对于每个询问输出1行,一个实数v,表示修建防线的花费,保留两位小数

Sample Input

4 2 1
2
1 2
3 2
5
2
1 1
2
1 2
2

Sample Output

6.47
5.84
4.47

HINT

m<=100000,q<=200000,n>1
所有点的坐标范围均在10000以内, 数据保证没有重点

题解:容易想到离线,然后本题就变成了加点和维护凸包,可以用平衡树来维护凸包来搞(以前一直以为必须用Splay,结果发现好像set就行)。

具体地,我们在set中的点是按x值从小到大排序的。在加入一个点时,先找到凸包中它的前驱和后继节点,用叉积判一下当前点是否已经在凸包里面。如果在则不用加点,否则加入当前点,并不断执行以下操作:

如果当前点和左侧两个点形成的不是上凸包,则将当前点的前驱删掉,否则退出。

右面同理。

然后就做完了,细节也不是特别多。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <set>
#include <cmath>
using namespace std;
const int maxn=200010;
struct point
{int x,y;point() {}point(int a,int b) {x=a,y=b;}bool operator != (const point &a) const {return (x!=a.x)||(y!=a.y);}bool operator < (const point &a) const {return (x==a.x)?(y<a.y):(x<a.x);}point operator - (const point &a) const {return point (x-a.x,y-a.y);}int operator * (const point &a) const {return x*a.y-y*a.x;}
}p[maxn];
double dis(const point &a,const point &b) {return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
set<point> s;
set<point>::iterator it;
double ans;
int n,m,q;
int del[maxn],op[maxn],qi[maxn];
double qa[maxn];
inline int rd()
{int ret=0,f=1;   char gc=getchar();while(gc<'0'||gc>'9')  {if(gc=='-')    f=-f;  gc=getchar();}while(gc>='0'&&gc<='9')  ret=ret*10+(gc^'0'),gc=getchar();return ret*f;
}
inline void add(point x)
{point a,b,c;it=s.upper_bound(x),b=*it,it--,a=*it;if(x.x==b.x) return ;if((x-a)*(b-x)>=0)  return ;ans-=dis(a,b);while(a!=point(0,0)){it--,c=*it;if((a-c)*(x-a)>=0) ans-=dis(a,c),s.erase(s.find(a)),a=c;else break;}it=s.find(b);while(b!=point(n,0)){it++,c=*it;if((b-x)*(c-b)>=0) ans-=dis(b,c),s.erase(s.find(b)),b=c;else break;}ans+=dis(a,x)+dis(x,b);s.insert(x);
}
int main()
{int i,a,b;n=rd(),a=rd(),b=rd(),m=rd();s.insert(point(0,0)),s.insert(point(n,0)),s.insert(point(a,b));ans=dis(point(0,0),point(a,b))+dis(point(n,0),point(a,b));for(i=1;i<=m;i++)  p[i].x=rd(),p[i].y=rd();q=rd();for(i=1;i<=q;i++){op[i]=rd();if(op[i]==1)   qi[i]=rd(),del[qi[i]]=1;}for(i=1;i<=m;i++) if(!del[i]) add(p[i]);for(i=q;i>=1;i--){if(op[i]==2) qa[i]=ans;else add(p[qi[i]]);}for(i=1;i<=q;i++) if(op[i]==2)  printf("%.2lf\n",qa[i]);return 0;
}

转载于:https://www.cnblogs.com/CQzhangyu/p/7965400.html

【BZOJ2300】[HAOI2011]防线修建 set维护凸包相关推荐

  1. bzoj2300 [HAOI2011]防线修建 离线凸包

    这个题可以在线做,splay维护凸包的同时维护区间最值线段树,当删除一个点就找跨过这个点的不在凸包上的点的y最值,然后分下去 由于一次确定一个点,所以复杂度还是nlogn的,只是比较难写 如果一个操作 ...

  2. [bzoj2300] [HAOI2011]防线修建

    Description 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上层现在还犹豫不决,到底该把哪些城市作为保护 ...

  3. BZOJ 2300: [HAOI2011]防线修建|set维护凸壳

    因为只支持插入操作,所以需要倒过来做.. 第一次写维护凸壳的东西-不抄题解..自己YY斜率乱搞(都忘了写叉乘)代码又长常数又大 1A一刻赛艇~ 还是找一个靠谱的的模板吧... #include< ...

  4. [luogu P2521] [HAOI2011]防线修建

    [luogu P2521] [HAOI2011]防线修建 题目描述 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国 ...

  5. P2521 [HAOI2011]防线修建

    传送门 一眼看出可以倒着做转为加点维护凸包,然后--然后我就不会了-- 看了一眼题解,大概是这样的,我们先把所有点都读进来,然后按极角排序,也就是说定义点的大小为他们极角的大小(本题里实际上直接按x坐 ...

  6. [haoi2011]防线修建

    动态加点维护凸包. 论STL的熟练运用. #include<cstdio> #include<algorithm> #include<cstring> #inclu ...

  7. 【bzoj2300】【Luogu P2521】 [HAOI2011]防线修建 动态凸包,平衡树,Set

    一句话题意:给你一个凸包,每次可以插入一个点或者询问周长. 动态凸包裸题嘛,用\(Set\)实现.最初每个点坐标做乘三处理,便于取初始三角形的重心作为凸包判定原点. #include <bits ...

  8. Bzoj2300 / 洛谷P2521 [HAOI2011]防线修建

    题目描述 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上层现在还犹豫不决,到底该把哪些城市作为保护对象呢?又由于 ...

  9. cogs 547:[HAOI2011] 防线修建

    ★★★☆   输入文件:defense.in   输出文件:defense.out   简单对比 时间限制:1 s   内存限制:128 MB 题目描述: 近来A国和B国的矛盾激化,为了预防不测,A国 ...

最新文章

  1. DOS call 中的%cd%,当前文件夹演示
  2. 数据中心胶体电池的使用寿命
  3. 《系统集成项目管理》第十三章 项目合同管理
  4. LINUX的“脏奶牛”
  5. Vue computed参数与各生命周期关系(主要是异步的时候)
  6. 调整自己--提升自己
  7. 列级触发器 SQL Server
  8. 保送清华成博士,华为 12 年搞通信,他为何如此看待 5G ?| 人物志
  9. string转int/float/double、int/float/double转string、转字符串数组的方法:stoi、stringstream、scanf、to_string、sprintf
  10. 电脑亮度多少对眼睛好_黄江办公文员学费大概是多少,黄江附近哪个电脑学校比较好一点...
  11. [2018.10.25 T1] 妹子
  12. 电商支付-使用Restful api接口集成Paypal支付方式(一)
  13. Android混淆介绍
  14. Resilient Distributed Datasets: A Fault-Tolerant Abstraction for In-Memory Cluster Computing
  15. Spring Mvc - 用户注册服务端生成默认头像实现
  16. android身高控件_Android 滑动选择身高体重控件——RulerView
  17. android自定义六边形,Android实现六边形图像
  18. 为什么我们需要新颖的身份验证方案?
  19. 数学建模----聚类分析
  20. ORACLE XE在centos平台下安装方法

热门文章

  1. php网页布局边框,用CSS来设置网页当中的边框
  2. 电脑运行adb闪退_adb+python进阶使用
  3. 【WebRTC---入门篇】(九)WebRTC网络基础:P2P/STUN/TURN/ICE
  4. 多方法接口回调_啊?Java反射遇到接口
  5. 计算机兴趣小组意义,信息技术兴趣小组开展之我见
  6. mysql+误操作怎么恢复_Mysql误操作恢复流程
  7. 用火狐录制脚本为空_功能测试——链接测试amp;脚本功能
  8. 该虚拟机似乎正在使用
  9. 工业交换机常用术语及常见知识点汇总
  10. 【渝粤教育】电大中专电商运营实操 (1)作业 题库