NKOI 1945 餐巾计划

【线性规划与网络流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;
}