BZOJ 1221: [HNOI2001] 软件开发

1221: [HNOI2001] 软件开发

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

#include<stdio.h>#include<algorithm>#include<string.h>#include<queue>using namespace std;const int maxm = 10005;const int maxn = 1000005;const int INF = 1e9 + 7;struct node{    int u, v, flow, cost, next;}edge[maxn];int dis[maxm], pre[maxm], head[maxm];int n, m, s, t, cnt, val;void init(){    cnt = 0, s = 0, t = n * 2 + 1;    memset(head, -1, sizeof(head));}void add(int u, int v, int w, int cost){    edge[cnt].u = u, edge[cnt].v = v;    edge[cnt].flow = w, edge[cnt].cost = cost;    edge[cnt].next = head[u], head[u] = cnt++;    edge[cnt].u = v, edge[cnt].v = u;    edge[cnt].flow = 0, edge[cnt].cost = -cost;    edge[cnt].next = head[v], head[v] = cnt++;}int bfs(){    queue<int>q;    for (int i = s;i <= t;i++) dis[i] = INF;    memset(pre, -1, sizeof(pre));    dis[s] = 0;q.push(s);    while (!q.empty())    {        int u = q.front();q.pop();        //printf("%d\n", u);        for (int i = head[u];i != -1;i = edge[i].next)        {            int v = edge[i].v;            //printf("%d\n", v);            if (dis[v] > dis[u] + edge[i].cost && edge[i].flow)            {                dis[v] = dis[u] + edge[i].cost;                pre[v] = i;                q.push(v);            }        }    }    if (dis[t] == INF) return 0;    return 1;}int MCMF(){    int minflow, ans = 0;    while (bfs())    {        minflow = INF;        for (int i = pre[t];i != -1;i = pre[edge[i].u])            minflow = min(minflow, edge[i].flow);        //printf("%d\n", minflow);        for (int i = pre[t];i != -1;i = pre[edge[i].u])        {            edge[i].flow -= minflow;            edge[i ^ 1].flow += minflow;        }        ans += dis[t] * minflow;    }    return ans;}int main(){    int i, j, k, sum, a, b, fa, fb;    scanf("%d%d%d%d%d%d", &n, &a, &b, &val, &fa, &fb);    init();    for (i = 1;i <= n;i++)    {        scanf("%d", &k);        add(s, i, k, 0);        add(i + n, t, k, 0);        add(s, i + n, INF, val);         if (i + 1 <= n) add(i, i + 1, INF, 0);        if (i + a + 1 <= n) add(i, i + a + n + 1, INF, fa);        if (i + b + 1 <= n) add(i, i + b + n + 1, INF, fb);    }    printf("%d\n", MCMF());    return 0;}