Cod sursa(job #1966128)

Utilizator Horia14Horia Banciu Horia14 Data 14 aprilie 2017 21:54:57
Problema Flux maxim Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 2.25 kb
#include<cstdio>
#include<queue>
#define N_MAX 1000
using namespace std;

int c[N_MAX+1][N_MAX+1], f[N_MAX+1][N_MAX+1], T[N_MAX+1], n, m, flux, S, D;
struct nod
{
    int val;
    nod *next;
};

nod *L[N_MAX+1];

void Insert(int x, int y)
{
    nod *elem = new nod;
    elem->val = y;
    elem->next = L[x];
    L[x] = elem;
}

void Read()
{
    int i, a, b;
    FILE *fin = fopen("maxflow.in","r");
    fscanf(fin,"%d%d",&n,&m);
    for(i=1; i<=m; i++)
    {
        fscanf(fin,"%d%d",&a,&b);
        fscanf(fin,"%d",&c[a][b]);
        Insert(a,b);
        Insert(b,a);
    }
    fclose(fin);
}

bool BFS(int x)
{
    queue<int>Q;
    int i, Nod;
    nod *p;
    bool ok = 0;
    for(i=1; i<=n; i++)
        T[i] = 0;
    T[x] = -1;
    Q.push(x);
    while(!Q.empty())
    {
        Nod = Q.front();
        Q.pop();
        p = L[Nod];
        while(p != NULL)
        {
            if(!T[p->val] && c[Nod][p->val] > f[Nod][p->val])
            {
                if(p->val != D)
                {
                    T[p->val] = Nod;
                    Q.push(p->val);
                }
                else ok = 1;
            }
            p = p->next;
        }
    }
    return ok;
}

void MaxFlow()
{
    int drum = BFS(S), minim, Nod;
    nod *p;
    while(drum)
    {
        p = L[D];
        while(p != NULL)
        {
            if(T[p->val] && c[p->val][D] > f[p->val][D])
            {
                T[D] = p->val;
                minim = 0x7fffffff;
                for(Nod = D; Nod!=S; Nod = T[Nod])
                    if(minim > c[T[Nod]][Nod] - f[T[Nod]][Nod])
                        minim = c[T[Nod]][Nod] - f[T[Nod]][Nod];
                if(minim >= 0)
                {
                    flux += minim;
                    for(Nod = D; Nod!=S; Nod = T[Nod])
                    {
                        f[T[Nod]][Nod] += minim;
                        f[Nod][T[Nod]] -= minim;
                    }
                }
            }
            p = p->next;
        }
        drum = BFS(S);
    }
}

int main()
{
    FILE *fout = fopen("maxflow.out","w");
    Read();
    S = 1; D = n;
    MaxFlow();
    fprintf(fout,"%d\n",flux);
    fclose(fout);
    return 0;
}