BZOJ2957: 楼房重建

n<=100000栋楼房，第i栋在位置i，一开始全没建成高度0，m<=1e5个操作，每次操作把一栋楼x的高度修改为y，并问这次修改后从0位置往右边看能看到多少栋楼。他能看到一栋楼，必须这栋楼有高度，并且设楼高度Hi，那么(0,0)到(i,Hi)这条线段上没有其他的楼。1<=y<=1e9

说半天就是维护个斜率啦。。单点修改容易，不过这个信息有点难维护。。

可以用一个分治方法维护，现在是要求：从第一个有高度的楼房开始，斜率递增的最长序列有多长。把这个区间分成左右区间：

假如左右区间的要求的东西（即从该区间的第一个有高度的楼房开始的最长递增序列有多长）已知，那么在合并时，左区间的答案是可以直接算的，右区间收到左区间最大值的限制。也就是说，现在需要在右区间找：从该区间第一个大于左区间最大值的楼房开始的最长递增序列有多长。注意这里的递增都是在说斜率。

在右区间里面，关注一下左右子区间的信息和左区间最大值的关系，如果左子区间的最大值比左区间最大值小，那么接下来的在右子区间找大于左区间最大值的最长递增序列即可；反之，如果左子区间的最大值比左区间的大，那么右子区间在右区间里面表现出来的答案是不会变的，也就是：右区间的自己的以第一个有高度的楼房为起点的最长序列长度，减去左子区间的以……，剩下的部分就是右子区间收到左子区间最大值限制，以第一个大于左子区间最大值的楼房开始的最长序列，注意到这里说的还是“反之”的内容，这右子区间已知，那剩下的就是左子区间里面找大于左区间最大值的最长递增序列长度。

左右区间的操作？交给线段树或平衡树

#include<stdio.h>
#include<string.h>
#include<stdlib.h>
//#include<assert.h>
#include<math.h>
#include<algorithm>
//#include<iostream>
using namespace std;

int n,m;
#define maxn 200011
const double eps=1e-12;
struct SMT
{
    struct Node
    {
        double Max;int ans;
        int l,r;
        int ls,rs;
    }a[maxn<<1];
    int size;
    SMT() {size=0;}
    void build(int &x,int L,int R)
    {
        x=++size;
        a[x].l=L;a[x].r=R;
        a[x].Max=0;a[x].ans=0;
        if (L==R) {a[x].ls=a[x].rs=0;return;}
        const int mid=(L+R)>>1;
        build(a[x].ls,L,mid);
        build(a[x].rs,mid+1,R);
    }
    void build() {int x;build(x,1,n);}
    int findans(int x,double v)
    {
//        cout<<x<<endl;
        if (a[x].Max-v<eps) return 0;
        if (a[x].l==a[x].r) return 1;
        const int &p=a[x].ls,&q=a[x].rs;
        if (v-a[p].Max>-eps) return findans(q,v);
        return findans(p,v)+a[x].ans-a[p].ans;
    }
    void up(int x)
    {
        const int &p=a[x].ls,&q=a[x].rs;
        a[x].Max=max(a[p].Max,a[q].Max);
        a[x].ans=a[p].ans+findans(q,a[p].Max);
    }
    int q;double v;
    void besingle(int x,double v)
    {
        a[x].Max=v;
        a[x].ans=1;
    }
    void be(int x)
    {
//        cout<<x<<endl;
        if (a[x].l==a[x].r && a[x].l==q) besingle(x,v);
        else
        {
            const int mid=(a[x].l+a[x].r)>>1;
            if (q<=mid) be(a[x].ls);
            else be(a[x].rs);
            up(x);
        }
    }
    void be(int x,double v)
    {
        q=x;this->v=v;
        be(1);
    }
}t;
int main()
{
    scanf("%d%d",&n,&m);
    t.build();//cout<<":)";
    int x,y;
    for (int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
//        assert(x<=n);
        t.be(x,1.0*y/x);
        printf("%d
",t.a[1].ans);
    }
    return 0;
}