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

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

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
#include<set>
#include<algorithm>
#include<iostream>
#define N 500050
#define inf 2000000000
using namespace std;
int sc()
{int i=0,f=1; char c=getchar();while(c>'9'||c<'0'){if(c=='-')f=-1;c=getchar();}while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();return i*f;
}
set<pair<int,int> >st;
double ans,A[100002];
int a[100002],b[100002],d[100002];
int n,fx,fy,top,m,q;
bool c[100002];
double kk(int x,int y){return sqrt((double)x*x+y*y);}
double cal(int x,int y)
{if(!x)return inf;return (double)y/x;
}
void add(int x,int y)
{set<pair<int,int> >::iterator it,e;it=st.lower_bound(make_pair(x,y));if((*it).first==x)return;bool flag=0;while(++it!=st.end()){pair<int,int> b=*it,a=*--it;double l1=cal(a.first-x,a.second-y);double l2=cal(b.first-x,b.second-y);if(l2<l1)break;if(!flag){pair<int,int> c=*--it;it++;ans-=kk(c.first-a.first,c.second-a.second);flag=1;}ans-=kk(a.first-b.first,a.second-b.second);e=it,++it;st.erase(e);}it=st.lower_bound(make_pair(x,y));it--;while(it!=st.begin()){pair<int,int> a=*it,b=*--it;double l1=cal(a.first-x,a.second-y);double l2=cal(b.first-x,b.second-y);if(l2>l1)break;if(!flag){it++;it++;pair<int,int> c=*it;it--;it--;ans-=kk(c.first-a.first,c.second-a.second);flag=1;}ans-=kk(a.first-b.first,a.second-b.second);e=it,st.erase(++e);}if(flag) {it=st.lower_bound(make_pair(x,y));ans+=kk((*it).first-x,(*it).second-y);it--;ans+=kk((*it).first-x,(*it).second-y);st.insert(make_pair(x,y));}   else{it=st.lower_bound(make_pair(x,y));pair<int,int> a=*it,b=*--it;double l1=cal(a.first-x,a.second-y);double l2=cal(b.first-x,b.second-y);if(l1<l2){ans-=kk(a.first-b.first,a.second-b.second);ans+=kk((*it).first-x,(*it).second-y);it++;ans+=kk((*it).first-x,(*it).second-y);st.insert(make_pair(x,y));}}
}
int main()
{n=sc(),fx=sc(),fy=sc(),m=sc();for(int i=1;i<=m;i++)a[i]=sc(),b[i]=sc();q=sc();for(int i=1;i<=q;i++){d[i]=sc();if(d[i]==1)c[d[i]=sc()]=1;else d[i]=0;}st.insert(make_pair(0,0));st.insert(make_pair(n,0));st.insert(make_pair(fx,fy));ans=kk(fx,fy)+kk(fx-n,fy);for(int i=1;i<=m;i++)if(!c[i])add(a[i],b[i]);for(int i=q;i;i--){if(!d[i])A[++top]=ans;else add(a[d[i]],b[d[i]]);}while(top)printf("%.2lf\n",A[top--]);return 0;
}

BZOJ 2300: [HAOI2011]防线修建|set维护凸壳相关推荐

【BZOJ2300】[HAOI2011]防线修建 set维护凸包
[BZOJ2300][HAOI2011]防线修建 Description 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可 ...
[luogu P2521] [HAOI2011]防线修建
[luogu P2521] [HAOI2011]防线修建题目描述近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国 ...
[bzoj2300] [HAOI2011]防线修建
Description 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上层现在还犹豫不决,到底该把哪些城市作为保护 ...
Bzoj2300 / 洛谷P2521 [HAOI2011]防线修建
题目描述近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上层现在还犹豫不决,到底该把哪些城市作为保护对象呢?又由于 ...
时空旅行[线段树分治][维护凸壳]
文章目录前言题目思路代码前言肝了一上午-这是我才学线段树分治的例题-真舒服题目温馨提示:首先在UOJ做,LOJ挖数据,BZOJ终极评测... UOJ198 二手剽- 思路为什么不能用 ...
cogs 547:[HAOI2011] 防线修建
★★★☆ 输入文件:defense.in 输出文件:defense.out 简单对比时间限制:1 s 内存限制:128 MB 题目描述: 近来A国和B国的矛盾激化,为了预防不测,A国 ...
P2521 [HAOI2011]防线修建
传送门一眼看出可以倒着做转为加点维护凸包,然后--然后我就不会了-- 看了一眼题解,大概是这样的,我们先把所有点都读进来,然后按极角排序,也就是说定义点的大小为他们极角的大小(本题里实际上直接按x坐 ...
[haoi2011]防线修建
动态加点维护凸包. 论STL的熟练运用. #include<cstdio> #include<algorithm> #include<cstring> #inclu ...
bzoj2300 [HAOI2011]防线修建离线凸包
这个题可以在线做,splay维护凸包的同时维护区间最值线段树,当删除一个点就找跨过这个点的不在凸包上的点的y最值,然后分下去由于一次确定一个点,所以复杂度还是nlogn的,只是比较难写如果一个操作 ...

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

BZOJ 2300: [HAOI2011]防线修建|set维护凸壳相关推荐

最新文章

热门文章