Cod sursa(job #2693050)

Utilizator monica_LMonica monica_L Data 4 ianuarie 2021 17:20:28
Problema Critice Scor 100
Compilator cpp-64 Status done
Runda Arhiva de probleme Marime 2.11 kb
#include <fstream>
#include <vector>
#include <queue>
#include <string.h>
#define nmax 1005
using namespace std;
ifstream f("critice.in");
ofstream g("critice.out");
int n,m,x[nmax*10],y[nmax*10],c,C[nmax][nmax],F[nmax][nmax],t[nmax],T[nmax],FC,FL;
vector <int> v[nmax];
vector <int> sol;
int bf()
{
    int nod,vec,i;
    queue <int> q;
    memset(t,0,sizeof(t));
    t[1]=-1;
    q.push(1);
    while(!q.empty())
    {
        nod=q.front();
        q.pop();
        for(i=0;i<v[nod].size();i++)
        {
            vec=v[nod][i];
            if(t[vec] || C[nod][vec]==F[nod][vec]) continue;
            t[vec]=nod;
            q.push(vec);
        }
    }
    if(t[n]) return 1;
    return 0;
}
void dinic()
{
    int j,i,vec;
    while(bf())
    {
        for(i=0;i<v[n].size();i++)
        {
            vec=v[n][i];
            if(!t[vec] || C[vec][n]==F[vec][n]) continue;
            FC=C[vec][n]-F[vec][n];
            for(j=vec;j>1;j=t[j])
                FC=min(FC,C[t[j]][j]-F[t[j]][j]);
            if(FC==0) continue;
            F[vec][n]+=FC;
            F[n][vec]-=FC;
            for(j=vec;j>1;j=t[j])
            {
                F[t[j]][j]+=FC;
                F[j][t[j]]-=FC;
            }
        }
    }
}
int main()
{
    int i,nod,vec;
    f>>n>>m;
    for(i=1;i<=m;i++)
    {
        f>>x[i]>>y[i]>>c;
        v[x[i]].push_back(y[i]);
        v[y[i]].push_back(x[i]);
        C[x[i]][y[i]]=c;
        C[y[i]][x[i]]=c;
    }
    
    dinic();
    bf();
   
    queue <int> q;
    T[n]=-1;
    q.push(n);
    while(!q.empty())
    {
        nod=q.front();
        q.pop();
        for(i=0;i<v[nod].size();i++)
        {
            vec=v[nod][i];
            if(T[vec]||C[vec][nod]==F[vec][nod]||vec==n) continue;
            T[vec]=nod;
            q.push(vec);
        }
    }

    for(i=1;i<=m;i++)
        if((t[x[i]] && T[y[i]] && C[x[i]][y[i]]==F[x[i]][y[i]]) || (T[x[i]] && t[y[i]] && C[y[i]][x[i]]==F[y[i]][x[i]]))
        sol.push_back(i);
    g<<sol.size()<<'\n';
    for(i=0;i<sol.size();i++) g<<sol[i]<<'\n';
    return 0;
}