「BZOJ2300」[HAOI2011] 防线修建
Description
近来A国和B国的矛盾激化,为了预防不测,A国准备修建一条长长的防线,当然修建防线的话,肯定要把需要保护的城市修在防线内部了。可是A国上层现在还犹豫不决,到底该把哪些城市作为保护对象呢?又由于A国的经费有限,所以希望你能帮忙完成如下的一个任务:
- 给出你所有的A国城市坐标
- A国上层经过讨论,考虑到经济问题,决定取消对i城市的保护,也就是说i城市不需要在防线内了
- 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
数据范围:
30%的数据m<=1000,q<=1000
100%的数据m<=100000,q<=200000,n>1
所有点的坐标范围均在10000以内, 数据保证没有重点
题解
离线后问题变成维护动态凸壳,只有插入点的操作
诶用个set就可以了
插入点x,判断x是否在凸壳外,然后维护x两边的凸壳。。。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 |
#include<iostream> #include<cstdio> #include<cstring> #include<cstdlib> #include<map> #include<ctime> #include<vector> #include<set> #include<cmath> #include<algorithm> #define inf 1000000000 #define ll long long #define pa pair<int,int> using namespace std; int read() { int x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } double now,res[200005]; int T; int t1,t2; int n,m,Q; int b[200005]; bool mark[100005]; struct P{int x,y;}del[100005],a[200005]; set<P> q; double dis(P a,P b) { return sqrt((double)((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y))); } bool operator<(P a,P b) { if(a.x==b.x)return a.y<b.y; return a.x<b.x; } int operator*(P a,P b) { return a.x*b.y-a.y*b.x; } P operator-(P a,P b) { P t; t.x=a.x-b.x;t.y=a.y-b.y; return t; } void insert(int a,int b) { P x=(P){a,b}; set<P>::iterator r=q.lower_bound(x),l=r,t; l--; if((*r-*l)*(x-*l)<0)return; now-=dis(*l,*r); q.insert(x); while(1) { t=r;r++; if(r==q.end())break; if((*r-x)*(*t-x)>0)break; now-=dis(*t,*r); q.erase(t); } while(l!=q.begin()) { t=l;l--; if((*t-x)*(*l-x)>0)break; now-=dis(*t,*l); q.erase(t); } q.insert(x); l=r=t=q.find(x); l--;r++; now+=dis(*l,x)+dis(*r,x); } int main() { n=read(); q.insert((P){0,0});q.insert((P){n,0}); P bas; bas.x=read();bas.y=read();q.insert(bas); now+=dis((P){0,0},bas); now+=dis((P){n,0},bas); m=read(); for(int i=1;i<=m;i++) a[i].x=read(),a[i].y=read(); int typ,x; Q=read(); for(int i=1;i<=Q;i++) { typ=read(); if(typ==1){x=read();del[++t1]=a[x];mark[x]=1;} else b[++t2]=t1; } for(int i=1;i<=m;i++) if(!mark[i])insert(a[i].x,a[i].y); int T=t1; for(int i=t2;i;i--) { while(T>b[i]) { insert(del[T].x,del[T].y); T--; } res[i]=now; } for(int i=1;i<=t2;i++)printf("%.2lf\n",res[i]); return 0; } |
56行的代码和72行的代码重复,觉得应该删去56行的insert。
pascal选手怎么办
TAT
用treap代替