【bzoj1221】[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

第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

题解
把每天分为二分图两个集合中的顶点Xi,Yi，建立附加源S汇T。
1、从S向每个Xi连一条容量为ni，费用为0的有向边。
2、从每个Yi向T连一条容量为ni，费用为0的有向边。
3、从S向每个Yi连一条容量为无穷大，费用为f的有向边。
4、从每个Xi向Xi+1(i+1<=N)连一条容量为无穷大，费用为0的有向边。
5、从每个Xi向Yi+a+1(i+a+1<=N)连一条容量为无穷大，费用为fa的有向边。
6、从每个Xi向Yi+b+1(i+b+1<=N)连一条容量为无穷大，费用为fb的有向边。
求网络最小费用最大流，费用流值就是要求的最小总花费。

代码

#include<bits/stdc++.h>
#define ll long long
#define inf 0x7fffffff
#define mod 1000000007
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,f,a,b,fa,fb,ans,tot,S,T;
int from[2005],q[2005],dis[2005],Head[2005];
bool vis[2005];
struct data{int from,to,next,w,c;}e[1000001];
inline void ins(int u,int v,int w,int c)
{
    tot++;
    e[tot].to=v;e[tot].from=u;
    e[tot].w=w;e[tot].c=c;
    e[tot].next=Head[u];Head[u]=tot;
}
inline void insert(int u,int v,int w,int c)
{
    ins(u,v,w,c);ins(v,u,0,-c);
}
inline bool spfa()
{
    for (int i=0;i<=T;i++) dis[i]=inf;
    queue<int>q;q.push(S);dis[S]=0;vis[S]=1;
    while (!q.empty())
    {
        int now=q.front();q.pop();
        for (int i=Head[now];i;i=e[i].next)
            if (dis[e[i].to]>dis[now]+e[i].c&&e[i].w)
            {
                from[e[i].to]=i;
                dis[e[i].to]=dis[now]+e[i].c;
                if (!vis[e[i].to])
                {
                    vis[e[i].to]=1;
                    q.push(e[i].to);
                }
            }
        vis[now]=0;
    }
    if (dis[T]==inf)return 0;return 1;
}
void mcf()
{
    int i=from[T],x=inf;
    while (i)
    {
        x=min(x,e[i].w);
        i=from[e[i].from];
    }
    i=from[T];
    while (i)
    {
        ans+=e[i].c*x;
        e[i].w-=x;e[i^1].w+=x;
        i=from[e[i].from];
    }
}
int main()
{
    n=read();a=read();b=read();
    f=read();fa=read();fb=read();
    S=2*n+1;T=S+1;tot=1;
    for (int i=1;i<=n;i++)
    {
        int x=read();
        insert(S,i,x,0);
        insert(i+n,T,x,0);
        insert(S,i+n,inf,f);
        if (i+1<=n) insert(i,i+1,inf,0);
        if (i+a+1<=n) insert(i,i+n+a+1,inf,fa);
        if (i+b+1<=n) insert(i,i+n+b+1,inf,fb);
    }
    while (spfa()) mcf();
    printf("%d",ans);
    return 0;
}