二分图的最大独立集 最大匹配解题 Hopcroft-Karp算法

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<queue>
 6 using namespace std;
 7 
 8 const int N=1005,INF=0x3f3f3f3f;
 9 int bmap[N][N],cx[N],cy[N],dx[N],dy[N];
10 bool bmask[N];
11 int nx,ny,dis,ans;
12 bool searchpath()
13 {
14     queue<int> q;
15     dis=INF;
16     memset(dx,-1,sizeof(dx));
17     memset(dy,-1,sizeof(dy));
18     for(int i=1;i<=nx;i++)
19     {
20         if(cx[i]==-1){ q.push(i); dx[i]=0; }
21         while(!q.empty())
22         {
23             int u=q.front(); q.pop();
24             if(dx[u]>dis) break;
25             for(int v=1;v<=ny;v++)
26             {
27                 if(bmap[u][v]&&dy[v]==-1)
28                 {
29                     dy[v]= dx[u] + 1;
30                     if(cy[v]==-1) dis=dy[v];
31                     else
32                     {
33                         dx[cy[v]]= dy[v]+1;
34                         q.push(cy[v]);
35                     }
36                 }
37             }
38         }
39     }
40     return dis!=INF;
41 }
42 int findpath(int u)
43 {
44     for(int v=1;v<=ny;v++)
45     {
46         if(!bmask[v]&&bmap[u][v]&&dy[v]==dx[u]+1)
47         {
48             bmask[v]=1;
49             if(cy[v]!=-1&&dy[v]==dis) continue;
50             if(cy[v]==-1||findpath(cy[v]))
51             {
52                 cy[v]=u; cx[u]=v;
53                 return 1;
54             }
55         }
56     }
57     return 0;
58 }
59 void maxmatch()
60 {
61     ans=0;
62     memset(cx,-1,sizeof(cx));
63     memset(cy,-1,sizeof(cy));
64     while(searchpath())
65     {
66         memset(bmask,0,sizeof(bmask));
67         for(int i=1;i<=nx;i++)
68             if(cx[i]==-1) ans+=findpath(i);
69     }
70 }
71 int main()
72 {
73     //freopen("test.txt","r",stdin);
74     int m,i,j,k=1,cas;
75     while(scanf("%d%d%d",&nx,&ny,&m)!=EOF)
76     {
77         if(!nx) break;
78         for(i=1;i<=nx;i++)
79             for(j=1;j<=ny;j++)
80                 bmap[i][j]=1;
81         while(m--)
82         {
83             scanf("%d%d",&i,&j);
84             bmap[i][j]=0;
85         }
86         maxmatch();
87         printf("Case %d: ",k++);
88         printf("%d\n",nx+ny-ans);
89     }
90     return 0;
91 }