题目链接
http://acm.hdu.edu.cn/showproblem.php?pid=6214
题意
给定一有向图,求图中边数最少的最小割
思路
跑两次最大流:第一次正常建图,第二次在满流的边上跑。比赛时用的是这样方法,后来听大神说这种算法是错误的,因为最小割一定满流,满流的不一定是最小割。。
在网上看了第二种做法,在原图的基础上,每条边权都乘上一个较大的数E在加1,即
w = w*E+1 。这样在新图中,只有边权最少的最小割是最小割了,因为每条边都会贡献一个1.所以在最后求得的最大流%eE就是最小边数.
#include<cstdio>
#include<queue>
#include<iostream>
#include<vector>
#include<map>
#include<cstring>
#include<string>
#include<set>
#include<stack>
#include<algorithm>
#define cle(a) memset(a,0,sizeof(a))
#define inf(a) memset(a,0x3f,sizeof(a))
#define ll long long
#define Rep(i,a,n) for(int i=a;i<=n;i++)
using namespace std;
const int INF = ( 2e9 ) + 2;
const int maxn = 1010;
struct Dinic
{
struct edge
{
int from,to,c,f;
edge(int getu,int getv,int getc,int getf)
{
from=getu;
to=getv;
c=getc;
f=getf;
}
};
int s,t,n;
int d[maxn];
int cur[maxn];
vector<edge> e;
vector<int> g[maxn];
void addedge(int u,int v,int c)
{
e.push_back(edge(u,v,c,0));
e.push_back(edge(v,u,0,0));
int m=e.size();
g[u].push_back(m-2);
g[v].push_back(m-1);
}
void init()
{
e.clear();
for(int i=0; i<maxn; i++)
g[i].clear();
}
bool bfs()
{
queue<int> q;
q.push(s);
memset(d,0,sizeof(d));
d[s]=1;
while(!q.empty())
{
int u=q.front();
q.pop();
for(int i=0; i<g[u].size(); i++)
{
int v=e[g[u][i]].to;
int c=e[g[u][i]].c;
int f=e[g[u][i]].f;
if(!d[v]&&c-f>0)
{
d[v]=d[u]+1;
q.push(v);
}
}
}
return d[t];
}
int dfs(int u,int maxf,int t)
{
if(u==t)return maxf;
int ret=0;
for(int &i=cur[u]; i<g[u].size(); i++)
{
int c=e[g[u][i]].c-e[g[u][i]].f;
int v=e[g[u][i]].to;
int f;
if(d[u]+1==d[v]&&c)
{
f=dfs(v,min(maxf-ret,c),t);
e[g[u][i]].f+=f;
e[g[u][i]^1].f-=f;
ret+=f;
if(ret==maxf)return ret;
}
}
return ret;
}
int maxflow(int s,int t)
{
this->s=s;
this->t=t;
int flow=0;
while(bfs())
{
memset(cur,0,sizeof(cur));
int temp=dfs(s,INF,t);
flow+=temp;
}
return flow;
}
}dinic;
int main()
{
int n,m,w;
int T;
// freopen("in.txt","r",stdin);
// freopen("out1.txt","w",stdout);
scanf("%d",&T);
for(int cas=1;cas<=T;cas++)
{
dinic.init();
int s,t;
int u,v;
scanf("%d%d%d%d",&n,&m,&s,&t);
for(int i=0;i<m;i++)
{
scanf("%d%d%d",&u,&v,&w);
dinic.addedge(u,v,w*10000+1);
}
printf("%d\n",dinic.maxflow(s,t)%10000);
}
}