Pagini recente » Cod sursa (job #472601) | Cod sursa (job #1688199) | Cod sursa (job #1865561) | Cod sursa (job #183185) | Cod sursa (job #719320)
Cod sursa(job #719320)
#include <cstdio>
#include <algorithm>
#define N 1005
#define M 10005
#define inf 1999999999
using namespace std;
struct nod{ int val; nod *urm; }*G[N];
int n,m,cap[N][N],flux[N][N],tata[N],c[N],viz[N],muc[N][N];
bool fol[N];
int pas,k,sol[M];
void read()
{ int x,y,z,i;
nod *aux;
freopen("critice.in","r",stdin); scanf("%d %d\n",&n,&m);
for(i=1;i<=m;++i)
{
scanf("%d %d %d\n",&x,&y,&z);
cap[x][y]=cap[y][x]=z;
muc[x][y]=muc[y][x]=i;
aux=new nod; aux->val=y; aux->urm=G[x]; G[x]=aux;
aux=new nod; aux->val=x; aux->urm=G[y]; G[y]=aux;
}
fclose(stdin);
}
int bf()
{ int p,u;
nod *aux;
p=u=1; c[1]=1; viz[1]=pas; tata[1]=0;
while(p<=u)
{
if(c[p]!=n)
{
for(aux=G[c[p]];aux!=NULL;aux=aux->urm)
if(viz[aux->val]!=pas&&cap[c[p]][aux->val]>flux[c[p]][aux->val])
{
c[++u]=aux->val; tata[aux->val]=c[p]; viz[aux->val]=pas;
}
}
++p;
}
if(viz[n]==pas)return 1;
return 0;
}
inline int min(int x,int y)
{
if(x<y)return x;
return y;
}
void maxflow()
{ nod *aux;
int i,fm;
pas=1;
while(bf())
{
for(aux=G[n];aux!=NULL;aux=aux->urm)
if(viz[aux->val]==pas&&cap[aux->val][n]>flux[aux->val][n])
{
tata[n]=aux->val; fm=inf;
for(i=n;i!=1;i=tata[i])
fm=min(fm,cap[tata[i]][i]-flux[tata[i]][i]);
if(fm==0)continue;
for(i=n;i!=1;i=tata[i])
{
flux[tata[i]][i]+=fm;
flux[i][tata[i]]-=fm;
}
}
++pas;
}
}
void solve()
{ bool e,merge;
int i,x,y,p,u;
nod *aux;
for(i=1;i<=n;++i)viz[i]=0;
e=true; pas=1;
while(e)
{
e=false;
p=u=1; c[1]=1; viz[1]=pas; tata[1]=0; fol[1]=false;
while(p<=u)
{
if(c[p]!=n)
{
for(aux=G[c[p]];aux!=NULL;aux=aux->urm)
if(viz[aux->val]!=pas)
{
if(cap[c[p]][aux->val]>flux[c[p]][aux->val])
{
c[++u]=aux->val; tata[aux->val]=c[p]; viz[aux->val]=pas;
}
else if(cap[c[p]][aux->val]==flux[c[p]][aux->val]&&fol[c[p]]==false)
{
c[++u]=aux->val; tata[aux->val]=c[p]; viz[aux->val]=pas; fol[aux->val]=true;
}
}
}
++p;
}
if(viz[n]==pas)
{
e=true;
for(aux=G[n];aux!=NULL;aux=aux->urm)
{
merge=false;
if(viz[aux->val]==pas)
{
if(cap[aux->val][n]>flux[aux->val][n]&&fol[n]==false)merge=true;
else if(cap[aux->val][n]==flux[aux->val][n]&&fol[n]==true)merge=true;
}
if(merge)
{
tata[n]=aux->val; x=y=0;
for(i=n;i!=1;i=tata[i])
if(fol[i]==true){ x=tata[i]; y=i; }
sol[++k]=muc[x][y]; flux[x][y]+=inf; flux[y][x]-=inf;
}
}
}
++pas;
}
}
void write()
{ int i;
freopen("critice.out","w",stdout); printf("%d\n",k);
for(i=1;i<=k;++i) printf("%d\n",sol[i]);
fclose(stdout);
}
int main()
{
read();
maxflow();
k=0;
solve();
sort(sol+1,sol+k+1);
write();
return 0;
}