Pagini recente » Cod sursa (job #2358359) | Cod sursa (job #425298) | Monitorul de evaluare | Cod sursa (job #635603) | Cod sursa (job #3188999)
#include <iostream>
#include <fstream>
#include <vector>
#include <climits>
using namespace std;
int HungarianAlgorithm(vector<vector<int>>& costMatrix) {
int n = costMatrix.size();
vector<int> u(n), v(n), p(n), way(n);
for (int i = 1; i < n; ++i) {
p[0] = i;
int j0 = 0;
vector<int> minv(n, INT_MAX);
vector<char> used(n, false);
do {
used[j0] = true;
int i0 = p[j0], delta = INT_MAX, j1;
for (int j = 1; j < n; ++j) {
if (!used[j]) {
int cur = costMatrix[i0][j] - u[i0] - v[j];
if (cur < minv[j]) {
minv[j] = cur;
way[j] = j0;
}
if (minv[j] < delta) {
delta = minv[j];
j1 = j;
}
}
}
for (int j = 0; j < n; ++j) {
if (used[j]) {
u[p[j]] += delta;
v[j] -= delta;
} else {
minv[j] -= delta;
}
}
j0 = j1;
} while (p[j0] != 0);
do {
int j1 = way[j0];
p[j0] = p[j1];
j0 = j1;
} while (j0);
}
return -v[0];
}
int main() {
ifstream fin("cc.in");
ofstream fout("cc.out");
int N;
fin >> N;
vector<vector<int>> costMatrix(N + 1, vector<int>(N + 1));
for (int i = 1; i <= N; ++i) {
for (int j = 1; j <= N; ++j) {
fin >> costMatrix[i][j];
}
}
fout << HungarianAlgorithm(costMatrix) << endl;
fin.close();
fout.close();
return 0;
}