【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 1000000007using 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;}