poj2446 Chessboard 二分匹配

这个做法跟poj3020很相似，但是因为是用1*2的矩阵去铺盖，看铺盖是否能完成，所以做法有些不一样，想法讲一下，由于表达能力有限。所以复制了一个讲法比较好的一段话，这段话讲的很清晰，有时候自己做得出来 但是要讲却还是讲不好的，下面给出这段话以及他的代码，我自己的代码会放到最后

转自：http://blog.csdn.net/hackbuteer1/article/details/7398008，建议去看看，他的博客很不错的

//poj_2446
/*==================================================*\
| 二分图匹配（匈牙利算法DFS 实现）
| INIT: g[][]邻接矩阵;
| 优点：实现简洁容易理解，适用于稠密图，DFS找增广路快。
| 找一条增广路的复杂度为O（E），最多找V条增广路，故时间复杂度为O（VE）
==================================================*/
#include<stdio.h>
#include<memory.h>

#define MAX 1089 //33*33
bool g[MAX][MAX]; //邻接矩阵，true代表有边相连
bool flag,visit[MAX];    //记录V2中的某个点是否被搜索过
int match[MAX];   //记录与V2中的点匹配的点的编号
int cnt;   //二分图中左边、右边集合中顶点的数目
bool hole[MAX][MAX];
int id[MAX][MAX];

// 匈牙利算法
bool dfs(int u)
{
	for (int i = 1; i <= cnt; ++i)
	{
		if (g[u][i] && !visit[i])   //如果节点i与u相邻并且未被查找过
		{
			visit[i] = true;   //标记i为已查找过
			if (match[i] == -1 || dfs(match[i]))   //如果i未在前一个匹配M中，或者i在匹配M中，但是从与i相邻的节点出发可以有增广路径
			{
				match[i] = u;  //记录查找成功记录，更新匹配M（即“取反”）
				return true;   //返回查找成功
			}
		}
	}
	return false;
}
int MaxMatch()
{
	int i,sum=0;
	memset(match,-1,sizeof(match));
	for(i = 1 ; i <= cnt ; ++i)
	{
		memset(visit,false,sizeof(visit));   //清空上次搜索时的标记
		if( dfs(i) )    //从节点i尝试扩展
		{
			sum++;
		}
	}
	return sum;
}

int main(void)
{
    int i,j,k,m,n,ans,y,x;
    while (scanf("%d %d %d",&m,&n,&k)!=EOF)
    {
		  memset(g,false,sizeof(g));
		  memset(hole,false,sizeof(hole));
          for (i = 1; i <= k; ++i)
		  {
			  scanf("%d %d",&y,&x);
              hole[x][y] = true;
		  }
		  if((m*n-k)&1)   //奇偶剪枝
		  {
			  puts("NO");
			  continue;
		  }
          cnt = 0;

          for (i = 1; i <= m; ++i)
          {
			  for (j = 1; j <= n; ++j)
			  {
				  if(hole[i][j] == false)   //对没有涂黑的点进行标号
				  {
					  id[i][j] = ++cnt;
				  }
			  }
          }
		  for (i = 1; i <= m; ++i)
          {
			  for (j = 1; j <= n; ++j)
			  {
				  if(hole[i][j] == false)
				  {
					  if(i-1>0 && hole[i-1][j] == false)   //建图。。要注意边界问题
						  g[ id[i][j] ][ id[i-1][j] ] = true;
					  if(i+1<=m && hole[i+1][j] == false)
						  g[ id[i][j] ][ id[i+1][j] ] = true;
					  if(j-1>0 && hole[i][j-1] == false)
						  g[ id[i][j] ][ id[i][j-1] ] = true;
					  if(j+1<=n && hole[i][j+1] == false)
						  g[ id[i][j] ][ id[i][j+1] ] = true;
				  }
			  }
		  }

		  ans = MaxMatch();
		  if (ans == cnt)
			  puts("YES");
		  else
			  puts("NO");
	}
    
    return 0;
}

接下来是我自己写的代码，做法跟poj3020很相似，但是建议先做本题再去做poj3020，意义不同

#include<iostream>
#include<cstdio>
#include<list>
#include<algorithm>
#include<cstring>
#include<string>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<cmath>
#include<memory.h>
#include<set>

#define ll long long
#define LL __int64
#define eps 1e-8
const ll INF=9999999999999;

using namespace std;

#define M 400000100

#define inf 0xfffffff

//vector<pair<int,int> > G;
//typedef pair<int,int> P;
//vector<pair<int,int>> ::iterator iter;
//
//map<ll,int>mp;
//map<ll,int>::iterator p;

vector<int>G[1212];

int mp[1212][1212];
int tempmp[1212][1212];
int marry[1212];
int has[1212][1212];
bool vis[1212];

int dis[2][4]={0,-1,0,1,1,0,-1,0};

int w,h,k;
int cnt;

void clear()
{
	cnt=0;
	memset(marry,-1,sizeof(marry));
	memset(tempmp,0,sizeof(tempmp));
	memset(mp,0,sizeof(mp));
	memset(has,0,sizeof(has));
	for(int i=0;i<1012;i++)
		G[i].clear();
}

bool dfs(int x)
{
	for(int i=1;i<=cnt;i++)
	{
		if(tempmp[x][i] && !vis[i])
		{
			vis[i]=true;
			if(marry[i]==-1 || dfs(marry[i]))
			{
				marry[i]=x;
				return 1;
			}
		}
	}
	return 0;
}

int main(void)
{
	while(cin>>w>>h>>k)
	{
		clear();
		int u,v;
		for(int i=0;i<k;i++)
		{
			scanf("%d %d",&v,&u);
			mp[u][v]=1;
		}
		if((w*h-k)&1)
		{
			puts("NO");
			continue;
		}
		for(int i=1;i<=w;i++)
		{
			for(int j=1;j<=h;j++)
			{
				if(mp[i][j])
					continue;
				has[i][j]=++cnt;
			}
		}
		for(int i=1;i<=w;i++)
		{
			for(int j=1;j<=h;j++)
			{
				if(has[i][j])
				{
					for(int k=0;k<4;k++)
					{
						int dx=dis[0][k]+i;
						int dy=dis[1][k]+j;
						if(has[dx][dy])
							tempmp[has[i][j]][has[dx][dy]]=1;
					}
				}
			}
		}
		int ans=0;
		for(int i=1;i<=cnt;i++)
		{
			memset(vis,false,sizeof(vis));
			if(dfs(i))
				ans++;
		}
		if(ans==cnt)
			puts("YES");
		else
			puts("NO");
	}
}