【线性规划与网络流24题 10】餐巾计划
Time Limit:10000MS Memory Limit:65536K
Total Submit:36 Accepted:11
Case Time Limit:1000MS
Description
一个餐厅在相继的N 天里,每天需用的餐巾数不尽相同。假设第i天需要ri块餐巾(i=1,2,…,N)。餐厅可以购买新的餐巾,每块餐巾的费用为p分;或者把旧餐巾送到快洗部,洗一块需m天,其费用为f 分;或者送到慢洗部,洗一块需n 天(n>m),其费用为s<f 分。
每天结束时,餐厅必须决定将多少块脏的餐巾送到快洗部,多少块餐巾送到慢洗部,以及多少块保存起来延期送洗。但是每天洗好的餐巾和购买的新餐巾数之和,要满足当天的需求量。
试设计一个算法为餐厅合理地安排好N 天中餐巾使用计划,使总的花费最小。
编程任务:编程找出一个最佳餐巾使用计划.
Input
第1 行有6 个正整数N,p,m,f,n,s。N 是要安排餐巾使用计划的天数;p 是每块新餐巾的费用;m 是快洗部洗一块餐巾需用天数;f 是快洗部洗一块餐巾需要的费用;n是慢洗部洗一块餐巾需用天数;s是慢洗部洗一块餐巾需要的费用。
接下来的N 行是餐厅在相继的N 天里,每天需用的餐巾数。
Output
程序运行结束时,将餐厅在相继的N 天里使用餐巾的最小总花费输出
Sample Input
3 10 2 3 3 2
5
6
7
Sample Output
145
Hint
题目中所有数字都不超过1000
Source
感谢 Wo_ai_WangYuan 把原来的数据删掉又放上了和原来一模一样的数据
经典构图题。将每一天拆成两个点i, i’,一个记录一天的脏毛巾,一个记录一天的干净毛巾 ,加如下6条边:
(s, i, ri, p)——在第i天可以买至多ri个餐巾,每块p分;
(i, t, ri, 0)——第i天要用ri块餐巾;
(s, i’, ri, 0)——第i天用剩的ri块旧餐巾;
(i’, i+m, ∞, f)——第i天的旧餐巾送到快洗部,每块f分;
(i’, i+n, ∞, s)——第i天的旧餐巾送到慢洗部,每块s分;
(i’, i’+1, ∞, 0)——第i天的旧餐巾可以留到第i+1天再处理;
求一次最小费用流即为结果。
#include<cstdio>
#include<iostream>
#include<queue>
#include<cstring>
using namespace std;
const int maxn=2005,inf=1e9,maxm=2005;
int n,m,mincost,maxflow,op,ed,k,N,p,f,d;
int dis[maxn],path[maxn],last[maxm],t=3;
int need[maxn];
bool flag[maxn];
struct wr{
int a,b,c,w,NEXT;
wr(int a,int b,int w,int c,int NEXT):a(a),b(b),w(w),c(c),NEXT(NEXT){}
};
vector<wr>s;
void insert(int a,int b,int w,int c){
s.push_back(wr(a,b,w,c,last[a]));
last[a]=s.size()-1;
s.push_back(wr(b,a,-w,0,last[b]));
last[b]=s.size()-1;
}
struct wk{
int n;
void init(int n){
this->n =n;
memset(last,-1,sizeof(last));
}
bool find(int st,int nd){
int i;
queue<int>q;
memset(flag,0,sizeof(flag));
for(i=1;i<=n;i++)dis[i]=inf,path[i]=-1;
dis[st]=0,path[st]=0,flag[st]=1;
q.push(st);
while(!q.empty()){
int x=q.front();
q.pop();flag[x]=0;
for(i=last[x];i!=-1;i=s[i].NEXT){
wr& e=s[i];
if(e.c&&dis[x]+e.w<dis[e.b]){
dis[e.b]=dis[x]+e.w;
path[e.b]=i;
if(!flag[e.b]){
q.push(e.b) ;
flag[e.b]=1;
}
}
}
}
return dis[nd]<inf;
}
}spfa;
void addflow(){
int flow=inf,i;
for(i=ed;i!=op&&s[path[i]].a!=0;i=s[path[i]].a)
flow=min(flow,s[path[i]].c);
maxflow+=flow;
mincost+=dis[ed]*flow;
for(i=ed;i!=op;i=s[path[i]].a){
int x=path[i];
s[x].c-=flow;
s[x^1].c+=flow;
}
}
int main(){
cin>>N>>p>>m>>f>>n>>d;
int i,j,x;
op=2*N+1,ed=op+1;
spfa.init(ed);
for(i=1;i<=N;i++)cin>>need[i];
for(i=1;i<=N;i++)insert(op,i,0,need[i]);
for(i=1;i<N;i++)insert(i,i+1,0,inf);
for(i=N+1;i<=(2*N);i++){
insert(op,i,p,inf);
insert(i,ed,0,need[i-N]);
}
for(i=1;i<=N;i++){
if(i+m<=N)insert(i,i+m+N,f,inf);
if(i+n<=N)insert(i,i+n+N,d,inf);
}
while(spfa.find(op,ed))addflow();
cout<<mincost;
}