[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+1le 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
",MinCost);
	return 0;
}