hdu 1507 记录路径的二分匹配

  1 #include<cstdio>
  2 #include<iostream>
  3 #include<algorithm>
  4 #include<cstring>
  5 #include<cmath>
  6 #include<queue>
  7 #include<map>
  8 using namespace std;
  9 #define MOD 1000000007
 10 const int INF=0x3f3f3f3f;
 11 const double eps=1e-5;
 12 typedef long long ll;
 13 #define cl(a) memset(a,0,sizeof(a))
 14 #define ts printf("*****\n");
 15 const int MAXN=510;
 16 int a[110][110];
 17 int b[100];
 18 int n,m,tt;
 19 /* ***********************************************************
 20 //二分图匹配（匈牙利算法的DFS实现）(邻接矩阵形式)
 21 //初始化：g[][]两边顶点的划分情况
 22 //建立g[i][j]表示i->j的有向边就可以了，是左边向右边的匹配
 23 //g没有边相连则初始化为0
 24 //uN是匹配左边的顶点数，vN是匹配右边的顶点数
 25 //调用：res=hungary();输出最大匹配数
 26 //优点：适用于稠密图，DFS找增广路，实现简洁易于理解
 27 //时间复杂度:O(VE)
 28 //************************************************************ */
 29 //顶点编号从0开始的
 30 int uN,vN;//u,v的数目，使用前面必须赋值
 31 int g[MAXN][MAXN];//邻接矩阵
 32 int linker[MAXN];
 33 bool used[MAXN];
 34 bool dfs(int u)
 35 {
 36     for(int v = 0; v < vN;v++)
 37     if(g[u][v] && !used[v])
 38     {
 39         used[v] = true;
 40         if(linker[v] == -1 || dfs(linker[v]))
 41         {
 42             linker[v] = u;
 43             return true;
 44         }
 45     }
 46     return false;
 47 }
 48 int hungary()
 49 {
 50     int res = 0;
 51     memset(linker,-1,sizeof(linker));
 52     for(int u = 0;u < uN;u++)
 53     {
 54         memset(used,false,sizeof(used));
 55         if(dfs(u))res++;
 56     }
 57     return res;
 58 }
 59 int main()
 60 {
 61     int i,j,k;
 62     #ifndef ONLINE_JUDGE
 63     freopen("1.in","r",stdin);
 64     #endif
 65     while(scanf("%d%d",&n,&m)!=EOF)
 66     {
 67         if(n==0&&m==0)  break;
 68         scanf("%d",&k);
 69         int u,v;
 70         cl(a);
 71         while(k--)
 72         {
 73             scanf("%d%d",&u,&v);
 74             u--,v--;
 75             a[u][v]=-1;
 76         }
 77         int tot=0;
 78         for(i=0;i<n;i++)
 79         {
 80             for(j=0;j<m;j++)
 81             {
 82                 if(a[i][j]!=-1)
 83                 {
 84                     b[tot]=i*m+j;
 85                     a[i][j]=tot++;
 86                 }
 87             }
 88         }
 89         cl(g);
 90         for(i=0;i<n;i++)
 91         {
 92             for(j=0;j<m;j++)
 93             {
 94                 if(a[i][j]!=-1&&((i+j)%2==1))
 95                 {
 96                     int w=a[i][j];
 97                     if(i>0&&a[i-1][j]!=-1)
 98                         g[w][a[i-1][j]]=1;
 99                     if(i<n-1&&a[i+1][j]!=-1)
100                         g[w][a[i+1][j]]=1;
101                     if(j>0&&a[i][j-1]!=-1)
102                         g[w][a[i][j-1]]=1;
103                     if(j<m-1&&a[i][j+1]!=-1)
104                         g[w][a[i][j+1]]=1;
105                 }
106             }
107         }
108         vN=uN=tot;
109         printf("%d\n",hungary());
110         for(i=0;i<tot;i++)
111         {
112             if(linker[i]!=-1)
113             {
114                 int x1=b[i]/m;
115                 int y1=b[i]%m;
116                 int x2=b[linker[i]]/m;
117                 int y2=b[linker[i]]%m;
118                 printf("(%d,%d)--(%d,%d)\n",x1+1,y1+1,x2+1,y2+1);
119             }
120         }
121         printf("\n");
122     }
123 }