Pagini recente » Cod sursa (job #1374066) | Cod sursa (job #2873553) | Cod sursa (job #203752) | Cod sursa (job #1608238) | Cod sursa (job #671874)
Cod sursa(job #671874)
#include<fstream>
#include<vector>
#include<queue>
#include<bitset>
#define DIM 1001
#define pb push_back
#define INF 0x3f3f3f3f
using namespace std;
ifstream in("critice.in");
ofstream out("critice.out");
int N, M;
int C[DIM][DIM], F[DIM][DIM], T[DIM], edges[DIM][DIM];
int T1[DIM], T2[DIM];
vector< int > G[DIM], sol;
bitset<DIM> viz;
void Read()
{
int i, x, y, c;
in >> N >> M;
for(i = 1; i <= M; i++)
{
in >> x >> y >> c;
G[x].pb(y);
G[y].pb(x);
C[x][y] = C[y][x] = c;
edges[x][y] = edges[y][x] = i;
}
}
int BFs(int S, int D, int *T )
{
queue<int> Q;
viz.reset();
Q.push(S);
viz.set(S);
while( !Q.empty() )
{
int nod = Q.front();
Q.pop();
for(vector< int >:: iterator it = G[nod].begin(); it != G[nod].end(); ++it)
if( F[nod][*it] < C[nod][*it] && !viz[*it] )
{
Q.push(*it);
viz.set(*it);
T[*it] = nod;
}
}
return viz[D];
}
void Solve()
{
int fmin, i, nod ;
while( BFs(1, N, T) )
{
fmin = INF;
for(nod = N; nod != 1; nod = T[nod] )
fmin = min( fmin, C[T[nod]][nod] - F[T[nod]][nod] );
for(nod = N; nod != 1; nod = T[nod] )
{
F[T[nod]][nod] += fmin;
F[nod][T[nod]] -= fmin;
}
}
BFs(1, N, T1);
BFs(N, 1, T2);
for(i = 1; i <= N; i++)
for(vector< int > :: iterator it = G[i].begin(); it != G[i].end(); ++it)
if( T1[i] && T2[*it] )
sol.pb( edges[i][*it] );
}
int main()
{
Read();
Solve();
out << sol.size() << '\n';
for(vector<int>::iterator it = sol.begin(); it != sol.end(); ++it)
out << *it << '\n';
return 0;
}