1 typedef long long ll;
2 const ll inf=0x3f3f3f3f;
3 struct edge{
4 int x;
5 int y;
6 ll v;//流量
7 ll c;//费用
8 int nex;
9 }e[1005];
10 ll flow,res;
11 int cnt,head[505],n;
12 void adde(int x1,int x2,ll v1,ll v2){
13 e[cnt].c=v2;
14 e[cnt].v=v1;
15 e[cnt].x=x1;
16 e[cnt].y=x2;
17 e[cnt].nex=head[x1];
18 head[x1]=cnt++;
19 }
20 int dist[505],prevv[505],preve[1005]; // 最短路的前驱节点 和 对应的边
21 inline void min_cost_flow(int s/*源点*/,int t/*汇点*/)
22 {
23 ll f=inf;
24 while(f > 0)
25 {
26 memset(dist,inf,sizeof dist);
27 dist[s] = 0;
28 bool update = true;
29 while(update)
30 {
31 update = false;
32 for(int i=0;i<=n+1;i++)
33 {
34 if(dist[i] == inf) continue;
35 for(int j=head[i];j!=-1; j=e[j].nex)
36 {
37 if(e[j].v > 0 && dist[e[j].y] > dist[i] + e[j].c)
38 {
39 dist[e[j].y] = dist[i] + e[j].c;
40 prevv[e[j].y] = i;
41 preve[e[j].y] = j;
42 update = true;
43 }
44 }
45 }
46 }
47 // 此时的 INF 表示再也没有 增广路径,于是就是最后的答案了!
48 if(dist[t] == inf) break;
49 ll d = f;
50 for(int v0 = t; v0 != s; v0 = prevv[v0])
51 {d = min(d,e[preve[v0]].v);}
52 // d 就是 增广路的 流量!
53 f -= d;
54 flow += d;
55 res += d * dist[t];
56 for(int v0=t;v0!=s;v0=prevv[v0])
57 {
58 e[preve[v0]].v -= d;
59 e[preve[v0]^1].v += d; //反向边流量加相应的值
60 }
61 }
62 }