Cod sursa(job #3188999)

Utilizator ilieflorin2110Ilie Florin Alexandru ilieflorin2110 Data 4 ianuarie 2024 12:47:25
Problema Cc Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 1.77 kb
#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;
}