Cod sursa(job #1914924)

Utilizator ajeccAjechiloae Eugen ajecc Data 8 martie 2017 18:57:40
Problema Flux maxim Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 4.19 kb
#include <bits/stdc++.h>
#define for0(i,n) for(int i=0; i<n; i++)
#define for1(i,n) for(int i=1; i<=n; i++)
#define pb push_back
#define mp make_pair
#define all(v) v.begin(), v.end()
#define V vector<int>
#define VP vector<pair<int, int> >
#define clr(A,x) memset(A, x, sizeof(A))
#define cpy(A,B) memcpy(A, B, sizeof(B))
#define g(stream, s) getline(stream, s)
#define FASTIO ios_base::sync_with_stdio(0)
typedef long long ll;
typedef unsigned long long ull;
using namespace std;

const int CHARMAX = (1 << 10);
const int INFINT = 2000000000 + 100;
const ll INFLL = (1LL << 62);

class InParsare
{
private:
    FILE *fin;
    char *buffer;
    size_t index_of_buffer;

    char read_character()
    {
        index_of_buffer++;
        if(index_of_buffer == CHARMAX)
        {
            fread(buffer, 1, CHARMAX, fin);
            index_of_buffer = 0;
        }
        return buffer[index_of_buffer];
    }

public:
    InParsare(const char *name)
    {
        fin = fopen(name, "r");
        buffer = new char[CHARMAX]();
        index_of_buffer = CHARMAX - 1;
    }

    template<class T>
    InParsare &operator >> (T &n)
    {
        char c;
        while(!isdigit(c = read_character()) && c != '-');

        int sgn = 1;
        if(c == '-')
        {
            n = 0;
            sgn = -1;
        }
        else n = c - '0';

        while(isdigit(c = read_character()))
            n = n * 10 + c - '0';
        n *= sgn;

        return *this;
    }
};

class OutParsare
{
private:
    FILE *fout;
    char *buffer;
    size_t index_of_buffer;

    void write_character(char character)
    {
        if(index_of_buffer == CHARMAX)
        {
            fwrite(buffer, 1, CHARMAX, fout);
            index_of_buffer = 0;
            buffer[index_of_buffer++] = character;
        }
        else buffer[index_of_buffer++] = character;
    }

public:

    OutParsare(const char *name)
    {
        fout = fopen(name, "w");
        buffer = new char[CHARMAX]();
        index_of_buffer = 0;
    }

    ~OutParsare()
    {
        fwrite(buffer, 1, index_of_buffer, fout);
        fclose(fout);
    }
    template<class T>
    OutParsare &operator << (T n)
    {
        if(typeid(T).name() == typeid(char).name())
        {
            write_character(n);
            return *this;
        }

        if(n <= 9)
            write_character(n + '0');
        else
        {
            (*this) << (n / 10);
            write_character(n % 10 + '0');
        }

        return *this;
    }

};

InParsare fin("maxflow.in"); /// merg numai pt numere
OutParsare fout("maxflow.out"); /// merg numai pt numere + caractere individuale, nu stringuri
//ifstream fin("maxflow.in");
//ofstream fout("maxflow.out");

const int NMAX = 1e3 + 1;
const int MOD = 666013;
const double EPS = 0.000000001;

int n, m;
V graf[NMAX];
int cap[NMAX][NMAX];
int flx[NMAX][NMAX];
queue<int> coada;
bool viz[NMAX];
int tata[NMAX];


bool bfs()
{
    for1(i, n)
    viz[i] = 0;
    coada.push(1);
    viz[1] = 1;

    while(!coada.empty())
    {
        int nod = coada.front();
        coada.pop();
        if(nod != n)
            for(auto i: graf[nod])
                if(cap[nod][i] > flx[nod][i] && !viz[i])
                {
                    tata[i] = nod;
                    viz[i] = 1;
                    coada.push(i);
                }
    }

    return viz[n];
}

int sol;

int main()
{
    fin >> n >> m;
    for1(i, m)
    {
        int a, b;
        fin >> a >> b;
        graf[a].pb(b);
        graf[b].pb(a);
        fin >> cap[a][b];
    }

    while(bfs())
        for(auto i: graf[n])
            if(viz[i] && cap[i][n] > flx[i][n])
            {
                int mmin = INFINT;
                tata[n] = i;

                for(int nod = n; nod != 1; nod = tata[nod])
                    mmin = min(cap[tata[nod]][nod] - flx[tata[nod]][i], mmin);
                sol += mmin;
                for(int nod = n; nod != 1; nod = tata[nod])
                    flx[nod][tata[nod]] -= mmin,
                    flx[tata[nod]][nod] += mmin;
            }
    fout << sol;

    return 0;
}