bzoj2245[SDOI2011]工作安排 费用流

Posted GXZlegend

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj2245[SDOI2011]工作安排 费用流相关的知识,希望对你有一定的参考价值。

题目描述

你的公司接到了一批订单。订单要求你的公司提供n类产品,产品被编号为1~n,其中第i类产品共需要Ci件。公司共有m名员工,员工被编号为1~m员工能够制造的产品种类有所区别。一件产品必须完整地由一名员工制造,不可以由某名员工制造一部分配件后,再转交给另外一名员工继续进行制造。

我们用一个由0和1组成的m*n的矩阵A来描述每名员工能够制造哪些产品。矩阵的行和列分别被编号为1~m和1~n,Ai,j1表示员工i能够制造产品j,为0表示员工i不能制造产品j。

如果公司分配了过多工作给一名员工,这名员工会变得不高兴。我们用愤怒值来描述某名员工的心情状态。愤怒值越高,表示这名员工心情越不爽,愤怒值越低,表示这名员工心情越愉快。员工的愤怒值与他被安排制造的产品数量存在某函数关系,鉴于员工们的承受能力不同,不同员工之间的函数关系也是有所区别的。

对于员工i,他的愤怒值与产品数量之间的函数是一个Si+1段的分段函数。当他制造第1~Ti,1件产品时,每件产品会使他的愤怒值增加Wi,1,当他制造第Ti,1+1~Ti,2件产品时,每件产品会使他的愤怒值增加Wi,2……为描述方便,设Ti,0=0,Ti,si+1=+∞,那么当他制造第Ti,j-1+1~Ti,j件产品时,每件产品会使他的愤怒值增加Wi,j 1≤j≤Si+1。

你的任务是制定出一个产品的分配方案,使得订单条件被满足,并且所有员工的愤怒值之和最小。由于我们并不想使用Special Judge,也为了使选手有更多的时间研究其他两道题目,你只需要输出最小的愤怒值之和就可以了。

输入

第一行包含两个正整数m和n,分别表示员工数量和产品的种类数;

第二行包含n 个正整数,第i个正整数为Ci

以下m行每行n 个整数描述矩阵A;

下面m个部分,第i部分描述员工i的愤怒值与产品数量的函数关系。每一部分由三行组成:第一行为一个非负整数Si,第二行包含Si个正整数,其中第j个正整数为Ti,j,如果Si=0那么输入将不会留空行(即这一部分只由两行组成)。第三行包含Si+1个正整数,其中第j个正整数为Wi,jWi,j<Wi,j+1

输出

仅输出一个整数,表示最小的愤怒值之和。

样例输入

2 3
2 2 2
1 1 0
0 0 1
1
2
1 10
1
2
1 6

样例输出

24


题解

费用流

由于题目中限定了W(i,j)<W(i,j+1),因此可以直接拆边费用流。

那么建图很显然:S->每种产品,容量为Ci,费用为0;每种产品->能够生产它的人,容量为inf,费用为0;每个人->T连Si+1条边,第i条边容量为Ti - Ti-1,费用为Wi。

然后跑最小费用最大流即可。注意需要开long long

一个小优化:边上的费用只出现在与S/T之一相连的边中,这种情况下把带边权的边放到T一端能够大大减小时间复杂度,这也使得本题直接使用EK费用流即可AC。

#include <queue>
#include <cstdio>
#include <cstring>
#define N 510
#define M 1000010
using namespace std;
const int inf = 1 << 30;
queue<int> q;
int temp[10] , head[N] , to[M] , val[M] , cost[M] , next[M] , cnt = 1 , s , t , dis[N] , from[N] , pre[N];
void add(int x , int y , int v , int c)
{
	to[++cnt] = y , val[cnt] = v , cost[cnt] = c , next[cnt] = head[x] , head[x] = cnt;
	to[++cnt] = x , val[cnt] = 0 , cost[cnt] = -c , next[cnt] = head[y] , head[y] = cnt;
}
bool spfa()
{
	int x , i;
	memset(from , -1 , sizeof(from));
	memset(dis , 0x3f , sizeof(dis));
	dis[s] = 0 , q.push(s);
	while(!q.empty())
	{
		x = q.front() , q.pop();
		for(i = head[x] ; i ; i = next[i])
			if(val[i] && dis[to[i]] > dis[x] + cost[i])
				dis[to[i]] = dis[x] + cost[i] , from[to[i]] = x , pre[to[i]] = i , q.push(to[i]);
	}
	return ~from[t];
}
long long mincost()
{
	long long ans = 0;
	int i , k;
	while(spfa())
	{
		k = inf;
		for(i = t ; i != s ; i = from[i]) k = min(k , val[pre[i]]);
		ans += k * dis[t];
		for(i = t ; i != s ; i = from[i]) val[pre[i]] -= k , val[pre[i] ^ 1] += k;
	}
	return ans;
}
int main()
{
	int m , n , i , j , k , x;
	scanf("%d%d" , &m , &n) , s = 0 , t = n + m + 1;
	for(i = 1 ; i <= n ; i ++ ) scanf("%d" , &x) , add(s , i , x , 0);
	for(i = 1 ; i <= m ; i ++ )
	{
		for(j = 1 ; j <= n ; j ++ )
		{
			scanf("%d" , &x);
			if(x) add(j , i + n , inf , 0);
		}
	}
	for(i = 1 ; i <= m ; i ++ )
	{
		scanf("%d" , &k);
		for(j = 1 ; j <= k ; j ++ ) scanf("%d" , &temp[j]);
		temp[k + 1] = inf;
		for(j = 1 ; j <= k + 1 ; j ++ ) scanf("%d" , &x) , add(i + n , t , temp[j] - temp[j - 1] , x);
	}
	printf("%lld\n" , mincost());
	return 0;
}

 

 

以上是关于bzoj2245[SDOI2011]工作安排 费用流的主要内容,如果未能解决你的问题,请参考以下文章

BZOJ 2245[SDOI2011]工作安排

BZOJ 2245 SDOI 2011 工作安排 费用流

BZOJ2245 SDOI2011 工作安排 费用流

bzoj 2245 [SDOI2011]工作安排最小费用最大流

bzoj2245[SDOI2011]工作安排 费用流

BZOJ2245 [SDOI2011]工作安排 费用流