bzoj2957 楼房重建 线段树

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);
    }
}