Cod sursa(job #1653581)

Utilizator ZeBuGgErCasapu Andreas ZeBuGgEr Data 16 martie 2016 11:25:03
Problema Flux maxim Scor 100
Compilator cpp Status done
Runda Arhiva educationala Marime 1.87 kb
#include<cstdio>
#include<vector>
#include<queue>

using namespace std;

int n,m,x,y,z;
int flux[1010][1010],cap[1010][1010],tata[1010];
bool viz[1010];
vector<int> conn[1010];
queue<int> q;

int bfs()
{
    for(int i=1;i<=n;i++)
    {
        viz[i]=0;
    }
    viz[1]=1;
    q.push(1);
    int curr,next;
    while(!q.empty())
    {
        curr=q.front();
        for(vector<int>::iterator it=conn[curr].begin();it!=conn[curr].end();it++)
        {
            next=*it;
            if(flux[curr][next]<cap[curr][next]&&viz[next]==0)
            {
                viz[next]=1;
                q.push(next);
                tata[next]=curr;
            }
        }
        q.pop();
    }
    return viz[n];
}

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

    scanf("%d %d",&n,&m);

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

    int res=0;

    while(bfs())
    {
        for(vector<int>::iterator it=conn[n].begin();it!=conn[n].end();it++)
        {
            int curr=*it;
            if(flux[curr][n]<cap[curr][n]&&viz[curr]==1)
            {
                int mf=cap[curr][n]-flux[curr][n];
                while(tata[curr]!=0)
                {
                    mf=min(mf,cap[tata[curr]][curr]-flux[tata[curr]][curr]);
                    curr=tata[curr];
                }

                res+=mf;

                curr=*it;
                flux[curr][n]+=mf;
                flux[n][curr]-=mf;
                while(tata[curr]!=0)
                {
                    flux[tata[curr]][curr]+=mf;
                    flux[curr][tata[curr]]-=mf;
                    curr=tata[curr];
                }
            }
        }
    }

    printf("%d",res);
}