题目链接: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;
}