[BZOJ2245] [SDOI2011]工作安排

题目描述

Description

你的公司接到了一批订单。订单要求你的公司提供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,也为了使选手有更多的时间研究其他两道题目,你只需要输出最小的愤怒值之和就可以了。

Input

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

Output

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

Sample Input

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

Sample Output

24

HINT

题目分析

费用流
因为题目保证了 做的东西越多员工越生气 那么直接按照题意模拟即可

#include <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#include <algorithm>
using namespace std;
int n,m,s,t;
int head[10000],to[100000],net[100000],flo[100000],path[100000];
long long val[100000],dis[10000];
int tot=1;
void add(int x,int y,int v,long long c)
{
    net[++tot]=head[x],head[x]=tot,to[tot]=y,flo[tot]=v,val[tot]=c;
    net[++tot]=head[y];head[y]=tot,to[tot]=x,flo[tot]=0,val[tot]=-c;
}
bool vis[10000];
int spfa()
{
    memset(dis,0x3f,sizeof dis);
    memset(vis,0,sizeof vis);
    queue<int>q;
    vis[s]=1,q.push(s),dis[s]=0;
    while(q.size())
    {
        int nmp=q.front();
        q.pop(),vis[nmp]=0;
        for(int i=head[nmp];i;i=net[i])
            if(flo[i]>0&&dis[to[i]]>dis[nmp]+val[i])
            {
                dis[to[i]]=dis[nmp]+val[i],path[to[i]]=i^1;
                if(!vis[to[i]]) q.push(to[i]),vis[to[i]]=1;
            }
    }
    return dis[t]<0x3f3f3f3f;
}
long long mincost()
{
    long long ans=0;
    while(spfa())
    {
        int minflow=1<<30;
        for(int i=t;i!=s;i=to[path[i]]) minflow=min(minflow,flo[path[i]^1]);
        ans+=(long long)minflow*dis[t];
        for(int i=t;i!=s;i=to[path[i]])
            flo[path[i]]+=minflow,flo[path[i]^1]-=minflow;
    }
    return ans;
}
int T[10000];
int main()
{

    scanf("%d%d",&m,&n);
    s=0,t=m+n+1;
    for(int x,i=1;i<=n;i++)
        scanf("%d",&x),add(i+m,t,x,0);
    for(int i=1;i<=m;i++)
        for(int x,j=1;j<=n;j++)
        {
            scanf("%d",&x);
            if(x) add(i,j+m,1<<30,0);
        }
    for(int tmp,i=1;i<=m;i++)
    {
        scanf("%d",&tmp);
        for(int j=1;j<=tmp;j++) scanf("%d",&T[j]);
        for(int x,j=1;j<=tmp;j++) scanf("%d",&x),add(s,i,T[j]-T[j-1],x);
        scanf("%d",&tmp),add(s,i,1<<30,tmp);
    }
    printf("%lld",mincost());
    return 0;
}

发表评论

邮箱地址不会被公开。 必填项已用*标注