2957: 楼房重建
Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 2013 Solved: 962
[Submit][Status][Discuss]
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 4
2 4
3 6
1 1000000000
1 1
Sample Output
1
1
1
2
数据约定
对于所有的数据1<=Xi<=N,1<=Yi<=10^9
N,M<=100000
对线段树算是有了一个新的认识,以前基本上都是区间加减和01线段树维护…显然,这道题对线段树的运用说明了线段树实际上是十分灵活的.
首先题目描述有误,一栋楼房可见,应满足其最高点到(0,0)的连线不与其他楼房相交 .就是说这个连线的斜率大于之前的每一个斜率.
针对每一个线段,我们维护一个最大斜率cmax,和一个只针对当前线段能看到的房屋数目 的ans,再进行讨论情况.
首先第一个结论是,当前线段的ans一定包含左子树的ans.因为我们ans只针对当前线段,视线从0,0照过来左子树之前相当于是空的,所以整个线段的ans一定包含左子树的ans.
那是不是 当前线段ans=左子树ans+右子树ans呢?
不!
右子树的ans 只针对他自己而言 (加粗字体多次强调,只针对当前),是忽视了当前线段的左子树的房屋的.当时整个线段ans虽然只针对自己,但是是有这个左子树的,所以左子树的房屋有可能会遮到右子树的房屋,所以针对统计整个线段的ans,我们要递归计算右子树对**当前线段**ans的贡献,这是在修改(modify)return 的时候做的事.
第二个结论是:因为当前房屋的高度只会影响到当前房屋以后的房屋,所以我们在递归计算右子树的时候(不算左子树的原因在结论里说清楚了),就把左子树的cmax带到右子树统计一下影响,与左子树ans相加就是当前线段的,最后线段为整个1~n,答案呼之欲出.
分类讨论:
设右子树的左子树为左段,右段为右子树的右子树.
1.在递归计算右子树的时候,设左子树的cma为lc,左段的cmax为rlc,如果说lc>rlc,那么这个左段的房屋就都被遮盖了,对答案没有贡献…还是去计算这个右段看看有没有希望可以贡献ans吧.
2.但是要是lc<=rlc,那么对右段就没有影响(rlc比lc大,要影响右段也是这个左段的rlc去影响).右段对右子树的贡献就是右子树的ans-左段的ans(不直接加右段的ans是因为这是针对右子树而并非右段自己的,在前面有讲),我们再去找左段里面有多少大于lc的房屋,当然,递归处理.
可以配合代码理解一下.
/**************************************************************
Problem: 2957
User: MaxMercer
Language: C++
Result: Accepted
Time:1400 ms
Memory:11856 kb
****************************************************************/
#include<stdio.h>
inline char nc(){
static char buf[100000],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read(){
register int x=0;
register char ch=nc();
while(ch<'0'||ch>'9')ch=nc();
while(ch>='0'&&ch<='9')x=(x<<3)+(x<<1)+ch-'0',ch=nc();
return x;
}
inline double max(double x,double y){if(x>y) return x;return y;}
struct node{
node *ls,*rs;
double cmax;
int l,r,ans;
}pool[400005],*tail=pool,*root;
int n,m;
node *build(int lf,int rg){
node *bt=++tail;
if(lf==rg) {bt->l=bt->r=lf;return bt;}
int mid=(lf+rg)>>1;
bt->l=lf,bt->r=rg;
bt->ls=build(lf,mid);
bt->rs=build(mid+1,rg);
return bt;
}
int calc(node *bt,double val){
int lf=bt->l,rg=bt->r;
if(lf==rg) return bt->cmax>val;
if(bt->ls->cmax<=val) return calc(bt->rs,val);
return bt->ans-bt->ls->ans+calc(bt->ls,val);
}
void modify(node *bt,int lf,int rg,int pos,double val){
if(lf==rg) {bt->cmax=val,bt->ans=1;return;}
int mid=(lf+rg)>>1;
if(pos<=mid) modify(bt->ls,lf,mid,pos,val);
else modify(bt->rs,mid+1,rg,pos,val);
bt->cmax=max(bt->ls->cmax,bt->rs->cmax);
bt->ans=bt->ls->ans+calc(bt->rs,bt->ls->cmax);
}
int main(){
n=read(),m=read();
root=build(1,n);
while(m--){
int x=read(),y=read();
modify(root,1,n,x,(double)y/x);
printf("%d\n",root->ans);
}
}