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;
}