《挑战程序设计竞赛》3.2.5 常用技巧-坐标离散化 AOJ0531（1RE）

AOJ3061

http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=0531

题意

涂色：为了宣传信息竞赛，要在长方形的三合板上喷油漆来制作招牌。三合板上不需要涂色的部分预先贴好了护板。被护板隔开的区域要涂上不同的颜色，比如上图就应该涂上5种颜色。
请编写一个程序计算涂色数量，输入数据中，保证看板不会被护板全部遮住，并且护板的边一定是水平或垂直的。
输入：
第一个数是宽w（1 ≤ w ≤ 1000000），第二个数是高h（1 ≤ h ≤ 1000000）。
第二行是护板的数量n（1 ≤ n ≤ 1000），接着n行是每个护板的左下角坐标 (x1 , y1 )和右上角坐标 (x2 , y2 )，用空格隔开： x1 , y1 , x2 , y2 (0 ≤ x1< x2 ≤ w, 0 ≤ y1 < y2 ≤ h 都是整数)
招牌的坐标系如下，左下角是 (0, 0) ，右上角是(w, h) ， 测试集中的30%都满足w ≤ 100, h ≤ 100, n ≤ 100。

输出：
一个整数，代表涂色数量。
输入输出的例子
input.txt
15 6
10
1 4 5 6
2 1 4 5
1 0 5 1
6 1 7 5
7 5 9 6
7 0 9 2
9 1 10 5
11 0 14 1
12 1 13 5
11 5 14 6
0 0
output.txt
5

思路

坐标范围很大但是给定的坐标较少时，可以用坐标离散化处理。这里我按照书中例题写的程序最后报RE，不知道为什么，参考别人AC的程序检查了好几个小时也没检查出来。先贴出来以后再查吧。
另外能AC的程序参见别人的博客：AOJ 0531 Paint Color 题解 《挑战程序设计竞赛》

代码（RE，未通过）

#include <iostream>
#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;

const int N = 1000;

int w, h, n;
bool map[(N+1)*6][(N+1)*6];

int compose(int z[N][2], int t)
{
    vector<int> zs;
    for (int j = -1; j <= 1; j ++) {
        for (int i = 0; i < n; i ++) {
            for (int k = 0; k <= 1; k ++) {
                int zn = z[i][k]+j;
                if (0 <= zn && zn < t) zs.push_back(zn);
            }
        }
    }
    sort(zs.begin(), zs.end());
    zs.erase(unique(zs.begin(), zs.end()), zs.end());

    for (int i = 0; i < n; i ++) {
        for (int k = 0; k <= 1; k ++) {
            z[i][k] = find(zs.begin(), zs.end(), z[i][k]) - zs.begin();
        }
    }
    return zs.size();
}

void DFS(int x, int y)
{
    map[x][y] = false;
    int dir[4][2] = {{0, 1}, {1, 0}, {0, -1}, {-1, 0}};
    for (int k = 0; k < 4; k ++) {
        int nx = x + dir[k][0];
        int ny = y + dir[k][1];
        if (0 <= nx && nx < w && 0 <= ny && ny < h && map[nx][ny])
            DFS(nx, ny);
    }
}

int solve()
{
    int ans = 0;
    for (int i = 0; i < w; i ++) {
        for (int j = 0; j < h; j ++) {
            if (map[i][j]) {
                ans ++;
                //printf("%d, %d\n", i, j);
                DFS(i, j);
            }
        }
    }
    return ans;
}

int main(void)
{
    int x[N][2], y[N][2];

    while (cin >> w >> h) {
        if (!w && !h) break;

        cin >> n;
        for (int i = 0; i < n; i ++) {
            scanf("%d%d%d%d", &x[i][0], &y[i][0], &x[i][1], &y[i][1]);
            x[i][1]--, y[i][1]--;
        }

        w = compose(x, w);
        h = compose(y, h);

        for (int i = 0; i < w; i ++)
            fill(map[i], map[i]+h, true);
        for (int i = 0; i < n; i ++) {
            for (int j = x[i][0]; j <= x[i][1]; j ++) {
                for (int k = y[i][0]; k <= y[i][1]; k ++) {
                    map[j][k] = false;
                }
            }
        }

        printf("%d\n", solve());
    }

    return 0;
}