【网络流24题】餐巾计划（最小费用流）

最新推荐文章于 2017-08-28 19:22:57 发布

Etta19

最新推荐文章于 2017-08-28 19:22:57 发布

阅读量297

点赞数

分类专栏：网络流费用流

本文链接：https://blog.csdn.net/Etta19/article/details/77622327

版权

网络流同时被 2 个专栏收录

20 篇文章 0 订阅

订阅专栏

费用流

7 篇文章 0 订阅

订阅专栏

本文探讨了餐巾计划问题，通过引入最小费用最大流的概念来解决。文章设置了源汇点S和T，并对每天的餐巾需求进行建模。利用SPFA算法寻找最小费用的增广路径，以实现成本最低的餐巾分配方案。

摘要由CSDN通过智能技术生成

传送门

餐巾计划

I think

最小费用最大流。增设源汇点S，T，每i天分两个点xi,yi，分别表示每天应有餐巾与最终解决的餐巾数。记 $(x,y,u,v)$ 表示x连向y点容量为u费用为v点的边，连接以下边 $(S,x_i,a_i,0),(y_i,T,a_i,0),(x_i,x_{i+1},Inf,0),(x_i,y_i+qk,Inf,qc),(x_i,y_i+lk,Inf,lc);$ 最后用SPFA+增广路思想求出最小费用即可。
另算法函数名很无赖地用了SPFA。

Code

#include<cstdio>
#include<queue>
using namespace std;

const int sm = 2200;
const int sn = 12200;
const int Inf = 0x3f3f3f3f;

int N,nw,qk,qc,lw,lc,Ans,tot=1,S,T;
int to[sn],nxt[sn],hd[sm],c[sn],f[sn];
int pre[sm],vis[sm],cst[sm];

int Min(int x,int y) { return x<y?x:y; }

void Add(int u,int v,int w,int val) {
    to[++tot]=v,nxt[tot]=hd[u],hd[u]=tot,c[tot]=w,f[tot]=val;
    to[++tot]=u,nxt[tot]=hd[v],hd[v]=tot,c[tot]=0,f[tot]=-val;
}

void SPFA() {
    int t,df; queue<int>q;
    while(1) {
        for(int i=1;i<=T;++i) 
            pre[i]=vis[i]=0,cst[i]=Inf;
        cst[S]=0,vis[S]=1,q.push(S);
        while(!q.empty()) {
            t=q.front(),q.pop();
            vis[t]=0;
            for(int i=hd[t];i;i=nxt[i])
                if(cst[to[i]]>cst[t]+f[i]&&c[i]>0) {
                    cst[to[i]]=cst[t]+f[i];
                    pre[to[i]]=i;
                    if(!vis[to[i]]) {
                        vis[to[i]]=1;
                        q.push(to[i]);
                    }
                }
        }
        if(cst[T]==Inf)break;
        df=Inf;
        for(int i=T;i!=S;i=to[pre[i]^1]) 
            df=Min(df,c[pre[i]]);
        for(int i=T;i!=S;i=to[pre[i]^1]) 
            c[pre[i]]-=df,c[pre[i]^1]+=df;
        Ans+=df*cst[T];
    }
    printf("%d\n",Ans); 
}

int main() {
    int u;
    scanf("%d%d%d%d%d%d",&N,&nw,&qk,&qc,&lw,&lc);
    S=N<<1|1,T=S+1;
    for(int i=1;i<=N;++i) {
        scanf("%d",&u);
        Add(S,i,u,0),Add(i+N,T,u,0);
        Add(S,i+N,Inf,nw);
        if(i+1<=N) Add(i,i+1,Inf,0);
        if(i+qk<=N) Add(i,i+N+qk,Inf,qc);
        if(i+lw<=N) Add(i,i+N+lw,Inf,lc);
    }
    SPFA();
    return 0;
}