求最大匹配的一种显而易见的算法是:先找出全部匹配,然后保留匹配数最多的。但是这个算法的复杂度为边数的指数级函数。因此,需要寻求一种更加高效的算法。 增广路的定义(也称增广轨或交错轨): 若P是图G中一条连通两个未匹配顶点的路径,并且属M的边和不属M的边(即已匹配和待匹配的边)在P上交替出现,则称P为相对于M的一条增广路径。 由增广路的定义可以推出下述三个结论: 1-P的路径长度必定为奇数,第一条边和最后一条边都不属于M。 2-P经过取反操作可以得到一个更大的匹配M’。 3-M为G的最大匹配当且仅当不存在相对于M的增广路径。 用增广路求最大匹配(称作匈牙利算法,匈牙利数学家Edmonds于1965年提出) 算法轮廓: (1)置M为空 (2)找出一条增广路径P,通过取反操作获得更大的匹配M’代替M (3)重复(2)操作直到找不出增广路径为止 程序清单: #include<stdio.h> #include<string.h> bool g[201][201]; int n,m,ans; bool b[201]; int link[201]; bool init() { int _x,_y; memset(g,0,sizeof(g)); memset(link,0,sizeof(link)); ans=0; if(scanf("%d%d",&n,&m)==EOF)return false; for(int i=1;i<=n;i++) { scanf("%d",&_x); for(int j=0;j<_x;j++) { scanf("%d",&_y); g[i][_y]=true; } } return true; } bool find(int a) { for(int i=1;i<=m;i++) { if(g[a][i]==1&&!b[i]) { b[i]=true; if(link[i]==0||find(link[i])) { link[i]=a; return true; } } } return false; } int main() { while(init()) { for(int i=1;i<=n;i++) { memset(b,0,sizeof(b)); if(find(i))ans++; } printf("%d\n",ans); } } 下面是同宿舍的小小牛写的,一起贴上吧,呵呵: #include <iostream> #include <fstream> using namespace std;
const int MAXN = 100; int uN, vN; bool g[MAXN][MAXN];//g[i][j] 表示 xi与yj相连 int xM[MAXN], yM[MAXN]; // 输出量 bool chk[MAXN]; //辅助量 检查某轮 y[v]是否被check
bool SearchPath(int u) { int v; for (v=0; v<vN; v++) { if (g[u][v] && !chk[v]) { chk[v] = true; if (yM[v] == -1 || SearchPath(yM[v])) { yM[v] = u; xM[u] = v; return true; } } } return false; }
int MaxMatch() { int u; int ret = 0; memset(xM, -1, sizeof(xM)); memset(yM, -1, sizeof(yM)); for (u=0; u<uN; u++) { if (xM[u] == -1) { memset(chk, false, sizeof(chk)); if (SearchPath(u)) ret++; } } return ret; }
int main() { int i, k; int tU, tV; ifstream cin("test.txt"); cin >> uN >> vN >> k; memset(g, false, sizeof(g)); for (i=0; i<k; i++) { cin >> tU >> tV; g[tU][tV] = true; } cout << MaxMatch() << endl; system("pause"); return 0; } |