Luogu-2521 [HAOI2011]防线修建
倒过来处理所有询问,就变成了一道动态凸包的裸题
吐槽一下这道题只要维护上凸壳就好了,我zz了没好好看题打了两个2333
// luogu-judger-enable-o2
#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define rp (*r)
#define lp (*l)
#define rrp (*rr)
#define llp (*ll)
using namespace std;
const int maxn=1e6+100;
struct Vector{
double x,y;
Vector(double xx=0,double yy=0){
x=xx,y=yy;
}
};
struct Point{
double x,y;
Point(double xx=0,double yy=0){
x=xx,y=yy;
}
bool operator < (Point a) const{
return x==a.x?y<a.y:x<a.x;
}
}a[maxn],b[maxn],tmp[maxn];
int dcmp(double x){return fabs(x)<1e-9?0:x>0;}
bool operator == (Point a,Point b){return a.x==b.x&&a.y==b.y;}
Vector operator - (Point a,Point b){return Vector(a.x-b.x,a.y-b.y);}
double operator * (Vector a,Vector b){return a.x*b.y-a.y*b.x;}
double dot(Vector a,Vector b){return a.x*b.x+a.y*b.y;}
double len(Vector a){return sqrt(dot(a,a));}
int m,s,t,ms[maxn],err[maxn],c[maxn],n;
double pp,x,y,ans,an[maxn];
set<Point>u,d;
void tb(Point *p,int &n){
if(n==1) return;
sort(p+1,p+n+1);
tmp[1]=p[1];
int m=1;
for(int i=2;i<=n;i++){
while(m>1&&(tmp[m]-tmp[m-1])*(p[i]-tmp[m-1])<=0)
m--;
tmp[++m]=p[i];
}
for(int i=1;i<m;i++)
d.insert(tmp[i]),ans+=len(tmp[i+1]-tmp[i]);
d.insert(tmp[m]);
int k=m;
for(int i=n-1;i>=1;i--){
while(m>k&&(tmp[m]-tmp[m-1])*(p[i]-tmp[m-1])<=0)
m--;
tmp[++m]=p[i];
}
for(int i=k;i<m;i++)
u.insert(tmp[i]),ans+=len(tmp[i+1]-tmp[i]);
u.insert(tmp[m]);
ans-=pp;
}
void uinsert(Point p){
set<Point>::iterator l,r,ll,rr;
r=u.lower_bound(p);
l=--r,r++;
// printf("%lf %lf\n%lf %lf\n",lp.x,lp.y,rp.x,rp.y);
if((p-lp)*(rp-lp)>0) return;
ans=ans-len(rp-lp)+len(p-lp)+len(p-rp);
while(l!=u.begin()){
ll=--l,l++;
// printf("%lf %lf\n",llp.x,llp.y);
if((lp-llp)*(p-llp)<0) break;
// printf("%lf %lf %lf\n",len(p-lp),len(lp-llp),len(p-llp));
ans=ans-len(p-lp)-len(lp-llp)+len(p-llp);
u.erase(l);
l=ll;
}
rr=++r,r--;
while(rr!=u.end()){
// printf("%lf %lf\n",rrp.x,rrp.y);
if((rrp-p)*(rp-p)>0) break;
// printf("%lf %lf %lf\n",len(rp-p),len(rrp-rp),len(rrp-p));
ans=ans-len(rp-p)-len(rrp-rp)+len(rrp-p);
u.erase(r);
r=rr;
rr=++r,r--;
}
u.insert(p);
}
/*
void dinsert(Point p){
set<Point>::iterator l,r,ll,rr;
r=d.lower_bound(p);
l=--r,r++;
if((p-lp)*(rp-lp)<0) return;
ans=ans-len(rp-lp)+len(p-lp)+len(p-rp);
while(l!=d.begin()){
ll=--l,l++;
if((lp-llp)*(p-llp)>0) break;
ans=ans-len(lp-p)-len(llp-lp)+len(llp-p);
u.erase(l);
l=ll;
}
rr=++r,r--;
while(rr!=d.end()){
if((rrp-p)*(rp-p)>0) break;
ans=ans-len(rp-p)-len(rrp-rp)+len(rrp-p);
u.erase(r);
r=rr;
rr=++r,r--;
}
d.insert(p);
}
*/
int main(){
scanf("%lf%lf%lf",&pp,&x,&y);
scanf("%d",&m);
for(int i=1;i<=m;i++)
scanf("%lf%lf",&b[i].x,&b[i].y);
scanf("%d",&t);
for(int i=1;i<=t;i++){
scanf("%d",&ms[i]);
if(ms[i]==1){
scanf("%d",&c[i]);
err[c[i]]=1;
}
}
a[1].x=a[1].y=a[2].y=0;
a[2].x=pp,a[3].x=x,a[3].y=y;
n=3;
for(int i=1;i<=m;i++)
if(!err[i])
a[++n]=b[i];
tb(a,n);
for(int i=t;i>=1;i--){
if(ms[i]==1)
uinsert(b[c[i]]);//,dinsert(b[c[i]]);
else
an[i]=ans;
}
for(int i=1;i<=t;i++)
if(ms[i]==2)
printf("%.2lf\n",an[i]);
return 0;
}
Luogu-2521 [HAOI2011]防线修建的更多相关文章
- [luogu P2521] [HAOI2011]防线修建
[luogu P2521] [HAOI2011]防线修建 题目描述 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国 ...
- 【BZOJ 2300】 2300: [HAOI2011]防线修建 (动态凸包+set)
2300: [HAOI2011]防线修建 Description 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上 ...
- BZOJ 2300: [HAOI2011]防线修建( 动态凸包 )
离线然后倒着做就变成了支持加点的动态凸包...用平衡树维护上凸壳...时间复杂度O(NlogN) --------------------------------------------------- ...
- P2521 [HAOI2011]防线修建
题目链接:P2521 [HAOI2011]防线修建 题意:给定点集 每次有两种操作: 1. 删除一个点 (除开(0, 0), (n, 0), 与指定首都(x, y)) 2. 询问上凸包长度 至于为什么 ...
- bzoj千题计划236:bzoj2300: [HAOI2011]防线修建
http://www.lydsy.com/JudgeOnline/problem.php?id=2300 维护动态凸包,人懒用的set 用叉积判断,不要用斜率 #include<set> ...
- 【BZOJ2300】[HAOI2011]防线修建 set维护凸包
[BZOJ2300][HAOI2011]防线修建 Description 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可 ...
- 【题解】P2521 [HAOI2011]防线修建(动态凸包)
[题解]P2521 [HAOI2011]防线修建(动态凸包) 凸包是易插入不好删除的东西,按照剧情所以我们时光倒流 然后问题就是维护凸包的周长,支持加入 本来很简单,但是计算几何就是一些小地方经验不足 ...
- 【bzoj2300】【Luogu P2521】 [HAOI2011]防线修建 动态凸包,平衡树,Set
一句话题意:给你一个凸包,每次可以插入一个点或者询问周长. 动态凸包裸题嘛,用\(Set\)实现.最初每个点坐标做乘三处理,便于取初始三角形的重心作为凸包判定原点. #include <bits ...
- BZOJ2300[HAOI2011]防线修建——非旋转treap+凸包(平衡树动态维护凸包)
题目描述 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上层现在还犹豫不决,到底该把哪些城市作为保护对象呢?又由于 ...
- LG2521 [HAOI2011]防线修建
题意 题目描述 近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了.可是A国上层现在还犹豫不决,到底该把哪些城市作为保护对象呢? ...
随机推荐
- svn hooks使用
最近要将某个目录做samba共享出去,而想通过svn同步文档到svn,然后通过svn hooks 同步到共享目录,实现自动化 现在svn服务器和samba server再同一台机器上: 在svn路径下 ...
- CentOs yum源安装 nginx
1 更新源 [root@server ~]#rpm -Uvh http://download.fedora.redhat.com/pub/epel/5/i386/epel-release-5-4.no ...
- Linux下redis安装与使用 (转)
尊重原创:https://www.cnblogs.com/codersay/p/4301677.html,并更正如下红字 redis官网地址:http://www.redis.io/ 最新版本:2.8 ...
- spark 的一些常用函数 filter,map,flatMap,lookup ,reduce,groupByKey
定义不带参数也不带返回值的函数(def :定义函数的关键字 printz:方法名称) scala> def printz = print("scala hello") ...
- JAVA问题之泛型数组
java中类似下面的代码编译器是会报错的: LinkedList<LinkedList<String>>[] li=new LinkedList<LinkedList ...
- Quartz.Net - Lesson2: 任务和触发器
Lesson 2: 任务和触发器 本系列文章是官方3.x文档的翻译,原文地址:https://www.quartz-scheduler.net/documentation/quartz-3.x/tut ...
- java -cp 命令 java jar 命令和 hadoop jar 命令
-cp 和 -classpath 一样,是指定类运行所依赖其他类的路径,通常是类库,jar包之类,需要全路径到jar包,window上分号“;” java -cp .;myClass.jar pa ...
- MHA常用命令
.查看ssh登陆是否成功 masterha_check_ssh --conf=/etc/masterha/app1.cnf .查看复制是否建立好 masterha_check_repl --conf= ...
- Mysql----mysql启动服务时提示"服务名无效"
1,首先说明一下我的环境,我刚开始是用的XAMPP这个集成的软件,里面安装了apache, mysql,tomcat这些软件,然后通过控制面板对其进行启动关闭的操作,这些操作很方便,但是我就用net ...
- 用代码构造PreferenceScreen
在PreferenceFregment中构造界面,简单省事的方法就是使用findPreference然后在xml里把UI写好.在代码中动态的添加UI内容也是需要的.核心代码是: PreferenceS ...