Pagini recente » Cod sursa (job #563192) | Cod sursa (job #1209088) | Cod sursa (job #1421655) | Cod sursa (job #2551069) | Cod sursa (job #927902)
Cod sursa(job #927902)
#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;
}