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>
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网址
- 1272 小希的迷宫
- 1272 小希的迷宫
- hdu 1250 大数相加并用数组储存
- 矩阵的乘法操作
- 蚂蚁爬行问题
- 蚂蚁爬行问题
- 求两个数的最大公约数【ACM基础题】
- 打印出二进制中所有1的位置
- 杭电题目---一只小蜜蜂
- HDOJ 1002 A + B Problem II (Big Numbers Addition)
- POJ ACM 1002
- POJ 2635 The Embarrassed Cryptographe
- POJ 3292 Semi-prime H-numbers
- POJ 2773 HAPPY 2006
- POJ 3090 Visible Lattice Points
- HDU 1568
- HDU1290