[BZOJ1221/Luogu2223][HNOI2001]软件开发

题目链接：

BZOJ1221

Luogu2223

不是很好理解的费用流。

把每个点拆成两个点，一个早上（没用毛巾），一个晚上（用完毛巾）。

连边：

• \(1.\)

早上与汇点连边，容量\(n_i\)，费用\(0\)，表示早上提供\(n_i\)条毛巾以供使用。

• \(2.\)

源点与晚上连边，容量\(n_i\)，费用\(0\)，表示早上用剩下的毛巾。

• \(3.\)

源点与早上连边，容量\(Inf\)，费用\(f\)，表示可以随便新买毛巾，一条\(f\)元。

• \(4.\)

若\(i\)不为最后一天，晚上与第二天晚上连边，容量\(Inf\)，费用\(0\)，表示可以不消毒，留着。

• \(5.\)

若\(i+a+1\le n\)（晚上用完第二天才开始消毒，所以\(+1\)），\(i\)天晚上向\(i+a+1\)天早上连边，容量\(Inf\)，费用\(fa\)，表示用\(A\)种消毒方式，可以在\(i+a+1\)天早上使用。

• \(6.\)

与\(5\)类似，消毒方式改为\(B\)。

话说我的费用流好慢。。

#include <queue>
#include <cstdio>
#include <cstring>

inline int Min(int a,int b){return a<b?a:b;}

int n,a,b,f,fa,fb,St,Ed,MaxFlow,MinCost;
int Head[2005],Next[40005],To[40005],Val[40005],Cos[40005],En=1;
int Dis[2005],Pre[2005],Ref[2005],Inq[2005];
const int Inf=0x3f3f3f3f;

inline void Add(int x,int y,int z,int w)
{
	Next[++En]=Head[x],To[Head[x]=En]=y,Val[En]=z,Cos[En]=+w;
	Next[++En]=Head[y],To[Head[y]=En]=x,Val[En]=0,Cos[En]=-w;
}

bool SPFA()
{
	std::queue<int> q;
	memset(Dis,0x3f,sizeof Dis);
	q.push(St),Dis[St]=0,Ref[St]=1<<30;
	for(int x,y;!q.empty();q.pop(),Inq[x]=0)
		for(int i=Head[x=q.front()];i;i=Next[i])
			if(Val[i]&&Dis[y=To[i]]>Dis[x]+Cos[i])
			{
				Dis[y]=Dis[x]+Cos[Pre[y]=i];
				Ref[y]=Min(Ref[x],Val[i]);
				if(!Inq[y])Inq[y]=1,q.push(y);
			}
	if(Dis[Ed]==Inf)return false;
	MaxFlow+=Ref[Ed],MinCost+=Ref[Ed]*Dis[Ed];
	for(int x=Ed,i;x!=St;x=To[i^1])
		Val[i=Pre[x]]-=Ref[Ed],Val[i^1]+=Ref[Ed];
	return true;
}

int main()
{
	scanf("%d%d%d%d%d%d",&n,&a,&b,&f,&fa,&fb);
	St=n<<1|1,Ed=St+1;
	for(int i=1,x;i<=n;++i)
	{
		scanf("%d",&x);
		Add(i,Ed,x,0);
		Add(St,i+n,x,0);
		Add(St,i,Inf,f);
		if(i<n)Add(i+n,i+1+n,Inf,0);
		if(i+a+1<=n)Add(i+n,i+a+1,Inf,fa);
		if(i+b+1<=n)Add(i+n,i+b+1,Inf,fb);
		//6种连边如上所述
	}
	while(SPFA());
	printf("%d\n",MinCost);
	return 0;
}