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

题目

某软件公司正在规划一项n天的软件开发计划，根据开发计划第i天需要ni个软件开发人员，为了提高软件开发人员的效率，公司给软件人员提供了很多的服务，其中一项服务就是要为每个开发人员每天提供一块消毒毛巾，这种消毒毛巾使用一天后必须再做消毒处理后才能使用。消毒方式有两种，A种方式的消毒需要a天时间，B种方式的消毒需要b天（b>a），A种消毒方式的费用为每块毛巾fA, B种消毒方式的费用为每块毛巾fB，而买一块新毛巾的费用为f（新毛巾是已消毒的，当天可以使用）；而且f>fA>fB。公司经理正在规划在这n天中，每天买多少块新毛巾、每天送多少块毛巾进行A种消毒和每天送多少块毛巾进行B种消毒。当然，公司经理希望费用最低。你的任务就是：为该软件公司计划每天买多少块毛巾、每天多少块毛巾进行A种消毒和多少毛巾进行B种消毒，使公司在这项n天的软件开发中，提供毛巾服务的总费用最低。

输入格式

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

输出格式

最少费用

输入样例

4 1 2 3 2 1

8 2 1 6

输出样例

38

题解

除了送洗时间多算一天
同 餐厅计划问题

#include<iostream>
#include<cstdio>
#include<cmath>
#include<queue>
#include<cstring>
#include<algorithm>
#define LL long long int
#define Redge(u) for (int k = h[u],to; k; k = ed[k].nxt)
#define REP(i,n) for (int i = 1; i <= (n); i++)
#define BUG(s,n) for (int i = 1; i <= (n); i++) cout<<s[i]<<' '; puts("");
using namespace std;
const int maxn = 3005,maxm = 1000005,INF = 1000000000;
inline int read(){
	int out = 0,flag = 1; char c = getchar();
	while (c < 48 || c > 57){if (c == '-') flag = -1; c = getchar();}
	while (c >= 48 && c <= 57){out = (out << 3) + (out << 1) + c - 48; c = getchar();}
	return out * flag;
}
int h[maxn],ne = 2;
struct EDGE{int from,to,nxt,f,w;}ed[maxm];
inline void build(int u,int v,int f,int w){
	ed[ne] = (EDGE){u,v,h[u],f,w}; h[u] = ne++;
	ed[ne] = (EDGE){v,u,h[v],0,-w}; h[v] = ne++;
}
int S,T,d[maxn],p[maxn],minf[maxn],vis[maxn];
queue<int> q;
int mincost(){
	int flow = 0,cost = 0;
	while (true){
		for (int i = S; i <= T; i++) d[i] = INF,vis[i] = false;
		q.push(S); d[S] = 0; vis[S] = true; minf[S] = INF;
		int u;
		while (!q.empty()){
			u = q.front(); q.pop();
			vis[u] = false;
			Redge(u) if (ed[k].f && d[to = ed[k].to] > d[u] + ed[k].w){
				d[to] = d[u] + ed[k].w; minf[to] = min(minf[u],ed[k].f); p[to] = k;
				if (!vis[to]) q.push(to),vis[to] = true;
			}
		}
		if (d[T] == INF) break;
		flow += minf[T]; cost += minf[T] * d[T];
		u = T;
		while (u != S){
			ed[p[u]].f -= minf[T];
			ed[p[u] ^ 1].f += minf[T];
			u = ed[p[u]].from;
		}
	}
	return cost;
}
int n,a,b,f,fa,fb;
int main(){
	n = read(); a = read(); b = read(); f = read(); fa = read(); fb = read();
	S = 0; T = 3 * n + 1;
	for (int i = 1; i <= n; i++){
		int x = read();
		build(i,T,x,0);
		build(S,i + n,INF,0);
		build(i + n,i,INF,f);
		build(i + n,i + 2 * n,x,0);
		if (i + a <= n) build(i + 2 * n,i + a + 1,x,fa);
		if (i + b <= n) build(i + 2 * n,i + b + 1,x,fb);
		if (i < n) build(i,i + 1,INF,0);
	}
	printf("%d\n",mincost());
	return 0;
}