Cod sursa(job #1831521)

Utilizator NicuBaciuNicu Baciu NicuBaciu Data 18 decembrie 2016 11:33:06
Problema Flux maxim Scor 10
Compilator cpp Status done
Runda Arhiva educationala Marime 2.39 kb
#include <fstream>
#include <vector>
#include <queue>
#include <climits>

using namespace std;

ifstream fin("maxflow.in");
ofstream fout("maxflow.out");

struct muchie
{
    int n1;
    int n2;
    int cap;
    int f=0;
    int vr;
    bool inv;
};

struct nod
{
    bool viz=false;
    vector <muchie> v;
};

struct aux
{
    int nod0;
    int poz;
};

nod g[1001];
queue <aux> q;

int n, m;
int cr=INT_MAX;

bool sem=false;

void dfs(int nod)
{
    if(nod!=n)
    {
        for(int i=0; i<g[nod].v.size(); i++)
        {
            int test1=g[nod].v[i].vr;
            int test2=g[nod].v[i].n2;
            if(g[nod].v[i].vr>0)
            {
                if(g[g[nod].v[i].n2].viz==false)
                {
                    cr=min(cr, g[nod].v[i].vr);

                    aux t;
                    t.nod0=nod; t.poz=i;

                    q.push(t);
                    g[g[nod].v[i].n2].viz=true;

                    dfs(g[nod].v[i].n2);

                    if(sem==false)
                    {
                        g[g[nod].v[i].n2].viz=false;
                        q.pop();
                    }
                    else
                        return;
                }
            }
        }
    }
    else
    {
        sem=true;
        return;
    }
}

void flux()
{
    while(!q.empty())
    {
        aux x;

        x=q.front(); q.pop();

        if(g[x.nod0].v[x.poz].inv==false)
        {
            g[x.nod0].v[x.poz].f+=cr;
            g[x.nod0].v[x.poz].vr=g[x.nod0].v[x.poz].cap-g[x.nod0].v[x.poz].f;
        }
        if(g[x.nod0].v[x.poz].inv==true)
        {
            g[x.nod0].v[x.poz].f-=cr;
            g[x.nod0].v[x.poz].vr=g[x.nod0].v[x.poz].f;
        }
    }
}

int fluxmaxim=0;

int main()
{
    fin >> n >> m;

    for(int i=1; i<=m; i++)
    {
        int x, y, c;
        muchie m0;

        fin >> x >> y >> c;

        m0.n1=x; m0.n2=y; m0.cap=c; m0.vr=m0.cap-m0.f; m0.inv=false;
        g[x].v.push_back(m0);

        m0.n1=y; m0.n2=x; m0.cap=0; m0.vr=m0.f; m0.inv=true;
        g[y].v.push_back(m0);
    }

    dfs(1);

    while(!q.empty())
    {
        flux();
        fluxmaxim+=cr;
        cr=INT_MAX;
        sem=false;
        for(int i=1; i<=n; i++)
            g[i].viz=false;

        dfs(1);
    }

    fout << fluxmaxim;

    return 0;
}