「LuoguP3381」【模板】最小费用最大流

Description

---

如题，给出一个网络图，以及其源点和汇点，每条边已知其最大流量和单位流量费用，求出其网络最大流和在最大流情况下的最小费用。

Input

---

第一行包含四个正整数N、M、S、T，分别表示点的个数、有向边的个数、源点序号、汇点序号。

接下来M行每行包含四个正整数ui、vi、wi、fi，表示第i条有向边从ui出发，到达vi，边权为wi（即该边最大流量为wi），单位流量的费用为fi。

Output

---

一行，包含两个整数，依次为最大流量和在最大流量情况下的最小费用。

Sample Input

---

4 5 4 3
4 2 30 2
4 3 20 3
2 3 20 1
2 1 30 9
1 3 40 5

Sample Output

---

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。

题解

---

学费用流的模板题 选了zkw

对着hzwer的模板敲了一遍网络流&费用流模板by hzwer

其实跟dinic差不多 就是把bfs换成spfa 以单位费用为dis从源点跑 返回dis[t]是否能到达

以及不知道为什么spfa要反着搜 正搜也A了2333

 #include<queue>
 #include<cstdio>
 #include<cstring>
 #include<iostream>
 using namespace std;
 #define R register
 const int INF=;
 inline int read()
 {
     char ch=getchar();
     int s=,x=;
     while(ch<''||ch>''){if(ch=='-')s=-;ch=getchar();}
     while(ch>=''&&ch<=''){x=x*+ch-;ch=getchar();}
     return s*x;
 }
 int n,m,s,t;
 struct emm{
     int e,f,v,c;
 }a[];
 int h[];
 int tot=;
 void con(int x,int y,int vv,int cc)
 {
     a[++tot].f=h[x];
     h[x]=tot;
     a[tot].e=y;
     a[tot].v=vv;
     a[tot].c=cc;
     a[++tot].f=h[y];
     h[y]=tot;
     a[tot].e=x;
     a[tot].c=-cc;
     return;
 }
 void scan()
 {
     n=read(),m=read(),s=read(),t=read();
     for(R int i=;i<=m;++i)
     {
         int x=read(),y=read(),v=read(),c=read();
         con(x,y,v,c);
     }
     return;
 }
 queue<int>q;
 bool sf[];
 int d[];
 bool spfa()
 {
     memset(sf,,sizeof(sf));
     //memset(d,127,sizeof(d));
     for(int i=;i<=n;++i)d[i]=INF;
     sf[s]=;d[s]=;q.push(s);
     while(!q.empty())
     {
         int x=q.front();q.pop();
         for(int i=h[x];i;i=a[i].f)
         if(d[x]+a[i].c<d[a[i].e]&&a[i].v)
         {
             d[a[i].e]=d[x]+a[i].c;
             if(!sf[a[i].e])
             {
                 sf[a[i].e]=;
                 q.push(a[i].e);
             }
         }
         sf[x]=;
     }
     //for(int i=1;i<=4;++i)
     //cout<<d[i]<<" ";
     //cout<<endl;
     return d[t]<INF;
 }
 int ans=;
 int dfs(int x,int al)
 {
     sf[x]=;
     if(x==t||!al)return al;
     int fl=;
     for(int i=h[x];i;i=a[i].f)
     if(d[a[i].e]==d[x]+a[i].c&&a[i].v&&!sf[a[i].e])
     {
         int f=dfs(a[i].e,min(al,a[i].v));
         fl+=f;
         al-=f;
         ans+=f*a[i].c;
         a[i].v-=f;
         a[i^].v+=f;
         if(!al)break;
     }
     if(!fl)d[x]=-;
     return fl;
 }
 int now=;
 void zkw()
 {
     while(spfa())
     {
         sf[t]=;
         while(sf[t])
         {
             memset(sf,,sizeof(sf));
             now+=dfs(s,INF);
         }
     }
     return;
 }
 int main()
 {
     scan();
     zkw();
     cout<<now<<" "<<ans;
     return ;
 }

csdn好丑啊 可是好看的lofter根本不是用来写博客的吧？！

UPD

现在是十一月了　我写了这么多次费用流都是正向跑的spfa　没出过一点问题qwq