Libre 6008 「网络流 24 题」餐巾计划 (网络流,最小费用最大流)

Libre 6008 「网络流 24 题」餐巾计划 (网络流,最小费用最大流)

Description

一个餐厅在相继的N天里,第i天需要Ri块餐巾(i=l,2,…,N)。餐厅可以从三种途径获得餐巾。
(1)购买新的餐巾,每块需p分;
(2)把用过的餐巾送到快洗部,洗一块需m天,费用需f分(f<p)。如m=l时,第一天送到快洗部的餐巾第二天就可以使用了,送慢洗的情况也如此。
(3)把餐巾送到慢洗部,洗一块需n天(n>m),费用需s分(s<f)。
在每天结束时,餐厅必须决定多少块用过的餐巾送到快洗部,多少块送慢洗部。在每天开始时,餐厅必须决定是否购买新餐巾及多少,使洗好的和新购的餐巾之和满足当天的需求量Ri,并使N天总的费用最小。

Input

共 n+1 行:
第1行为总天教和每块餐巾的新购费用p,快洗所需天数d1,快洗所需费用c1,慢洗所需天数d2,慢洗所需费用c2;
下n行为每天所需的餐巾块数;

Output

一行,最小的费用

Sample Input

3 10 2 3 3 2
5
6
7

Sample Output

145

Http

Libre:https://loj.ac/problem/6008

Source

网络流,最小费用最大流

解决思路

题目还算好理解,但是建图需要考虑一些东西。
餐巾不管在那一天洗价格都是一样的,所以没有必要提前把旧餐巾洗成新餐巾让后放在那里不用。
首先对于每一天i,我们将其拆成两个点i和i+n,我们把前面一个称为旧餐巾,后面一个称为新餐巾。另建立源点0和汇点n*2+1。
需要连接的边有:

1 源点到旧餐巾,容量为r[i](就是第i天的餐巾需求啦),代价为0,这条表表示这一天有r[i]条新的旧餐巾产生
2 源点到新餐巾,容量为无穷大,代价为新购餐巾的费用。这条边的意思就是新购餐巾
3 新餐巾到汇点,容量为R[i],代价为0。这条边使用来限制第i天只能花掉r[i]条餐巾,也代表着强制让剩余的餐巾向后流,即不让多余的餐巾提前洗掉
4 从第i天的旧餐巾到第i+1天的旧餐巾,容量为无穷大,代价为0。这条边代表的是延时送洗,即把餐巾送到后面几天去
5 从第i天的旧餐巾到第i+d1(即快洗所需时间)天的新餐巾,容量为无穷大,代价为快洗所需代价。这条边代表送到快洗部快洗
6 从第i天的旧餐巾到第i+d2(即慢洗所需时间)天的新餐巾,容量为无穷大,代价为慢洗所需代价。这条边代表送到慢洗部慢洗

然后跑一边spfa最小费用最大流即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

const int maxN=3100;
const int maxM=maxN*maxN*2;
const int inf=2147483647;

class Edge
{
public:
	int u,v,cost,flow;
};

int n,P,d1,c1,d2,c2;
int cnt=-1;
int Head[maxN];
int Next[maxM];
Edge E[maxM];
int Dist[maxN];
int path[maxN];
int Q[maxN*100];
bool inqueue[maxN];
int Flow[maxN];

void Add_Edge(int u,int v,int cost,int flow);
bool spfa();

int main()
{
	memset(Head,-1,sizeof(Head));
	scanf("%d%d%d%d%d%d",&n,&P,&d1,&c1,&d2,&c2);
	for (int i=1;i<=n;i++)
	{
		int r;
		scanf("%d",&r);
		Add_Edge(0,i,0,r);//源点到旧餐巾
		Add_Edge(0,i+n,P,inf);//从源点到新餐巾,即新购买餐巾
		Add_Edge(i+n,n*2+1,0,r);//从新餐巾到汇点,表示要用掉这么多餐巾
		if (i+1<=n)
			Add_Edge(i,i+1,0,inf);//延时送洗
		if (i+d1<=n)
			Add_Edge(i,i+d1+n,c1,inf);//送到快洗部
		if (i+d2<=n)
			Add_Edge(i,i+d2+n,c2,inf);//送到慢洗部
	}
	int Ans=0;
	while (spfa())
	{
		int now=n*2+1;
		int p=path[now];
		while (now!=0)
		{
			E[p].flow-=Flow[n*2+1];
			E[p^1].flow+=Flow[n*2+1];
			now=E[p].u;
			p=path[now];
		}
		Ans+=Flow[n*2+1]*Dist[n*2+1];
	}
	cout<<Ans<<endl;
	return 0;
}

void Add_Edge(int u,int v,int cost,int flow)
{
	cnt++;
	Next[cnt]=Head[u];
	Head[u]=cnt;
	E[cnt].u=u;
	E[cnt].v=v;
	E[cnt].cost=cost;
	E[cnt].flow=flow;

	cnt++;
	Next[cnt]=Head[v];
	Head[v]=cnt;
	E[cnt].u=v;
	E[cnt].v=u;
	E[cnt].cost=-cost;
	E[cnt].flow=0;
}

bool spfa()
{
	memset(Dist,-1,sizeof(Dist));
	memset(inqueue,0,sizeof(inqueue));
	int h=1,t=0;
	Q[1]=0;
	Dist[0]=0;
	inqueue[0]=1;
	Flow[0]=inf;
	do
	{
		t++;
		int u=Q[t];
		inqueue[u]=0;
		for (int i=Head[u];i!=-1;i=Next[i])
		{
			int v=E[i].v;
			if ((E[i].flow>0)&&((E[i].cost+Dist[u]<Dist[v])||(Dist[v]==-1)))
			{
				Dist[v]=Dist[u]+E[i].cost;
				path[v]=i;
				Flow[v]=min(Flow[u],E[i].flow);
				if (inqueue[v]==0)
				{
					h++;
					Q[h]=v;
					inqueue[v]=1;
				}
			}
		}
	}
	while (h!=t);
	if (Dist[n*2+1]==-1)
		return 0;
	return 1;
}
posted @ 2017-08-04 21:48  SYCstudio  阅读(...)  评论(...编辑  收藏