P3381 【模板】最小费用最大流
作者：互联网
题目描述
如题，给出一个网络图，以及其源点和汇点，每条边已知其最大流量和单位流量费用，求出其网络最大流和在最大流情况下的最小费用。
输入输出格式
输入格式：第一行包含四个正整数N、M、S、T，分别表示点的个数、有向边的个数、源点序号、汇点序号。
接下来M行每行包含四个正整数ui、vi、wi、fi，表示第i条有向边从ui出发，到达vi，边权为wi（即该边最大流量为wi），单位流量的费用为fi。
输出格式：一行，包含两个整数，依次为最大流量和在最大流量情况下的最小费用。
输入输出样例
输入样例#1：4 5 4 3 4 2 30 2 4 3 20 3 2 3 20 1 2 1 30 9 1 3 40 5输出样例#1：
50 280
说明
时空限制：1000ms,128M
（BYX：最后两个点改成了1200ms）
数据规模：
对于30%的数据：N<=10，M<=10
对于70%的数据：N<=1000，M<=1000
对于100%的数据：N<=5000，M<=50000
样例说明：
如图，最优方案如下：
第一条流为4-->3，流量为20，费用为3*20=60。
第二条流为4-->2-->3，流量为20，费用为（2+1）*20=60。
第三条流为4-->2-->1-->3，流量为10，费用为（2+9+5）*10=160。
故最大流量为50，在此状况下最小费用为60+60+160=280。
故输出50 280。
Solution:
费用流的模板，关于费用流该怎么算，大家可以去看下别人的博客，我这里也稍微讲一下：
一、求费用流的方法有很多，一般有两条途径：
第一种：先求出最大流，在保证最大流量不变的情况下，调整边的流量使得总费用减少，直到费用无法降低停止。显然，只要网络图中没有负圈（可以画图或想象一下），则当前的流就是最小费用流了。
第二种：每次寻找一条S到T的增广路，这条增广路必须是所有增广路中费用最小的一条，然后不停迭代知道无法增广。显然，此时的流为最小费用流。
二、具体实现：
由于上述的第二种方法是在不停增广，类似于求最大流算法，于是通常我们照此思路求最小费用流，即可转化为求S到T的最短路径问题。
1、费用网络只是在原网络中附加了边的费用cost(e)，那么建图时一切照常，只不过对于原有边<u,v>附加费用cost(e)，而它的反向边附加费用-cost(e)。
2、建图后，开始增广，注意到增广时只是累加边的费用而边的大小只要大于0就行了，但是因为边权有负值（反向边），所以选用SPFA求最短增广路。
3、因为负边会降低求最短路的效率，所以在每次建立增流网路求得最短路径后，将网络的权值进行一次修正，式再建的增流网络不含负边（要保证最短路径不会因此改变）。修改方法如下：当流量为0，首次建增流网络求最短路时，因为无负边权，所以可以采用标号法。为了使后面建立增流网络时不会出现负边，于是将网络中有流（flow>0）的边权w修改为0。
这样每次增流网络上求得最短路径后，新的边权 W"(u,v)=dis(u)-dis(v)+W(u,v) (dis(u)和dis(v)为新图中x到y的最短路径时u和v的标号值)。首次求最短路时，若(u,v)是增流路径上的边，则由最短路算法必有dis(v)=dis(u)+W'(u,v)=dis(u)+W(u,v) 将这个式子带入加黑式子有W"(u,v)=0。
若(u,v)不是增流路径上的边，则dis(v)<=dis(u)+W(u,v) 再代入加黑式子，则W"(u,v)>=0。显然，首次修正w后，对于网络中的任意边w都大于等于0，并且有流量的边（增流路径上的边）的w一定等于0。之后的迭代时，若f(u,v)>0，增流网络建立（v,u）边，新的边权W'(v,u)=-W(u,v)=0，所以不会出现负边了。那么每次迭代都用最开始的加黑的式子去修改边权w，显然对于每一条x到y的路径，它的路径长度都同样增加dis(x)-dis(y)。所以，x到y的最短路径就不会因为对w的修改而发生改变了。
这里有字
我这里用的是用双端队列优化过的SPFA的代码：
1 // luogu-judger-enable-o2
2 #include<bits/stdc++.h>
3 #define il inline
4 using namespace std;
5 bool vis[5005];int dist[5005];
6 int n,m,s,t,ans=0;
7 int nedge=-1,p[100002],c[100002],cc[100002],nex[100002],head[5005];
8 il int gi()
9 {
10 int a=0;char x=getchar();bool f=0;
11 while((x<'0'||x>'9')&&x!='-')x=getchar();
12 if(x=='-')x=getchar(),f=1;
13 while(x>='0'&&x<='9')a=a*10+x-48,x=getchar();
14 return f?-a:a;
15 }
16 il void addedge(int x,int y,int z,int zz){
17 p[++nedge]=y;c[nedge]=z;cc[nedge]=zz;nex[nedge]=head[x];head[x]=nedge;
18 }
19 il bool spfa(int s,int t){
20 memset(vis,0,sizeof vis);
21 for(int i=0;i<=n;i++)dist[i]=1e9;dist[t]=0;vis[t]=1;
22 deque<int>q;q.push_back(t);
23 while(!q.empty()){
24 int now=q.front();q.pop_front();
25 for(int k=head[now];k>-1;k=nex[k])if(c[k^1]&&dist[p[k]]>dist[now]-cc[k]){
26 dist[p[k]]=dist[now]-cc[k];
27 if(!vis[p[k]]){
28 vis[p[k]]=1;
29 if(!q.empty()&&dist[p[k]]<dist[q.front()])q.push_front(p[k]);else q.push_back(p[k]);
30 }
31 }
32 vis[now]=0;
33 }
34 return dist[s]<1e9;
35 }
36 il int dfs(int x,int low){
37 if(x==t){vis[t]=1;return low;}
38 int used=0,a;vis[x]=1;
39 for(int k=head[x];k>-1;k=nex[k])if(!vis[p[k]]&&c[k]&&dist[x]-cc[k]==dist[p[k]]){
40 a=dfs(p[k],min(c[k],low-used));
41 if(a)ans+=a*cc[k],c[k]-=a,c[k^1]+=a,used+=a;
42 if(used==low)break;
43 }
44 return used;
45 }
46 il int costflow(){
47 int flow=0;
48 while(spfa(s,t)){
49 vis[t]=1;
50 while(vis[t]){
51 memset(vis,0,sizeof vis);
52 flow+=dfs(s,1e9);
53 }
54 }
55 return flow;
56 }
57 int main()
58 {
59 memset(nex,-1,sizeof nex);memset(head,-1,sizeof head);
60 n=gi(),m=gi(),s=gi(),t=gi();
61 int x,y,z,zz;
62 for(int i=1;i<=m;i++){
63 x=gi(),y=gi(),z=gi(),zz=gi();
64 addedge(x,y,z,zz);addedge(y,x,0,-zz);
65 }
66 printf("%d ",costflow());printf("%d",ans);
67 return 0;
68 }
提供一种不会被普通SPFA卡极限数据或者zkw卡极限数据的代码（速度不是特别快，关键是稳）：
1 // luogu-judger-enable-o2
2 #include <bits/stdc++.h>
3 #define il inline
4 #define debug printf("%d %s\n",__LINE__,__FUNCTION__)
5 using namespace std;
6 const int N=100005,inf=0x3f3f3f3f;
7 int n,m,s,t,cnt=1,h[5005],dis[5005],ans,fee,flow;
8 bool inq[5005],vis[5005];
9 struct edge{
10 int to,net,w,c;
11 }e[N];
12 queue<int> q;
13 il int gi()
14 {
15 int a=0;char x=getchar();bool f=0;
16 while((x<'0'||x>'9')&&x!='-')x=getchar();
17 if(x=='-')x=getchar(),f=1;
18 while(x>='0'&&x<='9')a=a*10+x-48,x=getchar();
19 return f?-a:a;
20 }
21 il void add(int u,int v,int w,int c)
22 {
23 e[++cnt].to=v,e[cnt].net=h[u],e[cnt].c=c,e[cnt].w=w,h[u]=cnt;
24 e[++cnt].to=u,e[cnt].net=h[v],e[cnt].c=-c,e[cnt].w=0,h[v]=cnt;
25 }
26 il bool bfs(int s,int t)
27 {
28 memset(dis,0x3f,sizeof(dis));
29 memset(inq,0,sizeof(inq));
30 int x;
31 dis[t]=0;
32 q.push(t);
33 inq[t]=1;
34 while(!q.empty())
35 {
36 x=q.front();
37 q.pop();
38 inq[x]=0;
39 for(int i=h[x];i;i=e[i].net)
40 if(e[i^1].w&&dis[e[i].to]>dis[x]-e[i].c)
41 {
42 dis[e[i].to]=dis[x]-e[i].c;
43 if(inq[e[i].to])continue;
44 inq[e[i].to]=1;
45 q.push(e[i].to);
46 }
47 }
48 return dis[s]<inf;
49 }
50 il int dfs(int u,int flow)
51 {
52 vis[u]=1;
53 if(u==t||flow==0)
54 return flow;
55 int used=0,w;
56 for(int i=h[u];i;i=e[i].net)
57 if(!vis[e[i].to]&&e[i].w&&dis[e[i].to]==dis[u]-e[i].c)
58 {
59 w=dfs(e[i].to,min(e[i].w,flow-used));
60 e[i].w-=w;
61 e[i^1].w+=w;
62 used+=w;
63 if(used==flow)return used;
64 }
65 return used;
66 }
67 int main()
68 {
69 n=gi(),m=gi(),s=gi(),t=gi();
70 int u,v,w,c;
71 for(int i=1;i<=m;i++)u=gi(),v=gi(),w=gi(),c=gi(),add(u,v,w,c);
72 while(bfs(s,t))
73 {
74 vis[t]=1;
75 while(vis[t])
76 {
77 memset(vis,0,sizeof(vis));
78 flow=dfs(s,inf);
79 ans+=flow;
80 fee+=flow*dis[s];
81 }
82 }
83 printf("%d %d",ans,fee);
84 return 0;
85 }
标签：费用,dist,int,最小,流量,vis,P3381,模板,dis 来源： https://blog.51cto.com/u_15180869/2799032