Luogu 3381（最小费用流）

来源：互联网发布：美股行情交易软件编辑：程序博客网时间：2024/06/05 02:44

传送门

模板题。最小费用最大流，要求同时输出最大流以及最小费用，复杂度有点玄学，每增广一次，下一次SPFA的复杂度又会降低。所以这种题暂时只能这么处理：看出来有后效性的相互约束关系的就考虑用费用流解决。

#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>#include<queue>using namespace std;const int N=5004,M=50004,INF=0x3f3f3f3f;int n,m,source,sink;int head[N],etot=0,dis[N],pre[N];bool inq[N];struct EDGE {int u,v,nxt,c,r;}e[M<<1];inline int read() {int x=0;char c=getchar();while (c<'0'||c>'9') c=getchar();while (c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-'0',c=getchar();return x;}inline void adde(int u,int v,int r,int c) {e[etot].nxt=head[u],e[etot].u=u,e[etot].v=v,e[etot].r=r,e[etot].c=c,head[u]=etot++;e[etot].nxt=head[v],e[etot].u=v,e[etot].v=u,e[etot].r=0,e[etot].c=-c,head[v]=etot++;}inline bool SPFA() {queue<int > q;memset(dis,INF,sizeof(dis));memset(inq,false,sizeof(inq));pre[source]=-1,dis[source]=0,inq[source]=true;q.push(source);while (!q.empty()) {int p=q.front();q.pop();inq[p]=false;for (int i=head[p];~i;i=e[i].nxt) {int v=e[i].v;if (e[i].r>0&&dis[v]>dis[p]+e[i].c) {dis[v]=dis[p]+e[i].c;pre[v]=i;if (!inq[v]) q.push(v),inq[v]=true;}}}return dis[sink]^INF;}inline void MCMF() {int maxflow=0,mincost=0;while (SPFA()) {int temp=INF;for (int i=pre[sink];~i;i=pre[e[i].u]) temp=min(temp,e[i].r);maxflow+=temp,mincost+=temp*dis[sink];for (int i=pre[sink];~i;i=pre[e[i].u]) e[i].r-=temp,e[i^1].r+=temp;}printf("%d %d\n",maxflow,mincost);}int main() {//freopen("P3381.in","r",stdin);memset(head,-1,sizeof(head));n=read(),m=read(),source=read(),sink=read();for (register int i=0;i<m;++i) {int u=read(),v=read(),r=read(),c=read();adde(u,v,r,c);}MCMF();return 0;}

阅读全文

0 0