网络流24题餐巾计划问题
问题描述
一个餐厅在相继的 N N 天里,每天需用的餐巾数不尽相同。假设第天需要 ri r i 块餐巾( i=1,2,…,N)。餐厅可以购买新的餐巾,每块餐巾的费用为 p p 分;或者把旧餐巾送到快洗部,洗一块需 m 天,其费用为 f 分;或者送到慢洗部,洗一块需天(n>m),其费用为 s(s<f) s ( s < f ) 。
每天结束时,餐厅必须决定将多少块脏的餐巾送到快洗部,多少块餐巾送到慢洗部,以及多少块保存起来延期送洗。但是每天洗好的餐巾和购买的新餐巾数之和,要满足当天的需求量。
试设计一个算法为餐厅合理地安排好 N N 天中餐巾使用计划,使总的花费最小。编程找出一个最佳餐巾使用计划。
输入
由标准输入提供输入数据。文件第1行有1个正整数,代表要安排餐巾使用计划的天数。
接下来的 N N 行是餐厅在相继的天里,每天需用的餐巾数。
最后一行包含5个正整数 p,m,f,n,s p , m , f , n , s 。 p p 是每块新餐巾的费用;是快洗部洗一块餐巾需用天数; f f 是快洗部洗一块餐巾需要的费用; 是慢洗部洗一块餐巾需用天数; s s <script type="math/tex" id="MathJax-Element-15">s</script>是慢洗部洗一块餐巾需要的费用。
输出
将餐厅在相继的N天里使用餐巾的最小总花费输出
Solution
一道非常有意思的费用流,建图方式非常巧(gui)妙(chu)
建图:
1.将每天拆成两部分上午和下午(可以感性地理解)
2.建立超级源S与每天上午连一条容量为x,费用为0的边
3.建立超级源T,将每天上午向T连一条容量为x,费用为0的边
4.从每天上午向第二天上午连一条无限走不要钱的边
5.从每天上午向第m天连一条容量为inf,费用为f的边
6.同5
7.从S像每天下午连一条容量为inf,费用为p的边
然后跑一边费用流
洛谷上可以测
#include<cstdio>
#include<cstring>
#include<queue>
#include<algorithm>
using namespace std;
#define int long long//毒瘤?呵。
const int inf=10000000000000000ll;
#define maxn 1000010
inline int read(){
int ret=0,ff=1;
char ch=getchar();
while(ch<'0'||ch>'9'){
if(ch=='-') ff=-ff;
ch=getchar();
}
while(ch>='0'&&ch<='9'){
ret=ret*10+ch-'0';
ch=getchar();
}
return ret;
}
struct Edge{
int u,v,w,c,next;
}E[maxn];
int ecnt=1;
int head[maxn],dis[maxn];
bool vis[maxn];
void addedge(int u,int v,int w,int c){
E[++ecnt].u=u;
E[ecnt].v=v;
E[ecnt].w=w;
E[ecnt].c=c;
E[ecnt].next=head[u];
head[u]=ecnt;
}
void Addedge(int u,int v,int w,int c){
addedge(u,v,w,c);
addedge(v,u,0,-c);
}
int n,m,s,t;
bool spfa(){
queue<int> q;
for(int i=1;i<=t;++i) dis[i]=inf;
dis[t]=0;
vis[t]=1;
q.push(t);
while(!q.empty()){
int x=q.front();
q.pop();
vis[x]=0;
for(int i=head[x];i;i=E[i].next){
int v=E[i].v;
if(!E[i^1].w||dis[x]+E[i^1].c>=dis[v])continue;
dis[v]=dis[x]+E[i^1].c;
if(!vis[v]){
vis[v]=1;
q.push(v);
}
}
}
return dis[s]!=inf;
}
int bns=0;
int dfs(int x,int nar){
vis[x]=1;
if(x==t) return nar;
int used=0;
for(int i=head[x];i;i=E[i].next){
int v=E[i].v;
if(!E[i].w||dis[v]!=dis[x]-E[i].c||vis[v]) continue;
int tmp=nar-used;
int flow=dfs(v,min(tmp,E[i].w));
used+=flow;
bns+=flow*E[i].c;
E[i].w-=flow;
E[i^1].w+=flow;
if(used==nar) return nar;
}
return used;
}
void zkw(){
while(spfa()){
vis[t]=1;
while(vis[t]){
memset(vis,0,sizeof(vis));
dfs(s,inf);
}
memset(vis,0,sizeof(vis));
}
printf("%lld\n",bns);
}
main(){
n=read();
s=n<<1|1,t=s+1;
for(int i=1;i<=n;++i){
int x=read();
Addedge(s,i,x,0);
Addedge(i+n,t,x,0);
}
int t1=read(),t2=read(),t3=read(),t4=read(),t5=read();
for(int i=1;i<=n;++i){
if(i<n) Addedge(i,i+1,inf,0);
if(i+t2<=n) Addedge(i,n+i+t2,inf,t3);
if(i+t4<=n) Addedge(i,n+i+t4,inf,t5);
Addedge(s,i+n,inf,t1);
}
zkw();
return 0;
}