bzoj 2957 楼房重建 分块

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

Data Constraint

Hint
对于所有的数据1<=Xi<=N，1<=Yi<=10^9

测试点 N,M 测试点 N,M

1 <=100 6 <=50000

2 <=5000 7 <=70000

3 <=50000 8 <=80000

4 <=100000 9 <=90000

5 <=30000 10 <=100000

分析：
一个点被看到只有他的斜率大于前面的任何一个，等于也不行，样例可得。单点修改可以想到分块。每个块记录，假设站在这个块的最前面，这个块上的可以看到的房子的序列。每次修该暴力改变该块上的序列。假设第一个块序列的最大值为p1，那么第二个块的答案就是第二个序列的第一个大于p1的数到第二个块的最大值p2……，可以用二分求出第一个大于p1的数。

(强制类型转换(double)y/x，打成double(y/x)不会爆编译，以为是精度问题，调了一个小时)

代码：

/**************************************************************
    Problem: 2957
    User: beginend
    Language: C++
    Result: Accepted
    Time:4140 ms
    Memory:10328 kb
****************************************************************/

#include <iostream>
#include <cstdio>
#include <cmath>
#define N 100007
#define maxn 1002 
#define fo(i,a,b) for (int i=a;i<=b;i++)

using namespace std;

int block,num,belong[N],l[maxn],r[maxn],sum[maxn];
int n,m;
int x,y;
double T[maxn][maxn],a[N];

void build()
{
    block=trunc(sqrt(n))+1;
    num=n/block;
    if (n%block) num++;
    fo(i,1,num)
    {
        l[i]=(i-1)*block+1;
        r[i]=i*block;
    }
    r[num]=n;
    fo(i,1,n)
    {
        belong[i]=(i-1)/block+1;
    }
}

void change(int op,double x)
{
    a[op]=x;
    double p=0;
    int c=belong[op];
    sum[c]=0;
    fo(j,l[c],r[c])
    {
        if (a[j]>p)
        {
            p=a[j];
            T[c][++sum[c]]=p;
        }
    }
}

int ask()
{
    int l,r,k=0;
    int ans=0;
    double p=0; 
    fo(j,1,num)
    {
        l=1; r=sum[j];
        if ((r==0) || (T[j][r]<=p)) continue;
        while (l<=r)
        {
            int mid=(l+r)/2;
            if (T[j][mid]>p)
            {
                k=mid;
                r=mid-1;
            }
            else l=mid+1;
        }
        ans+=sum[j]-k+1;
        p=T[j][sum[j]]; 
    }
    return ans;
}

int main()
{
    scanf("%d%d",&n,&m);
    build();        
    fo(i,1,m)
    {       

        scanf("%d%d",&x,&y);                
        change(x,(double)y/x);
        printf("%d\n",ask());
    }
}