链接:https://www.luogu.org/problem/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。
思路:费用流模板题。
AC代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn=5050,maxm=50050;
const int INF=0x3f3f3f3f;
int n,m,s,t,maxflow,mincost;
int N;
int head[maxn],tot;
struct node
{
int to,next,w,f;
//w是容量,f是费用。
}edge[maxm*2];
inline void add_edge(int u,int v,int w,int f)
{
edge[tot]=(node){v,head[u],w,f};
head[u]=tot++;
edge[tot]=(node){u,head[v],0,-f};
head[v]=tot++;
}
inline void edge_init()
{
tot=0;N=n;
memset(head,-1,sizeof(head));
}
int dis[maxn],pre[maxn],minw[maxn];
bool vis[maxn];
bool spfa()
{
queue<int> que;
for(int i=0;i<=N;++i)
{
dis[i]=INF;pre[i]=-1;vis[i]=0;minw[i]=0;
}
que.push(s),dis[s]=0,vis[s]=1,minw[s]=INF;
while(!que.empty())
{
int u=que.front();que.pop();vis[u]=0;
for(int i=head[u];i!=-1;i=edge[i].next)
{
int v=edge[i].to;
if(edge[i].w==0) continue;
if(dis[u]+edge[i].f<dis[v])
{
dis[v]=dis[u]+edge[i].f,pre[v]=i,minw[v]=min(minw[u],edge[i].w);
if(!vis[v]) que.push(v),vis[v]=1;
}
}
}
return dis[t]!=INF;
}
void update()
{
int u=t;
maxflow+=minw[t],mincost+=dis[t]*minw[t];
while(u!=s) edge[pre[u]].w-=minw[t],edge[pre[u]^1].w+=minw[t],u=edge[pre[u]^1].to;
}
void solve()
{
while(spfa()) update();
}
int main()
{
scanf("%d%d%d%d",&n,&m,&s,&t);
edge_init();
for(int i=1;i<=m;i++)
{
int a,b,c,d;
scanf("%d%d%d%d",&a,&b,&c,&d);
add_edge(a,b,c,d);
}
solve();
printf("%d %d\n",maxflow,mincost);
return 0;
}