Cod sursa(job #2422014)

Utilizator Kln1000Ciobanu Bogdan Kln1000 Data 16 mai 2019 21:42:36
Problema Flux maxim Scor 40
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.85 kb
#include <iostream>
#include <fstream>
#include <queue>
#include <bitset>
#include <vector>

/// drumuri de augmentare
/// Ford Fulkerson algorithm
/// cuplaj maxim in graf bipartit
/// problema harta infoarena
using namespace std;

#define NMAX 1010

ifstream f ("maxflow.in");
ofstream t ("maxflow.out");

int main() {
    int C[NMAX][NMAX], n, m, total_flow = 0, parent[NMAX];
    bitset<NMAX> vaz;
    vector<int> v[NMAX];

    auto bfs = [&vaz, &v, &C, &parent, &n](int start, int finish) {
        queue<int> q;
        vaz.reset();

        vaz[start] = true;
        q.push(start);

        while (!q.empty()) {
            int current = q.front();
            q.pop();
            
            if (current == n)
                continue;
            
            for (auto i : v[current]) {
                if ((!vaz[i]) && C[current][i] > 0) {
                    vaz[i] = true;
                    parent[i] = current;
                    q.push(i);
                }
            }
        }

        return vaz[finish];
    };

    f >> n >> m;
    total_flow &= 0;
    for (int i = 0, x , y, c; i < m; ++i) {
        f >> x >> y >> c;
        C[x][y] = c;
        C[y][x] = 0;
        v[x].push_back(y);
        v[y].push_back(x);
    }
    
    for (int minflow; bfs(1, n);) {
        for(auto i : v[n]) {
            if (C[i][n] == 0 || !vaz[i])
                continue;
            parent[n] = i;

            minflow = 1e9;
            for (int nod = n; nod != 1; nod = parent[nod])
                minflow = C[parent[nod]][nod] < minflow ? C[parent[nod]][nod] : minflow;

            if (minflow == 0)
                continue;

            for (int nod = n; nod != 1; nod = parent[nod])
                C[parent[nod]][nod] -= minflow,
                C[nod][parent[nod]] += minflow;
            total_flow += minflow;
        }
    }

    t << total_flow;
    return 0;
}