- 因为含有负权边,所以使用spfa进行增广。
- 指定流量的最小费用流可以初始化一个f,然后每次一直循环到f<= 0为止。
- 函数返回的是最大流,当然经过少量修改可以返回最小费用,利用最小流量乘以相应的费用即可。
- prevv记录父节点,preve记录当前节点对应父节点的第几条边。
const int INF = 0x3fffffff;
struct Edge {
int to, cap, cost, rev;
Edge(int too, int capp, int costt, int revv) : to(too), cap(capp), cost(costt), rev(revv) {}
};
vector<vector<Edge> > G;
int st, en;
int dist[30];
int prevv[30], preve[30];
void add_edge(int from, int to, int cap, int cost) {
G[from].push_back(Edge(to, cap, cost, G[to].size()));
G[to].push_back(Edge(from, 0, -cost, G[from].size()-1));
}
int min_max_flow(int n) {
int ans = 0;
while(true) {
fill(dist, dist+n, INF);
dist[st] = 0;
queue<int> q; q.push(st);
while(!q.empty()) {
int u = q.front(); q.pop();
for(int i = 0; i < G[u].size(); i++) {
int to = G[u][i].to;
if(G[u][i].cap > 0 && dist[u] + G[u][i].cost < dist[to]) {
dist[to] = dist[u] + G[u][i].cost;
prevv[to] = u; preve[to] = i;
q.push(to);
}
}
}
if(dist[en] == INF)
break;
else {
int f = INF;
for(int pr = en; pr != st; pr = prevv[pr])
f = min(f, G[prevv[pr]][preve[pr]].cap);
ans += f;
for(int pr = en; pr != st; pr = prevv[pr]) {
Edge &e = G[prevv[pr]][preve[pr]];
e.cap -= f;
G[pr][e.rev].cap += f;
}
}
}
return ans;
}