Cod sursa(job #1166183)

Utilizator StefansebiStefan Sebastian Stefansebi Data 3 aprilie 2014 12:40:00
Problema Flux maxim Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.71 kb
#include<fstream>
#include<vector>
#include<queue>
using namespace std;
ifstream fin("maxflow.in");
ofstream fout("maxflow.out");
vector <int> v[1001];
queue <int> q;
int viz[1001], c[1001][1001], f[1001][1001], t[1001];
int n, m, i, j, a, b, co, flow, fmin;


int bfs(int st){
    int i, j;
    while (!q.empty())
        q.pop();
    for (i = 1; i <= n; i++)
        viz[i] = 0;
    viz[st] = 1; q.push(st);
    while (!q.empty()){
        int nod = q.front(); q.pop();
        for (i = 0; i < v[nod].size(); i++){
            int fiu = v[nod][i];
            if (viz[fiu] == 0 and c[nod][fiu] != f[nod][fiu]){
                viz[fiu] = 1;
                t[fiu] = nod;
                q.push(fiu);
                if (fiu == n)
                    return 1;
            }
        }
    }
    return 0;
}

int main(){
    fin >> n >> m;
    for (i = 1; i <= m; i++){
        fin >> a >> b >> co;
        v[a].push_back(b);
        v[b].push_back(a);
        c[a][b] = co;
    }
    while (bfs(1)){
        for (i = 0; i < v[n].size(); i++){
            if (c[v[n][i]][n] != f[v[n][i]][n] and viz[v[n][i]]){
                t[n] = v[n][i];
                int nod = n;
                int fmin = 99999;
                while (nod != 1){
                    fmin = min(fmin, c[t[nod]][nod] - f[t[nod]][nod]);
                    nod = t[nod];
                }
                nod = n;
                while (nod != 1){
                    f[t[nod]][nod] += fmin;
                    f[nod][t[nod]] -= fmin;
                    nod = t[nod];
                }
                flow += fmin;
            }
        }
    }
    fout << flow << " ";
    fin.close();
    fout.close();
}