Cod sursa(job #2974084)

Utilizator razvanalexrotaruRazvan Alexandru Rotaru razvanalexrotaru Data 3 februarie 2023 09:15:04
Problema Flux maxim Scor 70
Compilator cpp-64 Status done
Runda Arhiva educationala Marime 1.84 kb
#include <bits/stdc++.h>
#define MAX 1005
using namespace std;
map<pair<int,int>,int>mapp;
vector<int>v[MAX],dead_ends;
int n,prag,nvl[MAX];
ifstream cin("maxflow.in");
ofstream cout("maxflow.out");
void nivelare(int nod);
int dfs(int nod,int minn);
int main()
{
    ifstream cin("maxflow.in");
    ofstream cout("maxflow.out");
    int m,i,x,y,z,r=0;

    cin>>n>>m;
    for(i=1;i<=m;i++)
    {
        cin>>x>>y>>z;
        v[x].push_back(y);
        mapp[{x,y}]=z;
        if(mapp[{y,x}]==0)
            v[y].push_back(x);
    }
    nivelare(1);
    while(nvl[n]!=0)
    {
        while(prag!=0)
        {
            x=dfs(1,INT_MAX);
            if(x!=0)
                r+=x;
            else
            {
                prag/=2;
            }
        }
        nivelare(1);
    }
    cout<<r;
    return 0;
}
void nivelare(int nod)
{
    int maxx=0;
    queue<int>q;
    for(int i=1;i<=n;i++)
    {
        nvl[i]=0;
    }
    nvl[nod]=1;
    q.push(nod);
    while(!q.empty())
    {
        nod=q.front();
        q.pop();
        for(int it:v[nod])
        {
            if(nvl[it]==0)
                maxx=max(maxx,mapp[{nod,it}]);
            if(nvl[it]==0 && mapp[{nod,it}]>0)
            {
                nvl[it]=nvl[nod]+1;
                q.push(it);
            }
        }
    }
    prag=(1<<int(log2(maxx)));
}

int dfs(int nod, int minn)
{
    if(nod==n)
        return minn;
    for(int it:v[nod])
    if(nvl[it]==nvl[nod]+1 && mapp[{nod,it}]>=prag)
        {
            int x=dfs(it,min(minn,mapp[{nod,it}]));
            if(x>0)
            {
                mapp[{nod,it}]-=x;
                mapp[{it,nod}]+=x;
                return x;
            }
            else
            {
                dead_ends.push_back(it);
            }
        }
    return 0;
}