【分块】bzoj2957 楼房重建

http://www.cnblogs.com/wmrv587/p/3843681.html

ORZ 分块大爷。思路很神奇也很清晰。

把 块内最值 和 块内有序 两种良好的性质结合起来，非常棒地解决了这个问题。

图中黑色的楼房即为每个块内的“可视序列”，显而易见，在块内它们的K（斜率）是单增的。

由于上图中第一个块的maxK比后面两个块的maxK都要大，所以后两个块对答案没有贡献，这也是显然的。这就是维护maxK的意义所在。

否则，若一个块可以更新maxK的话，则其中的部分楼房是“可见的”，具体来说，就是在那个比之前的maxK要大的楼房的后面的在可视序列中的楼房数。<---请从方链接看原版题解。

另外，并不会像他说的，基本不会T，如果把块大小开得合适。

#include<cstdio>
#include<algorithm>
#include<vector>
#include<cmath>
#include<cstring>
using namespace std;
inline double max(const double &a,const double &b){return a>b?a:b;}
vector<double>See[320];
double k[100001];
int sz,sum,l[320],r[320],num[100001],n,m,x,y;
double maxv[320];
int Res,Num;char C,CH[12];
inline int G()
{
    Res=0;C='*'; 
    while(C<'0'||C>'9')C=getchar();
    while(C>='0'&&C<='9'){Res=Res*10+(C-'0');C=getchar();}
    return Res;
}
inline void P(int x)
{
    Num=0;if(!x){putchar('0');puts("");return;}
    while(x>0)CH[++Num]=x%10,x/=10;
    while(Num)putchar(CH[Num--]+48);
    puts("");
}
void makeblock()
{
    memset(maxv,0,sizeof(maxv));
    sz=sqrt((double)n*1.05);
    for(sum=1;sum*sz<n;sum++)
      {
        l[sum]=(sum-1)*sz+1;
        r[sum]=sum*sz;
        for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
      }
    l[sum]=sz*(sum-1)+1;
    r[sum]=n;
    for(int i=l[sum];i<=r[sum];i++) num[i]=sum;
}
inline void update()
{
    k[x]=(double)y/x;
    See[num[x]].clear();
    maxv[num[x]]=0.0;
    for(int i=l[num[x]];i<=r[num[x]];i++)
      if(k[i]>maxv[num[x]])
        {
          maxv[num[x]]=k[i];
          See[num[x]].push_back(k[i]);
        }
}
inline void query()
{
    int ans=0;double tmp=0.0;
    for(int i=1;i<=sum;i++)
      if(!See[i].empty())
        {
            ans+=See[i].end()-upper_bound(See[i].begin(),See[i].end(),tmp);
            tmp=max(tmp,maxv[i]);
        }
    P(ans);
}
int main()
{
    n=G();m=G();makeblock();
    for(int i=1;i<=m;i++){x=G();y=G();update();query();}
    return 0;
}