Cod sursa(job #1636137)

Utilizator arhivamanArhiva Man arhivaman Data 6 martie 2016 22:51:57
Problema Flux maxim Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 1.69 kb
#include <fstream>
#include <algorithm>
#include <queue>
#include <string.h>
#include <vector>
#define nMax 1005
#define INF 1 << 30
#define pb push_back
using namespace std;
ifstream fin("maxflow.in");
ofstream fout("maxflow.out");
vector<int> G[nMax];
queue<int> Q;
int n, m, C[nMax][nMax], F[nMax][nMax], TT[nMax], sink, source, mflow, viz[nMax];
void read()
{
    int a, b, c;
    fin>>n>>m;
    for(int i=1;i<=m;i++)
    {
        fin>>a>>b>>c;
        G[a].pb(b);
        G[b].pb(a);
        C[a][b]+=c;
    }
    source=1, sink=n;
}
bool dfs()
{
    memset(viz, 0, sizeof(viz));
    Q.push(source);
    viz[source]=1;
    while(!Q.empty())
    {
        int node=Q.front();
        Q.pop();
        for(vector<int>::iterator it=G[node].begin();it!=G[node].end();it++)
            if(viz[*it]==0 && C[node][*it]>F[node][*it])
            {
                TT[*it]=node;
                viz[*it]=1;
                Q.push(*it);
            }
    }
    return viz[sink];
}
void flow()
{
    while(dfs())
    {
        for(vector<int>::iterator it=G[sink].begin();it!=G[sink].end();it++)
        {
            int node=*it;
            if(C[node][sink]<=F[node][sink] || !viz[node])
                continue;
            TT[sink]=node;
            int fmin=INF;
            for(node=sink;node!=1;node=TT[node])
                fmin=min(fmin, C[TT[node]][node]-F[TT[node]][node]);
            for(node=sink;node!=1;node=TT[node])
            {
                F[TT[node]][node]+=fmin;
                F[node][TT[node]]-=fmin;
            }
            mflow+=fmin;
        }
    }
    fout<<mflow;
}
int main()
{
    read();
    flow();
    return 0;
}