Cod sursa(job #514080)

Utilizator APOCALYPTODragos APOCALYPTO Data 17 decembrie 2010 18:46:34
Problema Flux maxim Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.75 kb
using namespace std;
#include<iostream>
#include<fstream>
#include<vector>
#include<cstring>
#define pb push_back
#define nmax 1005
#define oo 0x3f3f3f3f
vector<int> G[nmax];
int C[nmax][nmax];
int F[nmax][nmax];
int viz[nmax];
int N,M;
int ante[nmax];
int q[nmax];

ofstream fout("maxflow.out");

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

void solve()
{int flow,fmin,i;
vector<int>::iterator it;

    for(flow=0;BFS();)
    {
        fmin=oo;
        for(it=G[N].begin();it<G[N].end();it++)
        {
            if(viz[*it]==1 && C[*it][N]>0)
            {
                fmin=C[*it][N]-F[*it][N];
                for(i=*it;i!=1;i=ante[i])
                 fmin=min(fmin,C[ante[i]][i]-F[ante[i]][i]);
                F[*it][N]+=fmin;
                F[N][*it]-=fmin;
                for(i=*it;i!=1;i=ante[i])
                {
                    F[ante[i]][i]+=fmin;
                    F[i][ante[i]]-=fmin;

                }

                flow+=fmin;
            }
        }
    }
    fout<<flow<<"\n";

}

void cit()
{int i,x,y,z;
    ifstream fin("maxflow.in");
    fin>>N>>M;
    for(i=1;i<=M;i++)
    {
        fin>>x>>y>>z;
        G[x].pb(y);
        G[y].pb(x);
        C[x][y]+=z;
        F[x][y]=F[y][x]=0;
    }
    fin.close();
}

int main()
{
    cit();
    solve();
    fout.close();
    return 0;

}