poj2446 （二分图匹配模型）

题意：给一个m*n的矩形棋牌，一些格子有洞。用1*2的卡片来覆盖剩余的格子，有洞的格子不能被覆盖，每个格子只能被覆盖一次。问是否能用这种卡片将所有剩下的格子覆盖完。

                                             

解法：将格子分成两种，一种是（i+j）是奇数的，一种是偶数的。会发现匹配时总是奇数格子和偶数格子匹配。将两种格子标号并建图，由于是稀疏图，所以用了邻接表。匹配时仍然是匈牙利。如果最大匹配数*2=剩下的格子数，则YES，否则 NO;

代码：

/****************************************************
* author:xiefubao
*******************************************************/
#pragma comment(linker, "/STACK:102400000,102400000")
#include <iostream>
#include <cstring>
#include <cstdlib>
#include <cstdio>
#include <queue>
#include <vector>
#include <algorithm>
#include <cmath>
#include <map>
#include <set>
#include <stack>
#include <string.h>

using namespace std;

#define eps 1e-8
typedef long long LL;
int m,n,k;
int rem[40][40];
struct edge
{
    int v;
    int next;
} edges[30000];
int head[1000];
int match[1000];
int used[1000];
int ncase=0;
int count1=0;
void addedge(int u,int v)
{
    edges[count1].v=v;
    edges[count1].next=head[u];
    head[u]=count1++;
}


void make(int i,int j)
{
    if(i>1&&rem[i-1][j]!=-1)
    {
        if((i+j)&1)
        addedge(rem[i][j],rem[i-1][j]);
        else
        addedge(rem[i-1][j],rem[i][j]);
    }
    if(j>1&&rem[i][j-1]!=-1)
    {
        if((i+j)&1)
       addedge(rem[i][j],rem[i][j-1]);
       else
       addedge(rem[i][j-1],rem[i][j]);
    }
}
bool dfs(int t)
{
    for(int i=head[t];i!=-1;i=edges[i].next)
    {
        int a=edges[i].v;
        if(used[a]!=ncase)
        {
            used[a]=ncase;
            if(match[a]==-1||dfs(match[a]))
            {
                match[a]=t;
                return true;
            }
        }
    }
    return false;
}
int main()
{
  while(scanf("%d%d%d",&m,&n,&k)!=EOF)
  {
      memset(rem,0,sizeof rem);
  for(int i=0;i<k;i++)
  {
      int a,b;
      scanf("%d%d",&a,&b);
      rem[b][a]=-1;
  }
  if((m*n-k)&1){cout<<"NO\n";continue;}
  memset(head,-1,sizeof head);
  memset(match,-1,sizeof match);
  count1=0;
  int tool1=0,tool2=0;
  for(int i=1;i<=m;i++)
    for(int j=1;j<=n;j++)
  {
      if(rem[i][j]!=-1)
      {
          if((i+j)&1)
            rem[i][j]=tool1++;
          else
            rem[i][j]=tool2++;
          make(i,j);
      }
  }
  int ans=0;
  for(int i=0;i<tool1;i++)
  {
      ncase++;
      if(dfs(i))
        {
            ans++;
        }
  }
  if(n*m-k==ans*2)
    cout<<"YES\n";
  else
    cout<<"NO\n";
  }
   return 0;
}