BZOJ1221[HNOI2001] 软件开发【拆点费用流】

$$[HNOI2001]软件开发$$

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 Format:

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

Output Format:

• 最少费用

Sample Input:

• 4 1 2 3 2 1
  8 2 1 6

Sample Output:

• 38

TJ：

很容易看出来是费用流，但是如何建边
考虑拆点，把每天拆成两个点，一个出点，一个入点
先说建边方式：

• S向i点的入点建边，流量为n[i]，费用为0
• S向i点的出点建边，流量为n[i]，费用为f
• i的出点向T建边，流量为n[i]，费用为0
• i的入点向i+a+1的出点建边，流量为INF，费用为fa
• i的入点向i+b+1的出点建边，流量为INF，费用为fb
• i的入点向i+1的入点建边，流量为INF，费用为0

首先因为每天的流量都要流到T点，但是这个流量又可以重复使用，所以源点流向i点有n[i]的流量
要保证每天都有足够的流量流出去，那就要让每天的出点向T点建n流量为[i]的边，同时源点向i的出点建费用为f的边，用来保证满流
然后就是第i天向第i+i天的流量是让没用过的毛巾继续使用。

#include<bits/stdc++.h>
using namespace std;
#define S 0
#define T MAXN-1
const int MAXN = 2222;
const int INF = 0x3f3f3f3f;
struct EDGE{
    int to,cap,fee,rev;
    EDGE(){}
    EDGE(int to, int cap, int fee, int rev){
        this->to = to;
        this->cap = cap;
        this->fee = fee;
        this->rev = rev;
    }
};
int n,da,db,ff,fa,fb;
int developer[MAXN],flow[MAXN],pre[MAXN],pree[MAXN],vis[MAXN],dist[MAXN];
vector<EDGE> G[MAXN];
void ADD_EDGE(int u, int v, int cap, int fee){
    G[u].push_back(EDGE(v,cap,fee,(int)G[v].size()));
    G[v].push_back(EDGE(u,0,-fee,(int)G[u].size()-1));
}
void build(){
    ADD_EDGE(S,1,INF,ff);
    for(int i = 1; i <= n; i++){
        if(i!=n) ADD_EDGE(i,i+1,INF,0);
        ADD_EDGE(i+n,T,developer[i],0);
        ADD_EDGE(S,i+n,developer[i],ff);
        ADD_EDGE(S,i,developer[i],0);
        if(i+da+1<=n) ADD_EDGE(i,i+da+1+n,INF,fa);
        if(i+db+1<=n) ADD_EDGE(i,i+db+1+n,INF,fb);
    }
}
bool SPFA(){
    memset(vis,0,sizeof(vis));
    memset(dist,INF,sizeof(dist));
    pre[T] = 0;
    dist[S] = 0;
    flow[S] = INF;
    queue<int> que;
    que.push(S);
    while(!que.empty()){
        int u = que.front();
        que.pop();
        vis[u] = 0;
        for(int i = 0; i < (int)G[u].size(); i++){
            EDGE e = G[u][i];
            if(!e.cap||dist[e.to]<=dist[u]+e.fee) continue;
            dist[e.to] = dist[u] + e.fee;
            pre[e.to] = u;
            pree[e.to] = i;
            flow[e.to] = min(flow[u],e.cap);
            if(!vis[e.to]){
                vis[e.to] = 1;
                que.push(e.to);
            }
        }
    }
    return pre[T]!=0;
}
int mincost_maxflow(){
    build();
    int cost = 0;
    while(SPFA()){
        cost += flow[T] * dist[T];
        int rt = T;
        while(rt!=S){
            int p = pre[rt];
            int d = pree[rt];
            G[p][d].cap -= flow[T];
            G[rt][G[p][d].rev].cap += flow[T];
            rt = p;
        }
    }
    return cost;
}
int main(){
    scanf("%d %d %d %d %d %d",&n,&da,&db,&ff,&fa,&fb);
    for(int i = 1; i <= n; i++) scanf("%d",&developer[i]);
    printf("%d\n",mincost_maxflow());
    return 0;
}