您的位置:首页 > 产品设计 > UI/UE

zoj3792 Romantic Value 最大流<->最小割

2015-09-06 21:56 441 查看
一道以最大流求最小割的题目,最大流运用模板能够求出来,可是在求最少割边的数目却不会了,百度到两种方法:

1.每条边的权值乘以1000然后+1,求出最大流/1000就是原本的最大流,最少割边就是最大流%1000

2.求出最大流后,遍历每一条边,满流量就设置为1,否则为0,重新求一遍最大流,就是最少的割边

比赛的时候一直死扣第一种方法,结果怎么都不对,比赛结束后去找了代码,用的是第二种方法,然后把这份AC代码改成第一种方法结果却是WA,现在还是不知道为什么,有空一定要弄明白原因。

运用第二种方法的AC代码

<span style="font-family:KaiTi_GB2312;font-size:18px;">#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#define inf 0x3f3f3f3f
using namespace std;
struct node
{
int v,w,next;
} edge[500000];
int head[5000],cnt,n,m,l[5000],s,t,p,q,T;
void add(int u,int v,int w) //加边
{
edge[cnt].v=v;
edge[cnt].w=w;
edge[cnt].next=head[u];
head[u]=cnt++;
edge[cnt].v=u;
edge[cnt].w=0;
edge[cnt].next=head[v];
head[v]=cnt++;
}
int bfs()
{
queue<int> Q;
Q.push(s);
memset(l,-1,sizeof(l));
l[s]=0;
while(!Q.empty())
{
int u=Q.front();
Q.pop();
for(int i=head[u]; i!=-1; i=edge[i].next)
{
int v=edge[i].v;
if(l[v]==-1&&edge[i].w)
{
l[v]=l[u]+1;
Q.push(v);
}
}
}
return l[t]>0;
}
int dfs(int x,int f) //求源点到汇点的最大流
{
if(x==t)return f;
int a;
for(int i=head[x]; i!=-1; i=edge[i].next)
{
int v=edge[i].v;
if(edge[i].w&&l[v]==l[x]+1&&(a=dfs(v,min(edge[i].w,f))))
{
edge[i].w-=a;
edge[i^1].w+=a;
return a;
}
}
l[x]=-1;
return 0;
}
int main()
{
scanf("%d",&T);
while(T--)
{
int a,b,v,ans;
scanf("%d%d%d%d",&n, &m, &p, &q);
memset(head, -1, sizeof(head));
cnt = 0;
ans = 0;
s = p;
t = q;
int sum = 0;
for(int i = 0; i < m; i++)
{
scanf("%d%d%d",&a, &b, &v);
add(a,b,v);
add(b,a,v);
sum += v;
}
if(!bfs())
{
printf("Inf\n");
continue;
}
while(bfs())
while(a = dfs(p,inf))
ans+=a;
for(int i=0; i<cnt; i+=2)
{
if(edge[i].w==0)
edge[i].w=1;
else edge[i].w=inf;
edge[i^1].w=0;
}
int nnt=0;
while(bfs())
while(a=dfs(s,inf))
nnt+=a;
/*
//printf("ans == %d\n",ans);
int flow = sum - ans/1000;
int edge_num = ans%1000;
if(edge_num == 0|| ans == 0||flow == 0)
printf("Inf\n");
else*/
printf("%.2lf\n", (double)((sum - ans)*1.0)/nnt);
}
}
</span>
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签:  acm 最大流