ZOJ 3792 最小割+求割边

题目链接：http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5300

题意：有个人特烦他某个亲戚。。称他不在的时候，想炸掉一些路，让他过不来，他对每条路都有喜爱程度，要 求在损失喜爱程度相同的情况下炸掉最少路的方案，然后根据这个方案算一个数= 。=

思路：损失最少就是求最小割，套上最大流模板就可以了，然后把满流的边流量设为1，不满的设为inf，再跑一遍最大流就是割边的数目（  这里不是十分理解  ，求讲解。。），我的理解是不考虑哪些inf的边，因为割边一定在那些新设为1的边中，跑一下最大流之后就是求出来源点和汇点之间有几条独立的路（因为流量为1嘛），然后割掉每条路中的一条边就好了。。。

PS：  ZOJ  这道题咯坑爹啊。。应该是PE的给了我一个WA，各种找都找不到啊。。。

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int maxn = 210;
const int maxm = 5000;
const int INF = 0x7FFFFFF;
struct Side{
	int to,next,c;
}side[maxm];
int top,node[maxn];
void add_side(int u,int v,int c,int rc){
	side[top]=(Side){v,node[u],c};
	node[u]=top++;
	side[top]=(Side){u,node[v],rc};
	node[v]=top++;
}
int start,end,cnt,dis[maxn],gap[maxn];
int get_flow(int u,int flow){
	//printf("%d %d\n",u,flow);
	if(u==end)return flow;
	int ans=0;
	for(int i=node[u];i!=-1;i=side[i].next){
		int v=side[i].to,c=side[i].c;
		if(dis[u]>dis[v]&&c){
			int f=get_flow(v,min(flow-ans,c));
			ans+=f;
			side[i].c-=f;
			side[i^1].c+=f;
			if(ans==flow)return ans;
		}
	}
	if(!(--gap[dis[u]]))dis[start]=cnt+2;
	gap[++dis[u]]++;
	return ans;
}
int main(){//注意maxn maxm的修改
	int n,m;
	int cas;
	scanf("%d",&cas);
	while(cas--){
        scanf("%d%d%d%d",&n,&m,&start,&end);
		top=0;
		memset(node,-1,sizeof(node));
		memset(gap,0,sizeof(gap));
		memset(dis,0,sizeof(dis));
        int all = 0;
		while(m--){
			int u,v,w;
			scanf("%d%d%d",&u,&v,&w);
			add_side(u,v,w,0);
			add_side(v,u,w,0);
			all+=w;
		}
		int x=0;
		//start=1;//源点
		//end=n;//汇点
		cnt=n;//点的总个数
		gap[0]=cnt;
		while(dis[start]<cnt){
            x+=get_flow(start,INF);
		}
        if(x==0){printf("Inf\n");continue;}
        //all-=ans;
        //ans=0;
        memset(gap,0,sizeof(gap));
		memset(dis,0,sizeof(dis));
        //cnt = n;
        gap[0]=cnt;
		for(int i=0;i<top;i+=2){
            if(side[i].c==0)side[i].c=1;
            else side[i].c=INF;
            side[i^1].c=0;
		}
		int tot=0;
		while(dis[start]<cnt){
            tot+=get_flow(start,INF);
		}
		double ans = (all-x)*1.0/tot;
		printf("%.2f",ans);
		//printf("%d\n",ans);//最大流/最小割
	}
	return 0;
}