Cod sursa(job #2592520)

Utilizator popashtefan10Popa Stefan popashtefan10 Data 1 aprilie 2020 19:47:59
Problema Flux maxim Scor 0
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.77 kb
#include <iostream>
#include <cstdio>
#include <vector>
#include <queue>
#define NMAX 1000
#define MMAX 5000
#define INF 550000000

using namespace std;

bool dest;
int n, m, flux_tot, flux;
int f[NMAX + 5][NMAX + 5], c[NMAX + 5][NMAX + 5];
int from[NMAX + 5], mincf[NMAX + 5];
vector<int> v[NMAX + 5];
queue<int> q;

void update_flux(int nod, int flux) {
  if(!from[nod])
    return;
  update_flux(from[nod], flux);
  f[from[nod]][nod] += flux;
  f[nod][from[nod]] -= flux;
  mincf[nod] = min(c[from[nod]][nod] - f[from[nod]][nod], mincf[from[nod]]);
}

void bfs(int start) {
  q.push(start);

  while(!q.empty()) {
    int nod = q.front();
    q.pop();

    mincf[nod] = min(c[from[nod]][nod] - f[from[nod]][nod], mincf[from[nod]]);
    for(int x: v[nod]) {
      if(x == n && c[nod][n] - f[nod][n] > 0 && mincf[nod] != 0)
        dest = true;
      if(c[nod][x] - f[nod][x] == 0 || from[x] || x == n || x == 1)
        continue;
      from[x] = nod;
      q.push(x);
    }
  }
}

int main() {
  freopen("maxflow.in", "r", stdin);
  freopen("maxflow.out", "w", stdout);
  int x, y, z;

  scanf("%d %d", &n, &m);
  for(int i = 1; i <= m; i++) {
    scanf("%d %d %d", &x, &y, &z);
    c[x][y] = z;
    v[x].push_back(y);
    v[y].push_back(x);
  }

  dest = true;
  mincf[0] = c[0][1] = INF + 5;
  while(dest) {
    for(int i = 1; i <= n; i++)
      from[i] = 0;
    dest = true;
    //bfs(1);
    if(!dest)
      break;
    for(int nod: v[n]) {
      if(c[nod][n] == f[nod][n])
        continue;
      from[n] = nod;
      mincf[n] = min(mincf[nod], c[nod][n] - f[nod][n]);
      flux_tot += mincf[n];
      f[nod][n] += mincf[n];
      update_flux(n, mincf[n]);
    }
    dest = false;
  }

  printf("%d", flux_tot);

  return 0;
}