Cod sursa(job #2246095)

Utilizator RaduXD1Nicolae Radu RaduXD1 Data 26 septembrie 2018 17:30:23
Problema Critice Scor 0
Compilator cpp Status done
Runda Arhiva de probleme Marime 1.9 kb
#include<fstream>
#include<bitset>
#include<queue>
#define DIM 1001
#define inf 2000000000
using namespace std;
ifstream fin("critice.in");
ofstream fout("critice.out");
int n,m,i,mini,ve,a,b,cap, sol[DIM],nod,j,va[10001],vb[10001],nr;
int c[DIM][DIM], f[DIM][DIM];
int t[DIM];
bitset<DIM> d;
deque<int> z;
vector<int> v[DIM];

int bfs()
{
    int i, nod,vecin;
    d.reset();
    z.push_back(1);
    d[1]=1;
    while(!z.empty())
    {
        nod=z[0];z.pop_front();
        for(i=0;i<v[nod].size();i++)
        {
            vecin=v[nod][i];
            if(c[nod][vecin]>f[nod][vecin]&&d[vecin]==0)
            {
                z.push_back(vecin);
                d[vecin]=1;t[vecin]=nod;
            }
        }
    }
    return d[n];
}

void solve(int start, int val)
{
    int i, nod,vecin;
    d.reset();
    z.push_back(start);
    d[start]=1;
    while(!z.empty())
    {
        nod=z[0];z.pop_front();
        for(i=0;i<v[nod].size();i++)
        {
            vecin=v[nod][i];
            if(!((f[vecin][nod]==c[vecin][nod])||(f[nod][vecin]==c[nod][vecin]))&&d[vecin]==0)
            {
                z.push_back(vecin);
                d[vecin]=1;t[vecin]=nod;
            }
        }
    }
    for(i=1;i<=n;i++)
        if(d[i])
            sol[i]+=val;
}

int main()
{
    fin>>n>>m;
    for(i=1;i<=m;i++)
    {
        fin>>a>>b>>cap;
        va[i]=a;vb[i]=b;
        v[a].push_back(b);
        v[b].push_back(a);
        c[a][b]=cap;
        c[b][a]=cap;
    }
    while(bfs()==1)
     {
        mini=inf;nod=n;
        while(t[nod])
        {
            mini=min(mini,c[t[nod]][nod]-f[t[nod]][nod]);
            nod=t[nod];
        }
        nod=n;
        while(t[nod])
        {
            f[t[nod]][nod]+=mini;
            f[nod][t[nod]]-=mini;
            nod=t[nod];
        }
    }
    solve(1, 1);solve(n, 2);

    return 0;
}