这个做法跟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");
}
}