【题解】P2521 [HAOI2011]防线修建(动态凸包)
【题解】P2521 [HAOI2011]防线修建(动态凸包)
凸包是易插入不好删除的东西,按照剧情所以我们时光倒流
然后问题就是维护凸包的周长,支持加入
本来很简单,但是计算几何就是一些小地方经验不足容易WA和RE
然后代码注释里有一些经验
//@winlere
#include<iostream>
#include<set>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<cmath>
#define pf(x) ((x)*(x))
#define DEBUG(s,a) cerr<<#s" = "<<(s)<<" \n"[(a)==1]
#define getchar() (__c==__ed?(__ed=__buf+fread(__c=__buf,1,1<<18,stdin),*__c++):*__c++)
using namespace std; typedef long long ll; char __buf[1<<18],*__c=__buf,*__ed=__buf;
inline int qr(){
int ret=0,f=0,c=getchar();
while(!isdigit(c))f|=c==45,c=getchar();
while(isdigit(c)) ret=ret*10+c-48,c=getchar();
return f?-ret:ret;
}
const int maxn=1e5+5;
struct NODE{//向量 坐标用一个结构体代替就行了
ll x,y;//用ll吧 别小气
NODE operator - (NODE b)const{return {x-b.x,y-b.y};}
ll operator % (NODE b)const{return x*b.y-y*b.x;}
bool operator < (NODE b)const{return x<b.x;}
double operator ^ (NODE b)const{return sqrt(pf(x-b.x)+pf(y-b.y));}
}data[maxn],capt;
vector< pair<int,int> > ve;
double len,ans[maxn];
bool in[maxn];
set<NODE> s;
void insert(NODE x){
vector<NODE> ve;
if((*s.upper_bound(x)-x)%(x-*--s.upper_bound(x))<0) return;//先判断是否能够插入,这样以后好写一些
for(set<NODE>::iterator t=--s.upper_bound(x),temp;t!=s.begin();t=temp){
temp=prev(t);
if((x-*temp)%(*t-*temp)<0) ve.push_back(*t);//避免一边用iterator一边erase
else break;
}
reverse(ve.begin(),ve.end());
for(set<NODE>::iterator t=s.upper_bound(x),temp;t!=--s.end();t=temp){
temp=next(t);
if((*temp-x)%(*temp-*t)>0) ve.push_back(*t);//同上
else break;
}
if(ve.size()==0)//多写几个if比仔细比巧妙统一快得多,对于=0,=1都特判一下,这里只要=0特判
len+=(*s.upper_bound(x)^x)+(*--s.upper_bound(x)^x)-(*--s.upper_bound(x)^*s.upper_bound(x));
else {
auto l=*--s.find(ve.front()),r=*++s.find(ve.back());
len+=(l^x)+(r^x)-(l^ve.front())-(r^ve.back());
for(int t=1,ed=ve.size();t<ed;++t) len-=ve[t]^ve[t-1];//多写一个循环没问题
for(auto t:ve) s.erase(t);
}
s.insert(x);
}
int n,m,x,y,q;
int main(){
n=qr(); capt.x=qr(),capt.y=qr(); m=qr();
s.insert(capt); s.insert({n,0}); s.insert({0,0}); len=(capt^*s.begin())+(capt^*s.rbegin());
for(int t=1;t<=m;++t) data[t].x=qr(),data[t].y=qr();
q=qr();
for(int t=1;t<=q;++t){
int op=qr();
ve.push_back(op==1?(pair<int,int>){op,qr()}:(pair<int,int>){op,t});
if(op==1) in[ve.back().second]=1;
}
for(int t=1;t<=m;++t)
if(!in[t])
insert(data[t]);
reverse(ve.begin(),ve.end());
for(auto t:ve)
if(t.first==1) insert(data[t.second]);
else ans[t.second]=len;
reverse(ve.begin(),ve.end());
for(auto t:ve)
if(t.first==2) printf("%.2lf\n",ans[t.second]);
return 0;
}
【题解】P2521 [HAOI2011]防线修建(动态凸包)的更多相关文章
- BZOJ 2300: [HAOI2011]防线修建( 动态凸包 )
离线然后倒着做就变成了支持加点的动态凸包...用平衡树维护上凸壳...时间复杂度O(NlogN) --------------------------------------------------- ...
- 【bzoj2300】【Luogu P2521】 [HAOI2011]防线修建 动态凸包,平衡树,Set
一句话题意:给你一个凸包,每次可以插入一个点或者询问周长. 动态凸包裸题嘛,用\(Set\)实现.最初每个点坐标做乘三处理,便于取初始三角形的重心作为凸包判定原点. #include <bits ...
- [luogu P2521] [HAOI2011]防线修建
[luogu P2521] [HAOI2011]防线修建 题目描述 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国 ...
- P2521 [HAOI2011]防线修建
题目链接:P2521 [HAOI2011]防线修建 题意:给定点集 每次有两种操作: 1. 删除一个点 (除开(0, 0), (n, 0), 与指定首都(x, y)) 2. 询问上凸包长度 至于为什么 ...
- Bzoj2300 / 洛谷P2521 [HAOI2011]防线修建
题目描述 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上层现在还犹豫不决,到底该把哪些城市作为保护对象呢?又由于 ...
- 【BZOJ 2300】 2300: [HAOI2011]防线修建 (动态凸包+set)
2300: [HAOI2011]防线修建 Description 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上 ...
- 【BZOJ2300】[HAOI2011]防线修建 set维护凸包
[BZOJ2300][HAOI2011]防线修建 Description 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可 ...
- bzoj千题计划236:bzoj2300: [HAOI2011]防线修建
http://www.lydsy.com/JudgeOnline/problem.php?id=2300 维护动态凸包,人懒用的set 用叉积判断,不要用斜率 #include<set> ...
- BZOJ [HAOI2011]防线修建(动态凸包)
听说有一种很高端的东西叫动态凸包维护dp就像学一下,不过介于本人还不会动态凸包就去学了下,还是挺神奇的说,维护上下凸包的写法虽然打得有点多不过也只是维护复制黏贴的事情而已罢了. 先说下动态凸包怎么写吧 ...
随机推荐
- PHP中header头设置Cookie与内置setCookie的区别
首先声明,浏览的Cookie操作都是通过HTTP Header(俗称“Http头”) 来实现.所有的服务器与客户端之间Cookie数据传输都是通过Http请求头来操作. PHP中setCookie(函 ...
- 【codeforces 798D】Mike and distribution
[题目链接]:http://codeforces.com/contest/798/problem/D [题意] 让你选一个下标集合 p1,p2,p3..pk 使得2*(a[p1]+a[p2]+..+a ...
- Jquery常用方法汇总(转)
https://blog.csdn.net/lucky___star/article/details/87883888
- LA 4973 Ardenia (3D Geometry + Simulation)
ACM-ICPC Live Archive 三维几何,题意是要求求出两条空间线段的距离.题目难度在于要求用有理数的形式输出,这就要求写一个有理数类了. 开始的时候写出来的有理数类就各种疯狂乱套,TLE ...
- MATLAB常用函数, 常见问题
MATLAB常用函数 1.常用取整函数 round(x):四舍五入函数 floor(x) : 向下取整, 即 floor(1.2)=1, floor(1.8) = 1 ceil(x) : 向上取整, ...
- js利用select标签生成简易计算功能
html中使用select option作为运算符的承接容器,输入值,选择不同运算符,计算结果. 文章地址 https://www.cnblogs.com/sandraryan/ <!DOCTY ...
- Python--day65--母版和继承的基本使用
- HDU 2674
0 <= N<=10^9 看到这个数据范围知道常规方法肯定做不出来. 不过一想想既然是mod2009,是不是只要其中含有一个2009,那么其结果一定是0了呢 说了这里思路,就是看什么时候出 ...
- 深入理解 Embedding层的本质
继上文https://blog.csdn.net/weixin_42078618/article/details/82999906探讨了embedding层的降维效果,时隔一个月,分享一下嵌入层在NP ...
- 2018-8-10-win10-sdk-是否向下兼容
title author date CreateTime categories win10 sdk 是否向下兼容 lindexi 2018-08-10 19:16:53 +0800 2018-2-13 ...