hdu 1507 记录路径的二分匹配 **

题意：N*M的矩形，向其中填充1*2的小块矩形，黑色的部分不能填充，问最多可以填充多少块。
链接：点我
黑白棋最大匹配
将棋盘中i+j为奇数的做A集合，偶数的做B集合，相邻的则建立联系。于是便转换成寻找最大匹配的问题

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#include<map>
using namespace std;
#define MOD 1000000007
const int INF=0x3f3f3f3f;
const double eps=1e-5;
typedef long long ll;
#define cl(a) memset(a,0,sizeof(a))
#define ts printf("*****
");
const int MAXN=510;
int a[110][110];
int b[100];
int n,m,tt;
/* ***********************************************************
//二分图匹配（匈牙利算法的DFS实现）(邻接矩阵形式)
//初始化：g[][]两边顶点的划分情况
//建立g[i][j]表示i->j的有向边就可以了，是左边向右边的匹配
//g没有边相连则初始化为0
//uN是匹配左边的顶点数，vN是匹配右边的顶点数
//调用：res=hungary();输出最大匹配数
//优点：适用于稠密图，DFS找增广路，实现简洁易于理解
//时间复杂度:O(VE)
//************************************************************ */
//顶点编号从0开始的
int uN,vN;//u,v的数目，使用前面必须赋值
int g[MAXN][MAXN];//邻接矩阵
int linker[MAXN];
bool used[MAXN];
bool dfs(int u)
{
    for(int v = 0; v < vN;v++)
    if(g[u][v] && !used[v])
    {
        used[v] = true;
        if(linker[v] == -1 || dfs(linker[v]))
        {
            linker[v] = u;
            return true;
        }
    }
    return false;
}
int hungary()
{
    int res = 0;
    memset(linker,-1,sizeof(linker));
    for(int u = 0;u < uN;u++)
    {
        memset(used,false,sizeof(used));
        if(dfs(u))res++;
    }
    return res;
}
int main()
{
    int i,j,k;
    #ifndef ONLINE_JUDGE
    freopen("1.in","r",stdin);
    #endif
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        if(n==0&&m==0)  break;
        scanf("%d",&k);
        int u,v;
        cl(a);
        while(k--)
        {
            scanf("%d%d",&u,&v);
            u--,v--;
            a[u][v]=-1;
        }
        int tot=0;
        for(i=0;i<n;i++)
        {
            for(j=0;j<m;j++)
            {
                if(a[i][j]!=-1)
                {
                    b[tot]=i*m+j;
                    a[i][j]=tot++;
                }
            }
        }
        cl(g);
        for(i=0;i<n;i++)
        {
            for(j=0;j<m;j++)
            {
                if(a[i][j]!=-1&&((i+j)%2==1))
                {
                    int w=a[i][j];
                    if(i>0&&a[i-1][j]!=-1)
                        g[w][a[i-1][j]]=1;
                    if(i<n-1&&a[i+1][j]!=-1)
                        g[w][a[i+1][j]]=1;
                    if(j>0&&a[i][j-1]!=-1)
                        g[w][a[i][j-1]]=1;
                    if(j<m-1&&a[i][j+1]!=-1)
                        g[w][a[i][j+1]]=1;
                }
            }
        }
        vN=uN=tot;
        printf("%d
",hungary());
        for(i=0;i<tot;i++)
        {
            if(linker[i]!=-1)
            {
                int x1=b[i]/m;
                int y1=b[i]%m;
                int x2=b[linker[i]]/m;
                int y2=b[linker[i]]%m;
                printf("(%d,%d)--(%d,%d)
",x1+1,y1+1,x2+1,y2+1);
            }
        }
        printf("
");
    }
}