[HNOI2001][BZOJ1221] 软件开发

1221: [HNOI2001] 软件开发

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 868 Solved: 482
[Submit][Status][Discuss]

Description

某软件公司正在规划一项n天的软件开发计划，根据开发计划第i天需要ni个软件开发人员，为了提高软件开发人员的效率，公司给软件人员提供了很多的服务，其中一项服务就是要为每个开发人员每天提供一块消毒毛巾，这种消毒毛巾使用一天后必须再做消毒处理后才能使用。消毒方式有两种，A种方式的消毒需要a天时间，B种方式的消毒需要b天（b>a），A种消毒方式的费用为每块毛巾fA, B种消毒方式的费用为每块毛巾fB，而买一块新毛巾的费用为f（新毛巾是已消毒的，当天可以使用）；而且f>fA>fB。公司经理正在规划在这n天中，每天买多少块新毛巾、每天送多少块毛巾进行A种消毒和每天送多少块毛巾进行B种消毒。当然，公司经理希望费用最低。你的任务就是：为该软件公司计划每天买多少块毛巾、每天多少块毛巾进行A种消毒和多少毛巾进行B种消毒，使公司在这项n天的软件开发中，提供毛巾服务的总费用最低。

Input

第1行为n,a,b,f,fA,fB. 第2行为n1，n2，……，nn. （注：1≤f,fA,fB≤60，1≤n≤1000）

Output

最少费用

Sample Input

4 1 2 3 2 1
8 2 1 6

Sample Output

HINT

Source

最小费用最大流

首先无视费用，要求的就是每天的需求都要满足。而每天的毛巾都有这几个来源：昨天没洗，买了条新的，a天前消毒的，b天前消毒的。所以构造二分图：左边是今天用完的，右面是今天需要的。然后连边跑费用流就好了。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
#include<map>
#include<set>
#define INF 100000007
#define MAX 2000002
#define MAXN 2005
#define MOD 2002
using namespace std;
int next[MAX],list[MAX],from[MAX],flow[MAX],cost[MAX];
int head[MAXN],dis[MAXN],q[MAXN],path[MAXN];
int n,m,a,b,fa,fb,f,ans,cnt;
bool v[MAXN];
void insert(int x,int y,int w,int z)
{
    next[++cnt]=head[x];
    head[x]=cnt;
    from[cnt]=x;
    list[cnt]=y;
    flow[cnt]=w;
    cost[cnt]=z;
}
bool spfa()
{
    memset(v,0,sizeof(v));
    for (int i=0;i<=2*n+1;i++) dis[i]=INF;
    dis[0]=1; v[0]=1; q[1]=0;
    int t=0,w=1,x;
    while (t!=w)
    {
        t=(t+1)%MOD;
        x=q[t];
        for (int i=head[x];i;i=next[i])
            if (flow[i]&&dis[x]+cost[i]<dis[list[i]])
            {
                dis[list[i]]=dis[x]+cost[i];
                path[list[i]]=i;
                if (!v[list[i]])
                {
                    v[list[i]]=1;
                    w=(w+1)%MOD;
                    q[w]=list[i];
                }
            }
        v[x]=0;
    }
    return dis[2*n+1]!=INF;
}
void mcf()
{
    int x=INF,i;
    i=path[2*n+1];
    while (i)
    {
        x=min(x,flow[i]);
        i=path[from[i]];
    }
    i=path[2*n+1];
    while (i)
    {
        ans+=cost[i]*x;
        flow[i]-=x;
        flow[i^1]+=x;
        i=path[from[i]];
    }
}
int main()
{
    scanf("%d%d%d%d%d%d",&n,&a,&b,&f,&fa,&fb);
    cnt=1;
    for (int i=1;i<=n;i++) 
    {
        scanf("%d",&m);
        insert(0,i,m,0);
        insert(i,0,0,0);
        insert(i+n,2*n+1,m,0);
        insert(2*n+1,i+n,0,0);
        insert(0,n+i,INF,f);
        insert(n+i,0,0,-f);
    }
    for (int i=1;i<=n;i++) 
    {
        if (i+1<=n)
        {
            insert(i,i+1,INF,0);
            insert(i+1,i,0,0);
        }
        if (i+a+1<=n) 
        {
            insert(i,n+i+a+1,INF,fa);
            insert(n+i+a+1,i,0,-fa);
        }
        if (i+b+1<=n)
        {
            insert(i,n+i+b+1,INF,fb);
            insert(n+i+b+1,i,0,-fb);
        }
    }
    ans=0;
    while (spfa()) mcf();
    printf("%d",ans);
    return 0;
}