Cod sursa(job #1688246)

Utilizator TibixbAndrei Tiberiu Tibixb Data 13 aprilie 2016 12:49:08
Problema Flux maxim de cost minim Scor 50
Compilator cpp Status done
Runda Arhiva educationala Marime 1.63 kb
#include<fstream>
#include<vector>
#include<cstring>
#include<queue>
#define inf 0x7f7f7f7f
using namespace std;
vector<int> L[352];
bool u[352];
queue<int> q;
int T[352], n, m, s, d, i, x, y, c, z, C[352][352], F[352][352], Z[352][352], nod, vecin, cost, D[352], minim;
int bellmanford()
{
    memset(D, 127, sizeof(D));
    D[s]=0;
    memset(u, 0, sizeof(u));
    q.push(s);
    u[s]=1;
    while(!q.empty())
    {
        nod=q.front();
        for(i=0; i<L[nod].size(); i++)
        {
            vecin=L[nod][i];
            if(D[vecin]>D[nod]+Z[nod][vecin] && C[nod][vecin]-F[nod][vecin]>0)
            {
                D[vecin]=D[nod]+Z[nod][vecin];
                T[vecin]=nod;
                if(u[vecin]==0)
                {
                    q.push(vecin);
                    u[vecin]=1;
                }
            }
        }
        q.pop();
        u[nod]=0;
    }
    return D[d]!=inf;
}
void fmcm()
{
    while(bellmanford())
    {
        minim=C[T[d]][d]-F[T[d]][d];
        for(nod=T[d]; T[nod]!=0; nod=T[nod])
        {
            minim=min(minim, C[T[nod]][nod]-F[T[nod]][nod]);
        }
        for(nod=d; T[nod]!=0; nod=T[nod])
        {
            cost+=minim*Z[T[nod]][nod];
            F[T[nod]][nod]+=minim;
            F[nod][T[nod]]-=minim;
        }
    }
}
ifstream in("fmcm.in");
ofstream out("fmcm.out");
int main()
{
    in>>n>>m>>s>>d;
    for(i=1; i<=m; i++)
    {
        in>>x>>y>>c>>z;
        L[x].push_back(y);
        C[x][y]=c;
        Z[x][y]=z;
        L[y].push_back(x);
        Z[y][x]=-z;
    }
    fmcm();
    out<<cost;
    return 0;
}