XTU 二分图和网络流 练习题 B. Uncle Tom's Inherited Land*

B. Uncle Tom's Inherited Land*

Time Limit: 1000ms

Memory Limit: 32768KB

64-bit integer IO format: %I64d      Java class name: Main

Special Judge

 

Your old uncle Tom inherited a piece of land from his great-great-uncle. Originally, the property had been in the shape of a rectangle. A long time ago, however, his great-great-uncle decided to divide the land into a grid of small squares. He turned some of the squares into ponds, for he loved to hunt ducks and wanted to attract them to his property. (You cannot be sure, for you have not been to the place, but he may have made so many ponds that the land may now consist of several disconnected islands.)

Your uncle Tom wants to sell the inherited land, but local rules now regulate property sales. Your uncle has been informed that, at his great-great-uncle's request, a law has been passed which establishes that property can only be sold in rectangular lots the size of two squares of your uncle's property. Furthermore, ponds are not salable property.

Your uncle asked your help to determine the largest number of properties he could sell (the remaining squares will become recreational parks). 

 

Input

Input will include several test cases. The first line of a test case contains two integers N and M, representing, respectively, the number of rows and columns of the land (1 <= N, M <= 100). The second line will contain an integer K indicating the number of squares that have been turned into ponds ( (N x M) - K <= 50). Each of the next K lines contains two integers X and Y describing the position of a square which was turned into a pond (1 <= X <= N and 1 <= Y <= M). The end of input is indicated by N = M = 0.

 

Output

For each test case in the input your program should first output one line, containing an integer p representing the maximum number of properties which can be sold. The next p lines specify each pair of squares which can be sold simultaneity. If there are more than one solution, anyone is acceptable. there is a blank line after each test case. See sample below for clarification of the output format.

 

Sample Input

4 4
6
1 1
1 4
2 2
4 1
4 2
4 4
4 3
4
4 2
3 2
2 2
3 1
0 0

Sample Output

4
(1,2)--(1,3)
(2,1)--(3,1)
(2,3)--(3,3)
(2,4)--(3,4)

3
(1,1)--(2,1)
(1,2)--(1,3)
(2,3)--(3,3)

解题：根据奇偶性建图+最大匹配数+最大匹配数的记录。根据奇偶性建图：在一个平面坐标系中，某个格子的 纵横坐标和 与其相领格子的 纵横坐标和 的绝对值相差1，即相差一个数。由于奇数与偶数相间分布，故只选择纵横坐标和为奇数或者纵横坐标和为偶数的点进行建图（从始至终只能选择其中一种方式），这样保证了是一个二分图。

这题的写法，把原先的匈牙利算法的线段的端点的由的写成了二维的，算法还是那样的。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <climits>
#include <vector>
#include <queue>
#include <cstdlib>
#include <string>
#include <set>
#define LL long long
#define INF 0x3f3f3f3f
using namespace std;
const int maxn = 110;
struct{
    int x,y;
}link[maxn][maxn];
int n,m;
bool mp[maxn][maxn],vis[maxn][maxn];
bool isIn(int x,int y){
    if(x < 1 || x > n || y < 1 || y > m)
        return false;
    return true;
}
bool dfs(int x,int y){
    static const int dir[4][2] = {{1,0},{-1,0},{0,1},{0,-1}};
    for(int i = 0; i < 4; i++){
        int tx = x+dir[i][0];
        int ty = y+dir[i][1];
        if(isIn(tx,ty) && !vis[tx][ty] && !mp[tx][ty]){
            vis[tx][ty] = true;
            if(link[tx][ty].x == -1 || dfs(link[tx][ty].x,link[tx][ty].y)){
                link[tx][ty].x = x;
                link[tx][ty].y = y;
                return true;
            }
        }
    }
    return false;
}
int main(){
    int i,j,u,v,k;
    while(scanf("%d%d",&n,&m),n+m){
        scanf("%d",&k);
        memset(mp,false,sizeof(mp));
        while(k--){
            scanf("%d%d",&u,&v);
            mp[u][v] = 1;
        }
        int ans = 0;
        memset(link,-1,sizeof(link));
        for(i = 1; i <= n; i++){
            for(j = 1; j <= m; j++){
                memset(vis,false,sizeof(vis));
                if(((i+j)&1) && !mp[i][j] && dfs(i,j)) ans++;
            }
        }
        printf("%d\n",ans);
        for(i = 1; i <= n; i++){
            for(j = 1; j <= m; j++){
                if(link[i][j].x != -1)
                    printf("(%d,%d)--(%d,%d)\n",link[i][j].x,link[i][j].y,i,j);
            }
        }
        printf("\n");
    }
    return 0;
}