nysit 239 二分图最大匹配

题目链接：http://acm.nyist.net/JudgeOnline/problem.php?pid=239

题目大意：月老准备给n个女孩与n个男孩牵红线，成就一对对美好的姻缘。现在，由于一些原因，部分男孩与女孩可能结成幸福的一家，部分可能不会结成幸福的家庭。现在已知哪些男孩与哪些女孩如果结婚的话，可以结成幸福的家庭，月老准备促成尽可能多的幸福家庭，请你帮他找出最多可能促成的幸福家庭数量吧。

思路：即女孩作为一个顶点集，男孩作为一个顶点集，可以结成幸福家庭的男孩女孩之间连线，然后求这个二分图的最大匹配，很明显的匈牙利算法，但由于刚开始时采用邻接矩阵来存储导致TLE，很无语，在网上搜二分图匹配的其他算法，才发现原来匈牙利算法的时间复杂度确实够呛，O(n*E)，同时也发现了Hopcroft_Karp算法，复杂度O(n^0.5*E)，有所提高，虽然没仔细看算法的实现，就直接拷个模板过来，谁知交上还是一样的TLE，真的无语了，再后来，发现了网上有说邻接表比邻接矩阵省时的，就想着拿过来试试，谁知竟然A了，那个激动啊，哈哈~不多说了，直接上代码！

#include<stdio.h> #include<string.h> #define N 510 int map[N][N],used[N],match[N],g[N][N],n; int find(int k) { int i,x; for(i=1;i<=g[k][0];i++) /*相当于邻接表存储，这样就避免无效查找导致的超时*/ { x=g[k][i]; if(!used[x]) { used[x]=1; if(match[x]==0 || find(match[x])) { match[x]=k; return 1; } } } return 0; } int main() { int i,j,k,ncase,count,x,y; scanf("%d",&ncase); while(ncase--) { scanf("%d%d",&n,&k); memset(map,0,sizeof(map)); while(k--) { scanf("%d%d",&x,&y); map[x][y]=1; } for(i=1;i<=n;i++) { g[i][0]=0; /*记录与i相连的顶点数*/ for(j=1;j<=n;j++) if(map[i][j]) { g[i][0]++; g[i][g[i][0]]=j; /*存储与i相连的顶点编号*/ } } count=0; /*记录匹配数*/ memset(match,0,sizeof(match)); /*初始时匹配边为空*/ for(i=1;i<=n;i++) { memset(used,0,sizeof(used)); /*所有顶点置未访问标记*/ if(find(i)) count++; /*若找到i的匹配边，则计数器加一*/ } printf("%d/n",count); } return 0; }

下面是两个超时的算法

#include<stdio.h> #include<string.h> #define N 510 int map[N][N],used[N],match[N],n; int find(int k) { int i; for(i=1;i<=n;i++) /*此处循环次数较多，导致超时*/ { if(map[k][i] && !used[i]) { used[i]=1; if(match[i]==0 || find(match[i])) { match[i]=k; return 1; } } } return 0; } int main() { int i,k,ncase,count,x,y; scanf("%d",&ncase); while(ncase--) { scanf("%d%d",&n,&k); memset(map,0,sizeof(map)); while(k--) { scanf("%d%d",&x,&y); map[x][y]=1; } count=0; /*记录匹配数*/ memset(match,0,sizeof(match)); /*初始时匹配边为空*/ for(i=1;i<=n;i++) { memset(used,0,sizeof(used)); /*所有顶点置未访问标记*/ if(find(i)) count++; /*若找到i的匹配边，则计数器加一*/ } printf("%d/n",count); } return 0; }

#include <algorithm> #include <cstdio> #include <vector> using namespace std; const int N = 510; int distx[N], disty[N], matex[N], matey[N], q[N]; vector<int> g[N]; bool BFS(int n) { bool found = false; int fore = 0, rear = 0; for (int i = 0; i < n; ++i) { if (matex[i] == -1) q[rear++] = i;//加入队列 distx[i] = 0;//要匹配的一边 disty[i] = 0;//要匹配的另一边 } for (; fore < rear; ++fore)//队列不空 { int x = q[fore];//从队列从取出x，并扩展 for (vector<int>::const_iterator it = g[x].begin(); it != g[x].end(); ++it) {// if (!disty[*it])//选择未访问的点，因为是要求不相交的增广路径 { disty[*it] = distx[x]+1;//至少是1， if (matey[*it] == -1)//为被匹配的自由点，未被访问过，说明找到一条路径 found = true; else { distx[matey[*it]] = disty[*it]+1; q[rear++] = matey[*it]; } } } } return found; } bool DFS(int x) { for (vector<int>::const_iterator it = g[x].begin(); it != g[x].end(); ++it) if (disty[*it] == distx[x]+1) {//走了这条边，但是不一定是增广路径 disty[*it] = 0; if (matey[*it] == -1 || DFS(matey[*it])) return matex[x] = *it, matey[*it] = x, true; } return false; } int Hopcroft_Karp(int n) { int res = 0; fill_n(matex, n, -1); fill_n(matey, n, -1); while (BFS(n)) for (int i = 0; i < n; ++i)//可能存在多个增广路径 if (matex[i] == -1 && DFS(i))//找到了一个增广路径 ++res; return res; } int main() { int j,i,n,m,ncase,x,y; scanf("%d",&ncase); for(j=0;j<ncase;j++) { scanf("%d%d",&n,&m); for(i=0;i<m;i++) { scanf("%d%d",&x,&y); g[x].push_back(y); } printf("%d/n", Hopcroft_Karp(n)); } }