http://acm.hdu.edu.cn/showproblem.php?pid=6214
建边的时候每条边权 w = w * (E + 1) + 1;
这样得到最大流 maxflow / (E + 1) ,最少割边数 maxflow % (E + 1)
道理很简单,如果原先两类割边都是最小割,那么求出的最大流相等
但边权变换后只有边数小的才是最小割了
大佬的解释:
乘(E+1)是为了保证边数叠加后依然是余数,不至于影响求最小割的结果
因为假设最小割=k,那么现在新图的最小割为k*(E+1)+p,p为割的边数,本质上是,原来你割一条边,需要代价,
由于你要求边数最小 所以你多割一条边,就多一的代价,但是这个代价不足以影响到原来的代价。
原来割一条边,代价xi,现在割一条边,代价xi*A+1,只要让A>m+1,m为边数,即使割了所有的边,自己加上去的代价也就m。
用两遍dinic 一直错,不过好像先跑完,再用残量网络再跑一遍对了。。
#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define INF 1<<30
#define maxn 22100
#define maxm 500000
using namespace std;
int v[maxm],Next[maxm],w[maxm];
int ww[maxn];
int first[maxn],d[maxn],q[maxn];
int e,S,T;
int n;
void init(){
e = 0;
memset(first,-1,sizeof(first));
}
void add_edge(int a,int b,int c){
v[e] = b;Next[e] = first[a];w[e] = c*(n+1)+1;first[a] = e++;
v[e] = a;Next[e] = first[b];w[e] = 0;first[b] = e++;
ww[e]=w[e];
}
int bfs(){
int rear = 0;
memset(d,-1,sizeof(d));
d[S] = 0;q[rear++] = S;
for(int i = 0;i < rear;i++){
for(int j = first[q[i]];j != -1;j = Next[j])
if(w[j] && d[v[j]] == -1){
d[v[j]] = d[q[i]] + 1;
q[rear++] = v[j];
if(v[j] == T) return 1;
}
}
return 0;
}
int dfs(int cur,int maxflow){
if(cur == T) return maxflow;
for(int i = first[cur];i != -1;i = Next[i]){
if(w[i] && d[v[i]] == d[cur] + 1)
if(int t = dfs(v[i],min(maxflow,w[i]))){
w[i] -= t;w[i^1] += t;
return t;
}
}
return 0;
}
int dinic(){
int ans = 0;
while(bfs()){
ans += dfs(S,INF);
}
return ans;
}
int main()
{
int m,ncase,cas = 1;
scanf("%d",&ncase);
while(ncase--)
{
init();
//scanf("%d%d",&S,&T);
scanf("%d%d",&n,&m);
scanf("%d%d",&S,&T);
//S = 0,T = n-1;
for(int i = 0;i < m;i++){
int a,b,c,d;
scanf("%d%d%d",&a,&b,&c);
add_edge(a,b,c);
//add_edge(b,a,c);
}
//dinic();
printf("%d\n",dinic()%(n+1));
}
return 0;
}