「网络流 24 题」餐巾计划 费用流 新流新用vs旧流新用

题目链接：https://loj.ac/problem/6008

 

        中文的题面，所以我就不解释了，之前做了一道飞行员的题目和这道几乎是一个意思。我就自作主张的起了“新流新用vs旧流新用”这个标题，为什么呢，解释一下这道题的建图。

 

        这个题把每一天的餐巾情况进行了拆点，点i代表的是这天因为用过了之后留下来要洗的餐巾，点i’代表的是当天可用的干净餐巾的数量。

 

        假设我们不考虑有旧流补入的情况（不考虑可以洗餐巾的情况）

        先将源点到1'（我的代码中是2），连一条容量为inf费用为P的边，代表可以用P一条的价格买入新的餐巾。

        将i'->i+1'连容量为inf费用为0的边，代表当天如果有没用用完的干净的餐巾，那么可以留到下一天（每天都要买新餐巾的情况直接用源点到1'代表了，因为如果后面要买的话，直接从前面把流扩大就好（即第一天会买好所有需要的），就不用每天都从源点连边啦）。

        将i'向汇点连一条流量为need[i]，费用为0的边，表示每天都要提供这么多餐巾，直接把这些流给收掉。

 

        然后我们再考虑有旧流补入的情况，那么我们就要用到点i了；

        因为每天都会有need[i]那么多的旧餐巾产生，所以直接从源点向点i连边，费用为0。

        重点来了！然后就是要连边i->i+M,容量为inf，费用为F，i->i+N,容量为inf，费用为S。为什么呢，因为从第i天可以直接向i+M天通过洗涤得到干净的餐巾，即这些流可以直接变为新的流再往下跑。为什么是inf呢，其实只要大于need[i]就好了，因为在源点连边的时候就已经限制了流量。

---

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=6000;
const int maxm=600000;
const int inf=0x3f3f3f3f;
int dis[maxn];
int vis[maxn],pre[maxn];
int head[maxn],cnt;
int n,m,sp,tp;
struct node{
    int to,cap,cost,next;
}e[maxm];
void add(int from,int to,int cap,int cost){
    e[cnt].to=to; e[cnt].cap=cap;
    e[cnt].cost=cost; e[cnt].next=head[from];
    head[from]=cnt++;

    e[cnt].to=from; e[cnt].cap=0;
    e[cnt].cost=-cost; e[cnt].next=head[to];
    head[to]=cnt++;
}
bool spfa(int s,int t,int &flow,int &cost){
    queue<int> q;
    memset(dis,inf,sizeof(dis));
    memset(vis,0,sizeof(vis));
    memset(pre,-1,sizeof(pre));
    dis[s]=0;  q.push(s);
    vis[s]=1;
    int d=inf;
    while(!q.empty()){
        int u=q.front();
        q.pop();
        vis[u]=0;
        for(int i=head[u];~i;i=e[i].next){
            int v=e[i].to;
            if(e[i].cap>0&&dis[v]>dis[u]+e[i].cost){
                dis[v]=dis[u]+e[i].cost;
                pre[v]=i;
                if(!vis[v]){
                    vis[v]=1;
                    q.push(v);
                }
            }
        }
    }
    if(dis[t]==inf){
        return false;
    }
    for(int i=pre[t];~i;i=pre[e[i^1].to]){
        d=min(d,e[i].cap);
    }
    for(int i=pre[t];~i;i=pre[e[i^1].to]){
        e[i].cap-=d;
        e[i^1].cap+=d;
        cost+=e[i].cost*d;
    }
    flow+=d;
    return true;
}
int mcmf(int s,int t){
    int flow=0,cost=0;
    while(spfa(s,t,flow,cost)){
        //cout<<flow<<" "<<cost<<endl;
    }
    return cost;
}
int need[maxn];
int gainp(int x,int y){
    return (x-1)*2+y;
}
int main(){
    int P,M,F,N,S;
    memset(head,-1,sizeof(head));

    scanf("%d%d%d%d%d%d",&n,&P,&M,&F,&N,&S);
    cnt=0,sp=0,tp=2*n+1;
    for(int i=1;i<=n;i++)
        scanf("%d",&need[i]);
    add(sp,2,inf,P);
    for(int i=1;i<=n;i++){
        add(sp,gainp(i,1),need[i],0);
        add(gainp(i,2),tp,need[i],0);
        if(i!=n) add(gainp(i,2),gainp(i+1,2),inf,0);
        if(i+M<=n)add(gainp(i,1),gainp(i+M,2),inf,F);
        if(i+N<=n) add(gainp(i,1),gainp(i+N,2),inf,S);
    }
    printf("%d\n",mcmf(sp,tp));
    return 0;
}