BZOJ 1221: [HNOI2001] 软件开发

1221: [HNOI2001] 软件开发

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 1459  Solved: 809
[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

38

HINT

Source

分析：

感谢机智的叶子童鞋...

这题一开始就是用费用流写的...后来叶子说可以用三分+贪心过掉...跑的飞起...(然而我并没有用堆优化...)

第一次写三分...

以下来自叶子...

显然费用与买的毛巾数成单峰函数(如果不够用函数值就是+INF)......并且无论什么时候买毛巾都和第一天买是等价的......
所以可以三分第一天买的毛巾数，贪心求函数值，具体的贪心策略就是能用买的新毛巾就用，不能用就从第一天往后找用完没洗的慢洗，慢洗完了不够用就从昨天往前快洗.....正确性很显然是吧......

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
//by NeighThorn
#define inf 0x3f3f3f3f
using namespace std;
//眉眼如初，岁月如故
 
const int maxn=1000+5;
 
int n,a,b,f,fa,fb,sum,s[maxn],re[maxn],num[maxn];
 
inline int calc(int x){
    int ans=x*f;
    memcpy(s,num,sizeof(s));
    memset(re,0,sizeof(re));
    for(int i=1;i<=n;i++){
        if(x){
            int lala=min(x,s[i]);
            x-=lala,s[i]-=lala,re[i]+=lala;
        }
        for(int j=1;j<=i-b-1&&s[i];j++){
            int lala=min(re[j],s[i]);
            re[j]-=lala,s[i]-=lala,ans+=fb*lala,re[i]+=lala;
        }
        for(int j=i-a-1;j>=1&&s[i];j--){
            int lala=min(re[j],s[i]);
            re[j]-=lala,s[i]-=lala,ans+=fa*lala,re[i]+=lala;
        }
        if(s[i])
            return inf;
    }
    return ans;
}
 
signed main(void){
    scanf("%d%d%d%d%d%d",&n,&a,&b,&f,&fa,&fb);sum=0;
    for(int i=1;i<=n;i++)
        scanf("%d",&num[i]),sum+=num[i];
    int l=1,r=sum;
    while(l<r-1){
        int mid=(l+r)>>1,Mid=(mid+r)>>1;
        if(calc(mid)<calc(Mid))
            r=Mid;
        else
            l=mid;
    }
    printf("%d
",min(calc(l),calc(r)));
    return 0;   
}//Cap ou pas cap. Pas cap.

---

By NeighThorn