BZOJ_P1221 [HNOI2001] 软件开发(网络流+最小费用最大流)

BZOJ传送门
Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 1052 Solved: 576
[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

和网络流24题的餐巾问题基本一样，题解看前面的好了，辣个写的挺详细(其实是我懒得写)，注意几个细节
1.题目中没说可以遗留清洗，其实是可以的
2.这道题要清洗a天，所以在a+1天才能使用，b同理
3.前6个数的输入顺序不一样

#include<cstdio>
#include<climits>
#include<cstring>
#include<queue>
#include<vector>
#include<iostream>
using namespace std;
#define INF INT_MAX/3*2
#define N 1005*2
struct NetWork{
    struct Edge{int fr,to,cap,flow,cost;};
    vector<Edge> edge;vector<int> g[N];
    int fl[N],p[N],d[N];
    int n,a,b,f,fa,fb,x,m,s,t,cost,flow;

    void Add_Edge(int fr,int to,int cap,int cost){
        edge.push_back(Edge{fr,to,cap,0,cost});
        edge.push_back(Edge{to,fr,0,0,-cost});m=edge.size();
        g[fr].push_back(m-2);g[to].push_back(m-1);
    }

    int MaxFlowMinCost(){
        bool b[N];memset(fl,0,sizeof(fl));memset(b,0,sizeof(b));
        for(int i=1;i<N;i++) d[i]=INF;
        queue<int> q;q.push(s);fl[s]=INF,d[s]=0,b[s]=true;
        while(!q.empty()){
            int x=q.front();q.pop();b[x]=false;
            for(int i=0;i<g[x].size();i++){
                Edge &e=edge[g[x][i]];
                if(e.cap>e.flow&&d[e.to]>d[x]+e.cost){
                    d[e.to]=d[x]+e.cost,p[e.to]=g[x][i];
                    fl[e.to]=min(fl[x],e.cap-e.flow);
                    if(!b[e.to]){b[e.to]=true;q.push(e.to);}
                }
            }
        }
        if(d[t]==INF) return 0;
        flow+=fl[t];cost+=d[t]*fl[t];
        for(int x=t;x!=s;x=edge[p[x]].fr) edge[p[x]].flow+=fl[t],edge[p[x]^1].flow-=fl[t];
        return 1;
    }

    void solve(){
        scanf("%d%d%d%d%d%d",&n,&a,&b,&f,&fa,&fb);s=2001,t=2002;
        for(int i=1;i<=n;i++){
            scanf("%d",&x);
            Add_Edge(i,t,x,0);Add_Edge(s,i+1000,x,0);Add_Edge(s,i,INF,f);
            if(i<n) Add_Edge(i,i+1,INF,0);
            if(i+a<n) Add_Edge(i+1000,i+a+1,INF,fa);
            if(i+b<n) Add_Edge(i+1000,i+b+1,INF,fb);
        }
        while(MaxFlowMinCost());
        printf("%d",cost);
    }
}s;
int main(){
    s.solve();return 0;
}