[网络流24题]餐巾计划问题

传送门

一道建模稍微有难度的网络流，首先你需要考虑到拆点，确实有点难度，你需要把一天拆成干净餐巾和肮脏餐巾，因为脏餐巾是可以转化为干净餐巾的，你需要做到这个过程，所以你需要拆点，其次难想到的应该是如何解决快洗部和慢洗部每次只能洗一件并且需要洗多天，如何安排哪天的餐巾去洗就是个问题，所以可以考虑每天的餐巾都可以留到下一天，多连一条这样的边就可以解决这个问题了，别的都很好想了，至此解决。

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
#define min(a,b) (a<b?a:b)
#define rg register
int ss,tt,p,m,f,k,s,n,inf=1e9,cnt=1,sla[3000001],dis[30001],pre[300001],nxt[300001],v[300001],w[300001],h[300001],r[2001];
bool vis[10001];long long ans;
inline void add(int x,int y,int z,int u)
{
    pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt,v[cnt]=z,w[cnt]=u,
    pre[++cnt]=x,nxt[cnt]=h[y],h[y]=cnt,v[cnt]=0,w[cnt]=-u;
}
inline int dfs(int x,int flow)
{
    if(x==tt||!flow){ans+=1ll*dis[tt]*flow;return flow;}
    int f=flow;vis[x]=1;
    for(rg int i=h[x];i;i=nxt[i])
        if(v[i]&&!vis[pre[i]])
        {
            if(!(dis[x]+w[i]-dis[pre[i]]))
            {
                int y=dfs(pre[i],min(v[i],f));
                f-=y,v[i]-=y,v[i^1]+=y;
                if(!f)return flow;
            }
            else sla[pre[i]]=min(sla[pre[i]],dis[x]+w[i]-dis[pre[i]]);
        }
    return flow-f;
}
inline int aug()
{
    int mn=inf;
    for(rg int i=ss;i<=tt;i++)if(!vis[i])mn=min(mn,sla[i]),sla[i]=inf;
    if(mn==inf)return 1;
    for(rg int i=ss;i<=tt;i++)if(!vis[i])dis[i]+=mn;
    return 0;
}
int main()
{
    read(n),ss=0,tt=2*n+1;
    for(rg int i=1;i<=n;i++)read(r[i]);
    read(p),read(m),read(f),read(k),read(s);
    for(rg int i=1;i<=n;i++)
    {
        add(ss,i+n,r[i],0),add(ss,i,inf,p),add(i,tt,r[i],0);
        if(i<n)add(i+n,i+n+1,inf,0);
        if(i+k<=n)add(i+n,i+k,inf,s);
        if(i+m<=n)add(i+n,i+m,inf,f);
    }
    memset(sla,127,sizeof sla);
    while(1)
    {
        while(1){memset(vis,0,sizeof vis);if(!dfs(ss,inf))break;}
        if(aug())break;
    }
    printf("%lld
",ans);
}