[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 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
　　

---

咱是来学线段树黑科技的……
把这种奇特的技巧出到NOIP模拟赛里的人究竟是什么心态……
讲题人问大家为什么没有想到这个这么经典的方法，全场寂静……
因为貌似全场80多个各省dalao根本™没有几个人写过（本机房个数为0）……

(╯‵□′)╯︵┻━┻

正解是在转换模型后在这个东西的基础上维护一个DP最小值……
出题人预告只会一天比一天难，然而还有两天……
又要被各省dalao按在地上摩擦，绝望.jpg（毕竟集训队dalao）

---

思路:
可以发现，一个楼房能被看到当且仅当主角看这栋楼房时的视线斜率大于前面的任何一个楼房。
那么这就是在维护极长上升子序列的最大长度……
所谓极长，就是能选的都必须选的意思……可以理解为局部最优解吧……

考虑线段树，记录只考虑当前区间时的极长上升序列的长度和区间最大值。
那么这么做最大的难点在于合并子树信息。
考虑以下几种情况:

1.左区间最大值大于右区间最大值
那么由于这是极长的，左区间的最大值会被选入，从而导致右区间毫无贡献，直接获得左区间的值即可。
相当于是右边的楼房全都被左边最高的楼房挡住了，所以右边毫无贡献。

2.左区间最大值小于右区间最大值
首先由于这是极长所以一定要选完左区间，直接累加。
其次就是处理右区间了。
由于右区间咱并不知道具体的贡献情况，那么递归处理。

令右区间的儿子分别为左子区间、右子区间，考虑如下情况:

2.1左子区间的最大值小于等于左区间的最大值
那么右区间的最大值出现在右子区间，且左子区间完全被左区间挡住，直接递归右子区间查找有效个数。

2.2左子区间的最大值大于左区间的最大值
那么左子区间此时有效，递归左子区间查找有效个数。
而右子区间的贡献则为右区间原本的答案减去左子区间原本的答案。
因为右区间的原本答案中属于右子区间部分的内容已经考虑过了左子区间的约束，直接减去属于左子区间的部分即可得到右子区间满足约束下的答案。

那么这就结束了……

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<algorithm>

using namespace std;

const int N=100009;
typedef double db;

inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0' || '9'<ch)ch=getchar();
    while('0'<=ch && ch<='9')x=x*10+(ch^48),ch=getchar();
    return x;
}

int n,m,x,y;
int t[N<<2];
db tm[N<<2];

inline int calc(int x,int l,int r,db maxh)
{
    if(l==r)
        return tm[x]>maxh;
    int mid=l+r>>1;
    if(tm[x<<1]>maxh)
        return calc(x<<1,l,mid,maxh)+t[x]-t[x<<1];
    return calc(x<<1|1,mid+1,r,maxh);
}

inline void update(int x,int l,int r)
{
    int mid=l+r>>1;
    tm[x]=max(tm[x<<1],tm[x<<1|1]);
    t[x]=t[x<<1]+calc(x<<1|1,mid+1,r,tm[x<<1]);
}

inline void modify(int x,int l,int r,int p,db val)
{
    if(l==r)
    {
        t[x]=1,tm[x]=val;
        return;
    }
    int mid=l+r>>1;
    if(p<=mid)
        modify(x<<1,l,mid,p,val);
    else
        modify(x<<1|1,mid+1,r,p,val);
    update(x,l,r);
}

int main()
{
    n=read();
    m=read();
    while(m--)
    {
        x=read();
        y=read();
        db h=(db)y/(db)x;
        modify(1,1,n,x,h);
        printf("%d\n",t[1]);
    }

    return 0;
}