[bzoj2957]: 楼房重建——线段树

题目大意：

　　小A的楼房外有一大片施工工地，工地上有N栋待建的楼房。每天，这片工地上的房子拆了又建、建了又拆。他经常无聊地看着窗外发呆，数自己能够看到多少栋房子。
　　为了简化问题，我们考虑这些事件发生在一个二维平面上。小A在平面上(0,0)点的位置，第i栋楼房可以用一条连接(i,0)和(i,Hi)的线段表示，其中Hi为第i栋楼房的高度。如果这栋楼房上任何一个高度大于0的点与(0,0)的连线没有与之前的线段相交，那么这栋楼房就被认为是可见的。
　　施工队的建造总共进行了M天。初始时，所有楼房都还没有开始建造，它们的高度均为0。在第i天，建筑队将会将横坐标为Xi的房屋的高度变为Yi(高度可以比原来大—修建，也可以比原来小—拆除，甚至可以保持不变—建筑队这天什么事也没做)。请你帮小A数数每天在建筑队完工之后，他能看到多少栋楼房？

思路：

这好像也可以说是线段树维护区间最长连续上升的模板题。。。我们用斜率来表示高度，用一个线段树来维护某个区间内的高度最大值和最长连续上升。修改高度的时候只需要普通的线段树更新就好了，但是我们需要更新最长的长度，考虑一个区间[l,r],它的lc和rc的信息都是已知的，那么就是求rc这个区间内大于Max[lc]的连续上升的数的个数。
现在我们开始递归修改，记h为Max[lc]。
1、若右区间的左子区间最大值大于等于h,那么它的右子区间的个数就可以直接求出来，继续递归lc的左子区间即可。
2、若右区间的左子区间最大值小于h，那么他的左子区间的贡献必定为0，所以继续递归右子区间即可。
总的时间复杂度$O(n*\log^2n)$(比我的分块不知道快了多少。。。)

// luogu-judger-enable-o2
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<cstdio>
using namespace std;
void File(){
    freopen("bzoj2957.in","r",stdin);
    freopen("bzoj2957_Segment_tree.out","w",stdout);
}
template<typename T>
void read(T &x){
    T _=0;char __=getchar();
    while(!isdigit(__))__=getchar();
    while(isdigit(__))_=(_<<1)+(_<<3)+(__^'0'),__=getchar();
    x=_;
}
template<typename T>
void write(T x,char c){
    if(x==0){
        putchar('0');
        putchar(c);
        return;
    }
    T y=10;
    int len=1;
    while(y<=x)y=(y<<1)+(y<<3),++len;
    while(len--){
        y=y/10;
        putchar(x/y+48);
        x=x%y;
    }
    putchar(c);
}
template<typename T>void chckmax(T &_,T __){_=_>__ ? _ : __;}
template<typename T>void chckmin(T &_,T __){_=_<__ ? _ : __;}
#define REP(i,a,b) for(register int i=a;i<=b;++i)
#define DREP(i,a,b) for(register int i=a;i>=b;--i)
#define MREP(i,x) for(register int i=beg[x];i;i=E[i].last)
#define ll long long
const int maxn=1e5+10;
int n,m;
struct Segment_Tree{
#define mid ((l+r)>>1)
#define lson rt<<1,l,mid
#define rson (rt<<1)|1,mid+1,r
    double Max[maxn<<2];
    int cnt[maxn<<2];
    int cal(int rt,int l,int r,double h){
        if(l==r)return Max[rt]>h ? 1 : 0;
        if(Max[rt<<1]<=h)return cal(rson,h);
        return cal(lson,h)+cnt[rt]-cnt[rt<<1];
    }
    void update(int rt,int l,int r,int pos,double x){
        if(l==r){
            Max[rt]=x;
            cnt[rt]=1;
            return;
        }
        else if(pos<=mid)update(lson,pos,x);
        else update(rson,pos,x);
        Max[rt]=Max[rt<<1]>Max[(rt<<1)|1] ? Max[rt<<1] : Max[(rt<<1)|1];
        cnt[rt]=cnt[rt<<1]+cal(rson,Max[rt<<1]);
    }
}T;
int main(){
    //File();
    scanf("%d%d",&n,&m);
    REP(i,1,m){
        int x,y;
        scanf("%d%d",&x,&y);
        T.update(1,1,n,x,((double)y)/x);
        write(T.cnt[1],'\n');
    }
    return 0;
}