题目描述

给你(0,0)、(n,0)、(x,y)和另外m个点,除(0,0)(n,0)外每个点横坐标都大于0小于n,纵坐标都大于0。

输入

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

输出

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

样例输入

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

样例输出

6.47
5.84
4.47


题解

离线+STL-set维护凸包

很容易想到离线,然后转变为加点,维护凸壳周长——经典的动态凸包问题。

把所有凸包上的点按横坐标维护平衡树,插入一个点时,首先看它是否在凸包内。具体方法:找出其前驱后继的点,判断是否上凸。容易验证这样时正确的。

然后考虑加入这个点,需要弹掉什么样的点:左边:找该点的前驱以及前驱的前驱,判断是否上凸,不上凸则弹掉前驱,否则停止。右边同理。

由于一个点只被删除一次,因此时间复杂度时 $O(n\log n)$ 的。

判断上凸可以使用叉积来判断。

由于本题不需要在凸包上二分,因此平衡树只需要维护点坐标,使用STL-set即可。

具体还是看代码吧。

#include <set>
#include <cmath>
#include <cstdio>
#define N 100010
using namespace std;
struct data
{
int x , y;
data() {}
data(int a , int b) {x = a , y = b;}
bool operator<(const data &a)const {return x == a.x ? y < a.y : x < a.x;}
data operator-(const data &a)const {return data(x - a.x , y - a.y);}
int operator*(const data &a)const {return x * a.y - y * a.x;}
inline double calc() {return sqrt(x * x + y * y);}
}a[N];
set<data> s;
int del[N] , opt[N << 1] , v[N << 1];
double now , ans[N << 1];
inline void modify(data p)
{
data a , b;
set<data>::iterator it = s.lower_bound(p);
b = *it , a = *--it;
if((p - a) * (b - p) >= 0) return;
now -= (a - b).calc();
while(it != s.begin())
{
a = *it , b = *--it;
if((p - a) * (b - a) >= 0) now -= (a - b).calc() , s.erase(a);
else break;
}
it = s.lower_bound(p);
while(it != --s.end())
{
a = *it , b = *++it;
if((p - a) * (b - a) <= 0) now -= (a - b).calc() , s.erase(a);
else break;
}
it = s.lower_bound(p) , b = *it , a = *--it;
now += (p - a).calc() + (p - b).calc() , s.insert(p);
}
int main()
{
int k , x , y , n , m , i;
scanf("%d%d%d%d" , &k , &x , &y , &n);
s.insert(data(0 , 0)) , s.insert(data(k , 0)) , s.insert(data(x , y)) , now = data(x , y).calc() + data(x - k , y).calc();
for(i = 1 ; i <= n ; i ++ ) scanf("%d%d" , &a[i].x , &a[i].y);
scanf("%d" , &m);
for(i = 1 ; i <= m ; i ++ )
{
scanf("%d" , &opt[i]);
if(opt[i] == 1) scanf("%d" , &v[i]) , del[v[i]] = 1;
}
for(i = 1 ; i <= n ; i ++ )
if(!del[i])
modify(a[i]);
for(i = m ; i ; i -- )
{
if(opt[i] == 1) modify(a[v[i]]);
else ans[i] = now;
}
for(i = 1 ; i <= m ; i ++ )
if(opt[i] == 2)
printf("%.2lf\n" , ans[i]);
return 0;
}

