餐巾计划问题

问题描述

    一个餐厅在相继的N 天里， 每天需用的餐巾数不尽相同。 假设第i天需要ri块餐巾(i=1，2，…，N)。餐厅可以购买新的餐巾，每块餐巾的费用
为p分；或者把旧餐巾送到快洗部，洗一块需 m天，其费用为 f 分；或者送到慢洗部，洗一块需 n 天(n>m)，其费用为 s<f 分。每天结束时，餐厅必
须决定将多少块算法为餐厅合理地安排好N 天中餐巾使用计划，使总的花费最小。

编程任务

    编程找出一个最佳餐巾使用计划. 由文件input.txt提供输入数据。

数据输入

    文件第 1 行有6 个正整数 N,p,m,f,n,s。N 是要安排餐巾使用计划的天数；p 是每块新餐巾的费用；m 是快洗部洗一块餐巾需用天数；f 是快洗
部洗一块餐巾需要的费用；n是慢洗部洗一块餐巾需用天数；s是慢洗部洗一块餐巾需要的费用。

结果输出

    程序运行结束时，将餐厅在相继的 N 天里使用餐巾的最小总花费输出到文件 output.txt中。

输入文件示例

input.txt 
3 10 2 3 3 2
5 
6 
7 

输出文件示例

output.txt 
145 

题解

这道题是一道最小费用最大流的题目，然而有一些神犇写的是dinic，速度碾压众人QAQ
对于这道题，我们可以建立二分图（$x_i，y_i$）。
从源点到每个$x_i$、每个$y_i$到汇点均连一条容量为$r_i$、花费为0的边；
从源点到每个$y_i$，连一条容量为INF、花费为p的边；
对于每一个$x_i$(i < N),连一条到$x_i$$_+$$_1$、容量为INF、花费为0的边；
对于每一个$x_i$(i < N-m),连一条到$y_i$$_+$$_m$、容量为INF、花费为f的边；
对于每一个$x_i$(i < N-n),连一条到$y_i$$_+$$_n$、容量为INF、花费为s的边。
跑最小费用最大流，得到结果即为答案。

这道题很有代表性，建出的二分图虽然不像题目描述中那样，但是其中从源点直接流向每个$x_i$这种等量替代的方法很值得学习。

CODE：

#include<cstdio>
#include<cstring>
typedef long long ll;
const int INF=1e9;
struct queue
{
    int h,t;
    int a[20000];
    inline void clear(){h=1,t=0;}
    inline void push(int n){a[++t]=n;}
    inline int front(){return a[h];}
    inline void pop(){h++;}
    inline bool empty(){return h>t;}
}q;
struct edge
{
    int next,to,remain,cost;
}a[100000];
int head[20000];
int dis[20000];
int h[20000];
int f[20000];
int u[20000];
bool b[20000];
int n,x,num=1,S,T,c1,c2,c3,d1,d2;
inline int min(int a,int b){return a<b?a:b;}
inline void add(int x,int y,int z,int w)
{
    a[++num].next=head[x],a[num].to=y,a[num].remain=z,a[num].cost=w,head[x]=num;
    a[++num].next=head[y],a[num].to=x,a[num].cost=-w,head[y]=num;
}
inline int SPFA(ll &x)
{
    memset(dis,0x3f,sizeof(dis));
    dis[S]=0,b[S]=1,q.clear(),q.push(S),u[S]=INF;
    while(!q.empty())
    {
        int tmp=q.front();q.pop();
        b[tmp]=0;
        for(int i=head[tmp];i;i=a[i].next)
          if(a[i].remain&&dis[a[i].to]>dis[tmp]+a[i].cost)
          {
            dis[a[i].to]=dis[tmp]+a[i].cost;
            f[a[i].to]=i;
            u[a[i].to]=min(u[tmp],a[i].remain);
            if(!b[a[i].to]) q.push(a[i].to),b[a[i].to]=1;
          }
    }
    if(dis[T]>INF) return 0;
    x+=u[T]*dis[T];
    int p=T;
    while(p!=S)
    {
        a[f[p]].remain-=u[T];
        a[f[p]^1].remain+=u[T];
        p=a[f[p]^1].to;
    }
    return 1;
}
inline ll mcmf(){ll ans=0;while(SPFA(ans));return ans;}
int main()
{
    scanf("%d",&n);
    S=n<<1|1,T=S+1;
    for(int i=1;i<=n;i++)
      scanf("%d",&x),add(S,i,x,0),add(n+i,T,x,0);
    scanf("%d%d%d%d%d",&c1,&d1,&c2,&d2,&c3);
    for(int i=1;i<=n;i++)
      add(S,n+i,INF,c1);
    for(int i=1;i<n;i++)
      add(i,i+1,INF,0);
    for(int i=1;i<=n-d1;i++)
      add(i,n+i+d1,INF,c2);
    for(int i=1;i<=n-d2;i++)
      add(i,n+i+d2,INF,c3);
    printf("%lld",mcmf());
    return 0;
}