BZOJ 2957：楼房重建（分块）

http://www.lydsy.com/JudgeOnline/problem.php?id=2957

题意：……

思路：对于每一个块，维护一个单调递增的斜率（因为小于前面的斜率的话是肯定看不见的），每次修改都暴力重建这个块，然后求和的时候可以二分查找当前的块比之前的块的最大斜率大的元素个数（注意精度调大一点），然后就可以得到答案了。（这么暴力感觉好爽啊，upper_bound（）好方便啊）

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
using namespace std;
#define N 100010
#define INF 0x3f3f3f3f
typedef long long LL;
const double eps = 1e-10;
int n, m, sz, num, l[N], r[N], belong[N], cnt[N];
double a[N], block[1010][1010];

void Build() {
    memset(a, 0, sizeof(a));
    memset(cnt, 0, sizeof(cnt));
    memset(block, 0, sizeof(block));
    sz = sqrt(n);
    num = n / sz; if(n % sz) num++;
    for(int i = 1; i <= num; i++)
        l[i] = (i - 1) * sz + 1, r[i] = i * sz;
    r[num] = n;
    for(int i = 1; i <= n; i++)
        belong[i] = (i - 1) / sz + 1;
}

void Update(int x) { // 暴力重建块
    double tmp = 0.0; cnt[x] = 0;
    for(int i = l[x]; i <= r[x]; i++) // 维护单调递增序列
        if(tmp < a[i]) tmp = a[i], block[x][++cnt[x]] = a[i];
}

int Query() {
    int ans = 0; double tmp = 0.0;
    for(int i = 1; i <= num; i++) { // 二分查找,注意精度
        if(!cnt[i]) continue;
        int now = cnt[i] - (upper_bound(block[i] + 1, block[i] + 1 + cnt[i], tmp + eps) - block[i] - 1);
        if(now) tmp = block[i][cnt[i]];
        ans += now;
    }
    return ans;
}

int main() {
    while(~scanf("%d%d", &n, &m)) {
        Build();
        for(int i = 1; i <= m; i++) {
            int id, w;
            scanf("%d%d", &id, &w);
            a[id] = (double)w / id;
            Update(belong[id]);
            printf("%d
", Query());
        }
    }
    return 0;
}