POJ2446——匈牙利算法讲解

POJ2446——chessboard（匈牙利算法） 二分图最大匹配数。

今天第一次写博客，给大家讲讲一个比较著名的二分图算法——匈牙利算法。

先来普及一下二分图的概念，

如图中有两个集合，U和V。只要U、V同一集合内点不直接相交（可以间接相连）就属于二分图。

匹配：

一条边的两个端点分别在两个集合里，这条边就是匹配边。最大匹配数图中最多能存在的没有公共端点的匹配边的数目

上图最大匹配数是4，U集合有点没连上也没办法，已经没有空的V点了。

增广路

假设只确定了一条匹配边（红色的），那么两条蓝色和红色的组成一条增广路，增广路的定义是：起点和终点是非匹配点，路径由匹配边与非匹配边交替。 显然，一条增广路的非匹配边比匹配边多1条，我们将蓝色的当作匹配边，红色当作非匹配边，这四个点的匹配边数就可以加1。 匈牙利算法的原理就是通过找增广路径来更新最大匹配边数。

具体匈牙利算法模拟可以看下面链接。

嗯，里面用男女关系来讲解确实简单易懂
趣味算法系列之匈牙利算法

下面讲题

Alice and Bob often play games on chessboard. One day, Alice draws a board with size M * N. She wants Bob to use a lot of cards with size 1 * 2 to cover the board. However, she thinks it too easy to bob, so she makes some holes on the board (as shown in the figure below).

We call a grid, which doesn’t contain a hole, a normal grid. Bob has to follow the rules below:
1. Any normal grid should be covered with exactly one card.
2. One card should cover exactly 2 normal adjacent grids.

Some examples are given in the figures below:

A VALID solution


An invalid solution, because there exists a grid, which is not covered.
Your task is to help Bob to decide whether or not the chessboard can be covered according to the rules above.
input
There are 3 integers in the first line: m, n, k (0 < m, n <= 32, 0 <= K < m * n), the number of rows, column and holes. In the next k lines, there is a pair of integers (x, y) in each line, which represents a hole in the y-th row, the x-th column.
output
If the board can be covered, output “YES”. Otherwise, output “NO”.
Sample input
4 3 2
2 1
3 3
Sample output
YES

题意是要用多条纸带覆盖一幅图的格子，因为一条纸带只能覆盖两个相邻的格子，我们将有洞的排除掉，将好的分成奇数和偶数号格子（我这里是根据横纵坐标和的奇偶性来分）。红色是奇数格子的编号，蓝色是偶数格子的编号。

把奇数格子和偶数格子当成二分图的两个集合求最大匹配边即可。如果最大匹配边数乘2和剩下格子数相等，就可以全覆盖。如果剩下格子为单数必定不能全覆盖。

#include<iostream>#include<cstdio>#include<string.h>using namespace std;const int M=35*35;bool used[M];     //记录偶数格子是否被使用；int dan[M];        //记录偶数格子和哪个单数格子配对；bool match[M][M];int map1[35][35];int v1,v2;      // v1单数格子数，v2偶数格子数;bool dfs(int x){    for(int i=1;i<v2;i++)    {        if(match[x][i]&&!used[i])        {            used[i]=true;            if(dan[i]==0||dfs(dan[i]))            {                dan[i]=x;                return 1;            }        }    }    return 0;}int seek(int n){    memset(dan,0,sizeof(dan));    int ans=0;    for(int i=1;i<n;i++)    {   memset(used,0,sizeof(used));        if(dfs(i))            ans+=1;    }    return ans;}int main(){    int n,m,k;    while(~scanf("%d%d%d",&m,&n,&k))    {        int x1,y1;        memset(map1,0,sizeof(map1));        for(int i=0;i<k;i++)        {            scanf("%d%d",&x1,&y1);            map1[y1][x1]=-1;           //题目的陷阱         }        if ((n*m - k) % 2 != 0)        {            printf("NO\n");            continue;        }        v1=v2=1;        for(int i=1;i<=m;i++)        {            for(int j=1;j<=n;j++)            {                if(map1[i][j]==0)                {                    if((i+j)%2==0)                    {                        map1[i][j]=v2++;                    }                    else                    {                        map1[i][j]=v1++;                    }                }            }        }        for(int i=1;i<=m;i++)        {            for(int j=1;j<=n;j++)            {                if((i+j)%2==0||map1[i][j]==-1)                {                    continue;                }                if(i-1>0&&map1[i-1][j]>0)                {                    match[map1[i][j]][map1[i-1][j]]=true;                }                if(j-1>0&&map1[i][j-1]>0)                {                    match[map1[i][j]][map1[i][j-1]]=true;                }                if(i+1<=m&&map1[i+1][j]>0)                {                    match[map1[i][j]][map1[i+1][j]]=true;                }                if(j+1<=n&&map1[i][j+1]>0)                {                    match[map1[i][j]][map1[i][j+1]]=true;                }            }        }        int cnt=seek(v1);        if(cnt==(n*m-k)/2)            printf("YES\n");        else            printf("NO\n");    }}

第一次编写博客可能写的不太好，这markdown编辑器还不太会用，而且换行的时候上文还经常自己打乱，心都烦。以后可能没周更一篇，欢迎大家浏览。。我还只是个acm新手，如果讲的不够清楚，大家多提意见哈。