BZOJ2957: 楼房重建

楼房i被看见的条件是他和起点连线的斜率是1~i里最大的

可以用线段树做
维护每个区间的ans，max
合并区间的时候，
左区间的ans是可以直接加的，
对于右区间，当左区间的max>右区间max时，整个区间被挡住，没有贡献
当左区间的max<右区间max时，需要递归下右区间找
递归过程中，带一个这个区间前方的斜率max，
当max<左子区间的max时，这个带下来的max对右子区间的贡献无影响，通过整个区间的ans-左子区间的ans算粗右子区间的答案直接加上就行，对于左子区间再递归下去
当max>左子区间的max时，左子区间无贡献，递归下右子区间
这个做法..我总觉得是$log^2$的

然后….说了这么多…其实我没打线段树…..（憋打我）

我写了分块
分块的话，每个块内维护一个斜率单调增的子序列，修改一个值时，这个块内暴力重构，因为修改一个值只会对右方的贡献产生影响，向右顺着块扫一遍，记录之前的最大值，在每个块内二分找有多少个值比他大
code：

#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const int maxn = 110000;
const int maxk = 110;
const double eps=1e-10;
int id[maxn],L[maxk];
int len[maxk],N;
double a[maxn],q[maxk][1200];
int find_(int k,double x)
{
    int l=1,r=len[k];
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(q[k][mid]-x>eps) r=mid-1;
        else l=mid+1;
    }
    return r+1;
}
void rebuild_(int k)
{
    int llen=0;
    for(int i=L[k];i<L[k+1];i++)if(a[i]-q[k][llen]>eps) 
            q[k][++llen]=a[i];
    len[k]=llen;
}
int n,m;

int main()
{
    scanf("%d%d",&n,&m); N=(3.0*sqrt(n));
    for(int i=1;i<=n;i++) id[i]=(i-1)/N+1;
    for(int i=1;i<=id[n];i++) L[i]=(i-1)*N+1;
    L[id[n]+1]=n+1;
    while(m--)
    {
        int x,y; scanf("%d%d",&x,&y);
        a[x]=(double)y/x;
        rebuild_(id[x]);
        double las=0.0; int ret=0;
        for(int i=1;i<=id[n];i++)if(len[i])
        {
            int k=find_(i,las);
            if(k<=len[i]) ret+=len[i]-k+1,las=q[i][len[i]];
        }
        printf("%d\n",ret);
    }

    return 0;
}