Cod sursa(job #927902)

Utilizator vendettaSalajan Razvan vendetta Data 26 martie 2013 09:27:30
Problema Critice Scor 80
Compilator cpp Status done
Runda Arhiva de probleme Marime 3.62 kb
#include <iostream>
#include <fstream>
#include <algorithm>
#include <vector>

using namespace std;

//ifstream f("critice.in");
ofstream g("critice.out");

#define nmax 1005
#define Mmax 10005
#define Cifmax (1<<14)

int n, m, q[nmax+Mmax], capac[nmax][nmax], flux[nmax][nmax], t[nmax];
bool viz[nmax];
int alb[nmax], negru[nmax];
pair<int,int> edge[Mmax];
vector<int> gf[nmax];
int rez[Mmax];
int pz = 0;
char S[Cifmax];

void buf(int &nr){
    nr = 0;
    for(; S[pz]<'0' || S[pz] >'9';){
        ++pz; if (pz == Cifmax){
            fread(S, 1, Cifmax, stdin);
            pz = 0;
        }
    }

    for(; S[pz]>='0' && S[pz]<='9';){
        nr = nr * 10 + (S[pz] - '0');
        ++pz; if (pz == Cifmax){
            fread(S, 1, Cifmax, stdin);
            pz = 0;
        }
    }
}

void citeste(){
    buf(n); buf(m);
    int x, y;
    for(int i=1; i<=m; ++i){
        buf(x); buf(y);
        gf[x].push_back(y);
        gf[y].push_back(x);
        buf(capac[x][y]);
        capac[y][x] = capac[x][y];
        edge[i] = make_pair(x, y);
    }
}

inline bool bfs(int S, int D){
    for(int i=0; i<=D; ++i) viz[i] = 0;
    int st = 1; int dr = 0;
    q[++dr] = S; viz[S] = 1;

    for(; st<=dr; ){
        int nod = q[st]; ++st;
        for(int i=0; i<gf[nod].size(); ++i){
            int vc = gf[nod][i];
            if (viz[vc] == 0 && capac[nod][vc] > flux[nod][vc]){
                t[vc] = nod;
                viz[vc] = 1;
                if(vc == D) return 1;
                q[++dr] = vc;
            }
        }
    }
    return 0;
}

void bagaFlux(){
    int S = 1; int D = n;
    for(; bfs(S, D); ){
        int Min = (1<<30);
        for(int i=D; i!=S; i=t[i]){
            //muchia t[i] -> i
            Min = min(Min, capac[ t[i] ][i] - flux[ t[i] ][i]);
        }
        //bag fluxul
        for(int i=D; i!=S; i=t[i]){
            flux[ t[i] ][i] += Min;
            flux[i][ t[i] ] -= Min;
        }
    }
}

void coloreaza(int S, int color[], int cod){
    int st = 1; int dr = 0;
    q[++dr] = S; color[S] =1;

    for(; st<=dr; ){
        int nod = q[st]; ++st;
        for(int i=0; i<gf[nod].size(); ++i){
            int vc = gf[nod][i];
            if (cod == 0){// luate in ordinea in care vin
                if (color[vc] == 0 && flux[nod][vc] < capac[nod][vc]){
                    color[vc] = 1;
                    q[++dr] = vc;
                }
            }else {// trebe in ordine inversa ca vin din destinatie
                if (color[vc] == 0 && flux[vc][nod] < capac[vc][nod]){
                    color[vc] = 1;
                    q[++dr] = vc;
                }
            }
        }
    }
}



void rezolva(){
    // dupa ce bag flux ma intereseaza muchiile critice acestea vor fi acele care :
    // au fluxul == capac si sunt accesibile din sursa mergand dooar pe muchii nesaturate
    // la fel si pentru destinatie

    bagaFlux();
    coloreaza(1, alb, 0);
    coloreaza(n, negru, 1);

    for(int i=1; i<=m; ++i){
        int X = edge[i].first; int Y = edge[i].second;
        if (flux[X][Y] == capac[X][Y] || flux[Y][X] == capac[Y][X]){
            if ( (alb[X] == 1 && negru[Y] == 1) || (negru[X] == 1 && alb[Y] == 1) ){// e critica
                rez[++rez[0]] = i;
            }
        }
    }

    //cout << rez[0] << "\n";
    g << rez[0] << "\n";
    for(int i=1; i<=rez[0]; ++i){
        //cout << rez[i] << "\n";
        g << rez[i] << "\n";
    }
}

int main(){
    freopen("critice.in", "r", stdin);

    citeste();
    rezolva();

    //f.close();
    g.close();

    return 0;
}