Cod sursa(job #475710)

Utilizator dornescuvladVlad Eugen Dornescu dornescuvlad Data 8 august 2010 02:09:02
Problema Flux maxim Scor 40
Compilator cpp Status done
Runda Arhiva educationala Marime 1.59 kb
using namespace std;
#include<iostream>
#include<fstream>
#include<vector>
#include<queue>
#include<cstring>
#include<bitset>
#define Nmax 1005
#define pb push_back
vector<short> G[Nmax];
ofstream fout("maxflow.out");
int C[Nmax][Nmax];
int F[Nmax][Nmax];
int qu[1005];
int ante[Nmax];
int viz[Nmax];
int N,M;
int pz;


inline int BFS()
{
	
	memset(viz, 0, sizeof(viz));
	int u,i,V;
    viz[1]=1;
    qu[0]=1;
    qu[1]=1;
    int alfa=1;

    while(alfa<=qu[0])
    {
        u=qu[alfa];
        alfa++;
        for(i=0;i<G[u].size();i++)
        { V=G[u][i];
          if(viz[V]==0)
           if(F[u][V]<C[u][V])
            {
                viz[V]=1;
                ante[V]=u;
                qu[0]++;
                qu[qu[0]]=V;
                if(V==N)
                return 1;
            }
        }


    }
    return 0;

}

inline void solve()
{int i,flow,fmin;

  for(flow=0;BFS();flow+=fmin)
  {
      for(i=N;i!=1;i=ante[i])
      fmin=min(fmin,C[ante[i]][i]-F[ante[i]][i]);
      for(i=N;i!=1;i=ante[i])
      {
          F[ante[i]][i]+=fmin;
          F[i][ante[i]]-=fmin;
      }

  }
  fout<<flow<<"\n";

}

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

     {fin>>x>>y>>z;
         /*cit(x);
         cit(y);
         cit(z);*/
     C[x][y]+=z;
     F[y][z]=0;
     G[x].pb(y);
     G[y].pb(x);
     }
     fin.close();

}

int main()
{ //freopen("maxflow.in","r",stdin);
    cit();
    solve();
    fout.close();
    return 0;
}