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,j为1表示员工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,j。Wi,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]工作安排 费用流的主要内容,如果未能解决你的问题,请参考以下文章