HDU 1507 Uncle Tom's Inherited Land*(行列匹配||棋盘覆盖||黑白染色||奇偶匹配(1X2的矩形覆盖))

题意：n*m的棋盘里，有几个点是不可以放的，要求在剩下的可以放的点上放置1X2的矩形，问最多可以放多少个
思路：观察发现，如果把矩阵按坐标和的奇偶分为A,B两个集合，每一个点可以到达周围的四个点，正好是1X2的覆盖，也就是分为A,B集合后的最大匹配。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<set>
#include<map>
#include<string>
#include<cstring>
#include<stack>
#include<queue>
#include<vector>
#include<limits>
#include<ctime>
#include<cassert>
#include<cstdlib>
#define lson (rt<<1),L,M
#define rson (rt<<1|1),M+1,R
#define M ((L+R)>>1)
#define cl(a,b) memset(a,b,sizeof(a));
#define LL long long
#define P pair<int,int>
#define X first
#define Y second
#define pb push_back
#define fread(a) freopen(a,"r",stdin);
#define fwrite(a) freopen(a,"w",stdout);
using namespace std;
const int maxn=105;
const int inf=999999;

int a[maxn][maxn];
vector<int> G[maxn*maxn];
int matching[maxn*maxn];
void build_map(int n,int m){
    for(int i=0;i<n;i++){
        for(int j=0;j<m;j++)if(!a[i][j]){
            if((i+j)%2==0){//坐标和为偶数的作为二分图的左边，这个自己定
                if(i-1>=0&&!a[i-1][j]){//四个可达的点，作为二分图右边，连边
                    G[i*m+j].pb((i-1)*m+j);
                }
                if(i+1<n&&!a[i+1][j]){
                    G[i*m+j].pb((i+1)*m+j);
                }
                if(j-1>=0&&!a[i][j-1]){
                    G[i*m+j].pb(i*m+j-1);
                }
                if(j+1<m&&!a[i][j+1]){
                    G[i*m+j].pb(i*m+j+1);
                }
            }
        }
    }
}

bool vis[maxn*maxn];
int num;
bool dfs(int u){
    int N=G[u].size();
    for(int i=0;i<N;i++){
        int v=G[u][i];
        if(vis[v])continue;
        vis[v]=true;
        if(matching[v]==-1||dfs(matching[v])){
            matching[v]=u;
            return true;
        }
    }
    return false;
}

int hungar(){
    int ans=0;
    cl(matching,-1);
    for(int i=0;i<num;i++){
        cl(vis,false);
        if(dfs(i))ans++;
    }
    return ans;
}
int main(){
    int n,m,k;
    while(~scanf("%d%d",&n,&m)){
        if(n==0&&m==0)return 0;
        scanf("%d",&k);
        cl(a,0);
        for(int i=0;i<maxn*maxn;i++)G[i].clear();
        for(int i=0;i<k;i++){
            int x,y;
            scanf("%d%d",&x,&y);
            x--;y--;
            a[x][y]=1;
        }

        build_map(n,m);
        num=n*m;
        int ans=hungar();
        printf("%d\n",ans);
        for(int i=0;i<num;i++)if(matching[i]!=-1){
            printf("(%d,%d)--(%d,%d)\n",i/m+1,i%m+1,matching[i]/m+1,matching[i]%m+1);
        }
        printf("\n");
    }
    return 0;
}