Cod sursa(job #2267949)

Utilizator DavidLDavid Lauran DavidL Data 24 octombrie 2018 12:46:47
Problema Critice Scor 20
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 2.11 kb
#include <bits/stdc++.h>
#define pb push_back
using namespace std;
ifstream fi("critice.in");
ofstream fo("critice.out");

const int NMAX = 1005;
const int INF = (1e9);

int n, m;
vector <int> G[NMAX];
int C[NMAX][NMAX], F[NMAX][NMAX];
int cnt[NMAX][NMAX];
int p[NMAX];
bool viz[NMAX];
queue <int> Q;
vector <int> rez;
pair <int, int> muchie[10005];

int bfs()
{
    memset(viz, 0, sizeof(viz));
    memset(p, 0, sizeof(p));

    Q.push(1);
    viz[1] = 1;

    while (!Q.empty())
    {
        int curr = Q.front();
        Q.pop();

        if (curr == n)
            continue;

        for (auto v: G[curr])
        {
            if (F[curr][v] == C[curr][v] || viz[v])
                continue;

            Q.push(v);
            viz[v] = 1;
            p[v] = curr;
        }
    }
    return viz[n];
}

int getFlow()
{
    memset(F, 0, sizeof(F));

    int flow = 0;

    while (bfs())
        for (auto v: G[n])
        {
            if (F[v][n] == C[v][n] || !viz[v])
                continue;

            int minim = INF;
            for (int nod = n; nod != 1; nod = p[nod])
                minim = min(minim, C[p[nod]][nod] - F[p[nod]][nod]);

            for (int nod = n; nod != 1; nod = p[nod])
            {
                F[p[nod]][nod] += minim;
                F[nod][p[nod]] -= minim;
            }

            flow += minim;
        }
    return flow;
}

int main()
{
    fi >> n >> m;
    for (int i = 1; i <= m; i++)
    {
        int u, v, c;
        fi >> u >> v >> c;
        G[u].pb(v);
        G[v].pb(u);
        C[u][v] = C[v][u] = c;
        cnt[u][v] = cnt[v][u] = i;
        muchie[i] = {u, v};
    }

    int maxim = getFlow();

    for (int i = 1; i <= m; i++)
    {
        C[muchie[i].first][muchie[i].second]++;
        C[muchie[i].second][muchie[i].first]++;

        if (getFlow() > maxim)
            rez.pb(i);

        C[muchie[i].first][muchie[i].second]--;
        C[muchie[i].second][muchie[i].first]--;
    }

    fo << rez.size() << "\n";
    for (auto x: rez)
        fo << x << "\n";

    return 0;
}