【bzoj2300】[HAOI2011]防线修建 离线+STL-set维护凸包的更多相关文章

  1. bzoj千题计划236:bzoj2300: [HAOI2011]防线修建

    http://www.lydsy.com/JudgeOnline/problem.php?id=2300 维护动态凸包,人懒用的set 用叉积判断,不要用斜率 #include<set> ...

  2. 2019.02.21 bzoj2300: [HAOI2011]防线修建(set+凸包)

    传送门 题意:动态维护凸包周长. 思路: 见这篇求面积的吧反正都是一个套路. 代码: #include<bits/stdc++.h> #define int long long #defi ...

  3. BZOJ2300[HAOI2011]防线修建——非旋转treap+凸包(平衡树动态维护凸包)

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

  4. BZOJ2300: [HAOI2011]防线修建

    题目:http://www.lydsy.com/JudgeOnline/problem.php?id=2300 (我只是在发以前写过的题.. 因为题目没说强制在线,所以离线乱搞就可以了.先把点删掉然后 ...

  5. 【BZOJ2300】[HAOI2011]防线修建 set维护凸包

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

  6. 【BZOJ 2300】 2300: [HAOI2011]防线修建 (动态凸包+set)

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

  7. BZOJ 2300: [HAOI2011]防线修建( 动态凸包 )

    离线然后倒着做就变成了支持加点的动态凸包...用平衡树维护上凸壳...时间复杂度O(NlogN) --------------------------------------------------- ...

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

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

  9. P2521 [HAOI2011]防线修建

    题目链接:P2521 [HAOI2011]防线修建 题意:给定点集 每次有两种操作: 1. 删除一个点 (除开(0, 0), (n, 0), 与指定首都(x, y)) 2. 询问上凸包长度 至于为什么 ...

随机推荐

  1. WPF MVVM从入门到精通4:命令和事件

    原文:WPF MVVM从入门到精通4:命令和事件   WPF MVVM从入门到精通1:MVVM模式简介 WPF MVVM从入门到精通2:实现一个登录窗口 WPF MVVM从入门到精通3:数据绑定 WP ...

  2. #ifdef 支持Mac #ifndef 支持Windows #if defined (Q_OS_WIN) 应该可以再两个系统通用

    //mac qt可以运行 #ifdef Q_OS_MAC qDebug()<<QSysInfo::MacintoshVersion; #endif //Mac不运行 #ifndef Q_O ...

  3. Mybaits: MyBaits的xml文件中大于号和小于号的转义

    < 小于号  <     > 大于号  & & 和 & ' 单引号 &apos; " 双引号  "

  4. DSP28335做FFT傅里叶变换

    1. 看了一下例程,居然没有FFT的例程,难道这个DSP28335不能做FFT吗?对了C2000系列是有C2000 ware这个库的.方便很多,不过目前不确定在C5000上运行的FFT能直接迁移到DS ...

  5. python 内置模块(hash lib)

    用于加密相关的操作,代替了md5模块和sha模块,主要提供 SHA1, SHA224, SHA256, SHA384, SHA512 ,MD5 算法 MD5 import hashlib hash=h ...

  6. [cogs2314][HZOI 2015] Persistable Editor - 可持久化平衡树

    [cogs2314][HZOI 2015]Persistable Editor - 可持久化平衡树 题目链接 首先吐槽扯淡几句 [题目描述] 维护一种可持久化的文本编辑器,支持下列操作: 1 p st ...

  7. ConfigurationProperties cannot be resolved to a type

    pom.xml 中报错之前: <parent> <groupId>org.springframework.boot</groupId> <artifactId ...

  8. windows下如何将Python文件打包成.exe可执行文件

    在使用Python做开发的时候,时不时会给自己编写了一些小工具辅助自己的工作,但是由于开发依赖环境问题,多数只能在自己电脑上运行,拿到其它电脑后就没法运行了.这显得很不方便,不符合我们的初衷,那么有没 ...

  9. pytest使用笔记(三)——pytest+allure+jenkins配置使用

    按照pytest使用笔记(二)把pytest+allure配置好后,现在在jenkins配置好,先实现手动构建(立个小目标) 一,安装jenkins插件 首页->系统管理->插件管理,从“ ...

  10. google::protobuf 编译方法

    这两天用了一下Protobuf 感觉很方便, 记录一下编译过程, 以做务忘(需要安装cmake): 1: 下载地址: https://developers.google.com/protocol-bu ...