Cod sursa(job #547125)

Utilizator APOCALYPTODragos APOCALYPTO Data 5 martie 2011 21:43:21
Problema Flux maxim Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.86 kb
using namespace std;
#include<iostream>
#include<fstream>
#include<vector>
#define oo 0x3f3f3f3f
#define nmax 1010
#define pb push_back
ofstream fout("maxflow.out");
struct nod{
    int lg,c;
};
int viz[nmax],ante[nmax],q[nmax*nmax];
int C[nmax][nmax], F[nmax][nmax];
int N,M;
vector<nod> G[nmax];
bool BFS()
{
    int i;
    for(i=0;i<=N;i++)
    {
        viz[i]=0;
        q[i]=0;
        ante[i]=0;
    }
    viz[1]=1;
    q[++q[0]]=1;

    int u,front=1;
    vector<nod>::iterator it;
    while(front<=q[0])
    {
        u=q[front++];
        for(it=G[u].begin();it<G[u].end();it++)
        {
            if(C[u][it->lg]>F[u][it->lg] && !viz[it->lg])
            {
                viz[it->lg]=1;
                ante[it->lg]=u;
                q[++q[0]]=it->lg;
            }
        }
    }
    return viz[N];
}

void solve()
{
    int fmin=oo,flow=0,i;
    vector<nod>::iterator it;
    for(;BFS();)
    {
        for(it=G[N].begin();it<G[N].end();it++)
        {
          if(viz[it->lg])
          {
              fmin=C[it->lg][N]-F[it->lg][N];
              for(i=it->lg;i!=1;i=ante[i])
              {
                  fmin=min(fmin,C[ante[i]][i]-F[ante[i]][i]);
              }

              F[it->lg][N]+=fmin;
              F[N][it->lg]-=fmin;

              for(i=it->lg;i!=1;i=ante[i])
              {
                  F[ante[i]][i]+=fmin;
                  F[i][ante[i]]-=fmin;
              }
              flow+=fmin;
          }
        }
    }
    fout<<flow<<"\n";
}

void cit()
{
    ifstream fin("maxflow.in");
    int i,x,y,z;
    fin>>N>>M;
    for(i=1;i<=M;i++)
    {
        fin>>x>>y>>z;

        C[x][y]+=z;
        G[x].pb((nod){y,z});
        G[y].pb((nod){x,z});


    }

    fin.close();
}

int main()
{
    cit();

    solve();

    fout.close();
    return 0;

}