匈牙利算法，求二分图最大匹配。
若P是图G中一条连通两个未匹配顶点的路径，并且属于M的边和不属于M的边(即已匹配和待匹配的边)在P上交替出现，则称P为相对于M的一条增广路径。（M为一个匹配）
由增广路的定义可以推出下述三个结论：
- P的路径长度必定为奇数，第一条边和最后一条边都不属于M。所以Line 25-27从first part出发，不从二分图的另一部分出发。Line 12实现了交替出现的逻辑；node->neig匹配，当且仅当neig没有被其他点匹配，或者neig被first中的其他点matches[neig]匹配，并且从matches[neig]能够找到一条增广路径。这里就实现了交替的逻辑了。
- 将M和P进行异或操作(去同存异)可以得到一个更大的匹配M’。这是因为，属于M的边和不属于M的边交替出现，且第一和最后一条边都不属于M，所以增广路径中，不属于M的边比属于M的边多1，去同存异之后，一定会得到一个更大的匹配（加1了）。Line 13实现的是去同存异的逻辑。如果从node到neig存在一条增广路径，那么中间这些相同的部分直接省略。
- M为G的最大匹配当且仅当不存在M的增广路径。
#include <iostream>
#include <cstdio>
#include <vector> using namespace std; bool augment(vector<vector<int> > &adj, int node,
vector<bool> &visited, vector<int> &matches) {
for (auto neig : adj[node]) {
if (!visited[neig]) {
visited[neig] = true;
if (matches[neig] == - || augment(adj, matches[neig], visited, matches)) {
matches[neig] = node;
return true;
}
}
}
return false;
} int hungary(vector<vector<int> > &adj, vector<int> &first) {
vector<bool> visited;
vector<int> matches(adj.size(), -);
int count = ;
for (auto f : first) {
visited.assign(adj.size(), false);
if (augment(adj, f, visited, matches)) {
count++;
}
}
for (int i = ; i < adj.size(); ++i) {
cout << i << "<->" << matches[i] << endl;
}
return count;
} int main(int argc, char** argv) {
freopen("input.txt", "r", stdin);
int first, n, m;
cin >> n >> first >> m;
vector<vector<int> > adj(n);
vector<int> left;
for (int i = ; i < first; ++i) {
int l;
cin >> l;
left.push_back(l);
}
for (int i = ; i < m; ++i) {
int n1, n2;
cin >> n1 >> n2;
adj[n1].push_back(n2);
adj[n2].push_back(n1);
} cout << hungary(adj, left) << endl;
return ;
}
时间复杂度是O(VE)，空间复杂度感觉O(V)就行了啊，为什么其他人都说是O(V+E)？.