餐巾计划问题

一个餐厅在相继的N天里,第i天需要Ri块餐巾(i=l,2,…,N)。餐厅可以从三种途径获得餐巾。
(1)购买新的餐巾,每块需p分;
(2)把用过的餐巾送到快洗部,洗一块需m天,费用需f分(f<p)。如m=l时,第一天送到快洗部的餐巾第二天就可以使用了,送慢洗的情况也如此。
(3)把餐巾送到慢洗部,洗一块需n天(n>m),费用需s分(s<f)。
在每天结束时,餐厅必须决定多少块用过的餐巾送到快洗部,多少块送慢洗部。在每天开始时,餐厅必须决定是否购买新餐巾及多少,使洗好的和新购的餐巾之和满足当天的需求量Ri,并使N天总的费用最小
输入格式:
输入文件共3行,第1行为总天数;第2行为每天所需的餐巾块数;第3行为每块餐巾的新购费用p,快洗所需天数m,快洗所需费用f,慢洗所需天数n,慢洗所需费用s。
输出格式:
输出文件共1行为最小的费用。

这道题比起数字梯形要难上不少。。在方格取数里我们学到了把一个点拆成2个点,两个点之间连两条边一个限制费用次数,另外一个正常走。
这道题中,我们汇点显然连得是每天需求的盘子。但是中间怎么弄呢。。
把一个点拆成两个点后,我们来考虑,源点和汇点分别连入还是出(例如最小路径覆盖中源点连出)

在这题中,源点如果连入点(一天的开始),那么一天的结束连汇点。可是结束之后我还要去洗盘子啊。。
怎么办?似乎没法办。
我们把一天的结束连到汇点,一天的开始连到源点。那么流入源点的东西是什么呢?显然还是那些需求的盘子(等等我写到这的时候突然感觉,需求的盘子明明就是应该由一天的结束来决定的啊,所以应该流入汇点)。
那么从源点流入结束的是什么,还是那些盘子,不过是已经用过需要洗的。
所以结束可以向开始连边,就是洗盘子~费用就是洗盘子的费用,容量是INF(因为最大流是根据流出即使用的盘子来决定的)。
。。那么问题来了,这样的话盘子不够(达不到满流)怎么办(需要买盘子):直接由S连入一天开始,费用是买盘子,可以买INF个。
好的交上去,WA了8个点。。。为什么。
我们来找一下逻辑的漏洞。一天结束洗好的盘子流入了另一天的新盘子中,看似没什么问题啊。。但是,如果那一天所需的盘子小于洗好的盘子的,剩下的盘子怎么办。。还要流入别的地方呀。。所以,我们把这一天的盘子留到明天洗,向下一天连一条容量为INF,费用为0的边。(为什么不会先走这个?因为先走这个显然不能达成最大流)。

所以这道题就结束辣。

#include<bits/stdc++.h>
using namespace std;

#define ll long long

const ll INF=1e9;
const int MAXN=2e5+5;
const int base=1e5;

struct edge{
    ll u,to,next,w,c;
}e[MAXN<<1];

int cnt=1,head[MAXN];
inline void add(ll u,ll v,ll w,ll cost){e[++cnt]=(edge){u,v,head[u],w,cost},head[u]=cnt;}

ll dis[MAXN],pre[MAXN];
bool vis[MAXN];
queue<int>q;

ll n,p,s,t;
ll dayfast,flowfast,dayslow,flowslow;

bool SPFA(int x){
    memset(pre,0,sizeof(pre));
    memset(vis,0,sizeof(vis));
    for(int i=s;i<=t;i++)dis[i]=INF;
    q.push(x);dis[x]=0;vis[x]=1;
    while(q.size()){
        int u=q.front();q.pop();
        vis[u]=0;
        for(int i=head[u];i;i=e[i].next){
            int v=e[i].to,w=e[i].c;
            if(e[i].w){
                if(dis[u]+w<dis[v]){
                    dis[v]=dis[u]+w;pre[v]=i;
                    if(!vis[v]){
                        vis[v]=1;q.push(v);
                    }
                }
            }
        }
    }
    if(dis[t]==INF)return 0;
    return 1;   
}

ll MF(){
    ll cost=0;
    while(SPFA(s)){
        ll tem=INF;
        for(int i=pre[t];i;i=pre[e[i].u]){
            tem=min(e[i].w,tem);
        }
        for(int i=pre[t];i;i=pre[e[i].u]){
            e[i].w-=tem;e[i^1].w+=tem;cost+=tem*e[i].c;
        }
    }
    return cost; 
}

int main(){
    scanf("%lld",&n);
    int tem;
    s=0,t=n*2+1+base;
    for(int i=1;i<=n;i++){
        scanf("%lld",&tem);
        add(s,i,tem,0),add(i,s,0,0);
        add(i+base,t,tem,0),add(t,i+base,0,0);
    }
    scanf("%lld%lld%lld%lld%lld",&p,&dayfast,&flowfast,&dayslow,&flowslow);
    for(int i=1;i<=n;i++){
        add(s,i+base,INF,p),add(i+base,s,0,-p);
        if(i+dayslow<=n)add(i,i+dayslow+base,INF,flowslow),add(i+dayslow+base,i,0,-flowslow);
        if(i+dayfast<=n)add(i,i+dayfast+base,INF,flowfast),add(i+dayfast+base,i,0,-flowfast);
        if(i+1<=n)add(i,i+1,INF,0),add(i+1,i,0,0);
    }
    printf("%lld\n",MF());
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值