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