问题描述
一个餐厅在相继的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
对于这道题,我们可以建立二分图(
xi,yi
)。
从源点到每个
xi
、每个
yi
到汇点均连一条容量为
ri
、花费为0的边;
从源点到每个
yi
,连一条容量为INF、花费为p的边;
对于每一个
xi
(i < N),连一条到
xi
+
1
、容量为INF、花费为0的边;
对于每一个
xi
(i < N-m),连一条到
yi
+
m
、容量为INF、花费为f的边;
对于每一个
xi
(i < N-n),连一条到
yi
+
n
、容量为INF、花费为s的边。
跑最小费用最大流,得到结果即为答案。
这道题很有代表性,建出的二分图虽然不像题目描述中那样,但是其中从源点直接流向每个 xi 这种等量替代的方法很值得学习。
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;
}