BZOJ 2957 楼房重建
Posted Doggu
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了BZOJ 2957 楼房重建相关的知识,希望对你有一定的参考价值。
2957: 楼房重建
Description
小A的楼房外有一大片施工工地,工地上有N栋待建的楼房。每天,这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆,数自己能够看到多少栋房子。
为了简化问题,我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置,第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示,其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交,那么这栋楼房就被认为是可见的。
施工队的建造总共进行了M天。初始时,所有楼房都还没有开始建造,它们的高度均为0。在第i天,建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大---修建,也可以比原来小---拆除,甚至可以保持不变---建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后,他能看到多少栋楼房?
Input
第一行两个正整数N,M
接下来M行,每行两个正整数Xi,Yi
Output
M行,第i行一个整数表示第i天过后小A能看到的楼房有多少栋
Sample Input
2 4
3 6
1 1000000000
1 1
Sample Output
1
1
2
对于所有的数据1<=Xi<=N,1<=Yi<=10^9,N,M<=100000
Source
原以为不过是道水题,是道线段树裸题。不就是个长度为n的序列,m个操作,每个操作单点修改(delta=y*1.0/x,斜率)然后询问root罢了。
不用pushdown,但是update确实让人很迷。我最开始以为自己的update是O(lg n)的,但最后居然成了O(n)的(逻辑不严谨)。啊啊啊啊啊!
注意到题目要求,nd->sum=nd->ls->sum+find(nd->rs,nd->ls->max);原来居然还存了一段线段的min值,结果发现并不需要。find(nd,val)指一段线段中比val大的能看见的总数。递归分解,若nd->ls->max<=nd,那么左边就注定为0了,再去find(nd->rs,val)。否则右边没有影响,则return nd->num-nd->ls->num+find(nd->ls,val),这样就可以保证每次范围减半,就是O(lg n)的。
代码如下:
1 /************************************************************** 2 Problem: 2957 3 User: Doggu 4 Language: C++ 5 Result: Accepted 6 Time:1424 ms 7 Memory:4732 kb 8 ****************************************************************/ 9 10 #include <cstdio> 11 #include <cstring> 12 #include <algorithm> 13 14 template<class T>inline void readin(T &res) { 15 static char ch;T flag=1; 16 while((ch=getchar())<‘0‘||ch>‘9‘)if(ch==‘-‘)flag=-1; 17 res=ch-48;while((ch=getchar())>=‘0‘&&ch<=‘9‘)res=(res<<1)+(res<<3)+ch-48;res*=flag; 18 } 19 20 const int N = 100100; 21 struct Node { 22 double max; 23 int num;Node *ls, *rs; 24 }pool[N<<1], *tail=pool, *root; 25 26 Node *build(int lf,int rg) { 27 Node *nd=tail++; 28 if(lf==rg) { 29 nd->max=nd->num=0; 30 } else { 31 int mid=(lf+rg)>>1; 32 nd->ls=build(lf,mid); 33 nd->rs=build(mid+1,rg); 34 nd->max=nd->num=0; 35 } 36 return nd; 37 } 38 int find(Node *nd,int lf,int rg,double delta) { 39 if(lf==rg) { 40 return nd->max>delta; 41 } 42 int mid=(lf+rg)>>1, ans=0; 43 if(nd->ls->max<=delta) return find(nd->rs,mid+1,rg,delta); 44 return nd->num-nd->ls->num+find(nd->ls,lf,mid,delta); 45 } 46 void modify(Node *nd,int lf,int rg,int pos,double delta) { 47 if(lf==rg) { 48 nd->max=delta; 49 nd->num=1;return ; 50 } 51 int mid=(lf+rg)>>1; 52 if(pos<=mid) modify(nd->ls,lf,mid,pos,delta); 53 else modify(nd->rs,mid+1,rg,pos,delta); 54 nd->num=nd->ls->num+find(nd->rs,mid+1,rg,nd->ls->max), nd->max=std::max(nd->ls->max,nd->rs->max); 55 } 56 57 int main() { 58 int n, m, x, y; 59 readin(n);readin(m); 60 root=build(1,n); 61 for( int i = 1; i <= m; i++ ) { 62 readin(x);readin(y); 63 modify(root,1,n,x,y*1.0/x); 64 printf("%d\n",root->num); 65 } 66 return 0; 67 }
不过,天算不如人算。没有想到,居然可以直接上暴力!当然,前提是要分块。每一次把修楼点所属的块暴力重构。这是有策略的,扫一遍保证单调性,是O(sqrt(n))的,之前做O(sqrt(n)log n)便会TLE。之后枚举每个块,upper_bound之前的最大值,加上并更新最值,而这是O(sqrt(n)log n)的。
我用vector来存块,因为不熟悉,所以卡了很久。调试时间比线段树久,但确实是不可多得的绝佳经历。
代码如下(似乎比线段树短):
1 /************************************************************** 2 Problem: 2957 3 User: Doggu 4 Language: C++ 5 Result: Accepted 6 Time:4052 ms 7 Memory:3976 kb 8 ****************************************************************/ 9 10 #include <cstdio> 11 #include <cmath> 12 #include <algorithm> 13 #include <vector> 14 template<class T>inline void readin(T &res) { 15 static char ch;T flag=1; 16 while((ch=getchar())<‘0‘||ch>‘9‘)if(ch==‘-‘)flag=-1; 17 res=ch-48;while((ch=getchar())>=‘0‘&&ch<=‘9‘)res=(res<<1)+(res<<3)+ch-48;res*=flag; 18 } 19 const int N = 100100; 20 int belong[N], pos, bound, num[N], tnum; 21 double aa[N], tmp; 22 std::vector<double> bb[N]; 23 int main() { 24 int n, m, S, x, y; 25 readin(n);readin(m); 26 S=(int)floor(sqrt(n*log(n)/log(2)/2)); 27 pos=1-S;for( int i = 1; i <= n; i++ ) { 28 belong[i]=pos; 29 if(i-pos>=S) pos=i, num[i]=++tnum, belong[i]=pos; 30 } 31 for( int i = 1; i <= m; i++ ) { 32 readin(x);readin(y); 33 aa[x]=y*1.0/x; 34 bound=std::min(belong[x]+S-1,n); 35 tmp=0;bb[num[belong[x]]].clear();for( int i = belong[x]; i <= bound; i++ ) if(aa[i]>tmp) bb[num[belong[x]]].push_back(aa[i]),tmp=aa[i]; 36 tmp=0;int ans=0; 37 for( int i = 1; i <= tnum; i++ ) { 38 bound=bb[i].size()-1; 39 pos=std::upper_bound(bb[i].begin(),bb[i].end(),tmp)-bb[i].begin(); 40 ans+=bound+1-pos; 41 if(bound>=0) tmp=std::max(tmp,bb[i][bound]); 42 } 43 printf("%d\n",ans); 44 } 45 return 0; 46 }
恩。最后讨论一下速度。似乎暴力只慢了3倍左右。
线段树 | 4732 kb | 1424 ms | 1523 B |
暴力分块+vector+upper_bound | 3976 kb | 4052 ms | 1232 B |
以上是关于BZOJ 2957 楼房重建的主要内容,如果未能解决你的问题,请参考以下文章