[HDOJ5925]Coconuts（BFS，离散化，计数）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5925

题意：一个 R×C 的棋盘，有 n≤200 个格子是黑的，其他都是白的，问所有白色格子构成的四联通块有多大。

题解：离散化后 BFS。

这个离散化以后会有一个问题，因为这个题要统计每一个连通块里的点数，所以我在离散化的时候顺便对每一对相邻的障碍物之间有多少个空白做了计数，这样在bfs的时候只需要加上乘积就行了。

#include <bits/stdc++.h>
using namespace std;

typedef long long LL;
const int maxn = 100100;
const int maxm = 3048;
const int dx[5] = {0,0,1,-1};
const int dy[5] = {1,-1,0,0};
typedef struct Point {
    int x, y;
    Point(){}
    Point(int x,int y):x(x),y(y){}
}Point;
int vis[maxm][maxm];
Point p[maxm];
int K, N, M, n, m;
int tx[maxn], ty[maxn];
LL vx[maxn], vy[maxn];
vector<LL> ret;
queue<Point> q;
int d[maxn];

inline bool ok(int x, int y) {
    return x >= 1 && x <= n && y >= 1 && y <= m;
}

int gao(int* x, LL* vx, int& w, bool flag) {
    sort(x, x+w);
    w = unique(x, x+w) - x;
    int t = 1;
    d[0] = x[0];
    for(int i = 1; i < w; i++) {
        if(x[i] != x[i-1] + 1) d[t++] = x[i-1] + 1;
        d[t++] = x[i];
    }
    sort(d, d+t);
    int mx = t - 1;
    for(int i = 1; i < t; i++) {
        vx[i] = d[i] - d[i-1];
    }
    for(int i = 0; i < K; i++) {
        if(flag) p[i].x = lower_bound(d, d+t, p[i].x) - d;
        else p[i].y = lower_bound(d, d+t, p[i].y) - d;
    }
    return mx;
}

LL bfs(int x, int y) {
    while(!q.empty()) q.pop();
    q.push(Point(x, y)); vis[x][y] = 1;
    LL cnt = vx[x] * vy[y];
    while(!q.empty()) {
        Point t = q.front(); q.pop();
        for(int i = 0; i < 4; i++) {
            int tx = t.x + dx[i];
            int ty = t.y + dy[i];
            if(!ok(tx, ty)||vis[tx][ty]) continue;
            vis[tx][ty] = 1; cnt += vx[tx] * vy[ty];
            q.push(Point(tx, ty));
        }
    }
    return cnt;
}

int main() {
    // freopen("in", "r", stdin);
    int T, _ = 1;
    scanf("%d", &T);
    while(T--) {
        ret.clear(); N = 0, M = 0;
        memset(vis, 0, sizeof(vis));
        scanf("%d%d%d",&n,&m,&K);
        tx[N++] = 0; tx[N++] = n;
        ty[M++] = 0; ty[M++] = m;
        for(int i = 0; i < K; i++) {
            scanf("%d%d",&p[i].x,&p[i].y);
            tx[N++] = p[i].x; tx[N++] = p[i].x - 1;
            ty[M++] = p[i].y; ty[M++] = p[i].y - 1;
        }
        n = gao(tx, vx, N, 1); m = gao(ty, vy, M, 0);
        for(int i = 0; i < K; i++) {
            // BUG
            // if(p[i].x >= maxm || p[i].y >= maxm) while(1);
            vis[p[i].x][p[i].y] = 1;
        }
        for(int i = 1; i <= n; i++) {
            for(int j = 1; j <= m; j++) {
                if(!vis[i][j]) ret.push_back(bfs(i, j));
            }
        }
        sort(ret.begin(), ret.end());
        printf("Case #%d:
", _++);
        printf("%d
", ret.size());
        for(int i = 0; i < ret.size(); i++) {
            if(i) printf(" ");
            printf("%I64d", ret[i]);
        }
        printf("
");
    }
    return 0;
}