首页 > 代码库 > [bzoj2957][楼房重建] (线段树)
[bzoj2957][楼房重建] (线段树)
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
3 42 43 61 10000000001 1
Sample Output
1112
数据约定
对于所有的数据1<=Xi<=N,1<=Yi<=10^9
N,M<=100000
HINT
Source
中国国家队清华集训 2012-2013 第一天
Solution
对于每次修改后的区间,从上往下修改再从下往上回溯
更新区间斜率最大值,对于每个区间来说,只计算左区间对右区间的遮盖作用
若左边最大值大于右边,则右边贡献为0
否则,递归处理即可
#include <stdio.h>#include <string.h>#define MaxN 100010#define MaxBuf 1<<22#define RG register#define inline __inline__ __attribute__((always_inline))#define Blue() ((S == T && (T=(S=B)+fread(B,1,MaxBuf,stdin),S == T)) ? 0 : *S++)#define mid ((x>>1)+(y>>1)+(x&y&1))#define dmax(a,b) ((a)>(b)?(a):(b))char B[MaxBuf],*S=B,*T=B;inline void Rin(RG int &x) { x=0;RG int c=Blue(),f=1; for(; c<48||c>57; c=Blue()) if(c==45)f=-1; for(; c>47&&c<58; c=Blue()) x=(x<<1)+(x<<3)+c-48; x*=f; }double tr[MaxN<<2];int n,m,ans,sum[MaxN<<2];int calc(RG int x,RG int y,RG int k,RG double val) { if(x == y)return tr[k] > val; if(tr[k<<1] < val) return calc(mid+1,y,k<<1|1,val); return sum[k]-sum[k<<1]+calc(x,mid,k<<1,val);}void modify(RG int x,RG int y,RG int k,RG int pos,RG double val) { if(x == y) { tr[k]=val; sum[k]=1; return; } pos <= mid ? modify(x,mid,k<<1,pos,val) : modify(mid+1,y,k<<1|1,pos,val); tr[k]=dmax(tr[k<<1],tr[k<<1|1]); sum[k]=sum[k<<1]+calc(mid+1,y,k<<1|1,tr[k<<1]);}int main() { Rin(n),Rin(m); while(m--) { RG int pos,x; Rin(pos),Rin(x); modify(1,n,1,pos,(double)x/pos); printf("%d\n",sum[1]); } fclose(stdin); return 0; }
[bzoj2957][楼房重建] (线段树)
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。