[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

Sample Output

数据约定
　　对于所有的数据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=;RG int c=Blue(),f=;
    for(; c<||c>; c=Blue())
        if(c==)f=-;
    for(; c>&&c<; c=Blue())
        x=(x<<)+(x<<)+c-;
    x*=f; }

double tr[MaxN<<];

int n,m,ans,sum[MaxN<<];

int calc(RG int x,RG int y,RG int k,RG double val) {
    if(x == y)return tr[k] > val;
    if(tr[k<<] < val)
        return calc(mid+,y,k<<|,val);
    return sum[k]-sum[k<<]+calc(x,mid,k<<,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]=;
        return; }
    pos <= mid ? modify(x,mid,k<<,pos,val) : modify(mid+,y,k<<|,pos,val);
    tr[k]=dmax(tr[k<<],tr[k<<|]);
    sum[k]=sum[k<<]+calc(mid+,y,k<<|,tr[k<<]);
}

int main() {
    Rin(n),Rin(m);
    while(m--) {
        RG int pos,x;
        Rin(pos),Rin(x);
        modify(,n,,pos,(double)x/pos);
        printf("%d\n",sum[]); }
    fclose(stdin);
    return ; }