【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

Source

第一轮day2

题解
源向每个工作连流量c[i]的边，费用0
每个人向汇点连s[i]+1条边
流量t[i][j]-t[i][j-1],费用w[i][j]
然后工作和人根据A矩阵互相连边
然后直接最小费用最大流

代码

#include<bits/stdc++.h>
#define N 500005
#define ll long long
#define inf 0x7fffffff
#define T 1001
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if (ch=='-') f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n,m;
int tot=1;ll ans;
int c[251],a[251][251],s[251],t[251][11],w[251][11];
int from[5005],dis[5005],head[5005];
bool inq[5005];
struct data{int from,to,next,v,c;}e[500001];
inline void ins(int u,int v,int w,int c)
{
    tot++;
    e[tot].to=v;e[tot].from=u;
    e[tot].v=w;e[tot].c=c;
    e[tot].next=head[u];head[u]=tot;
}
void insert(int u,int v,int w,int c)
{ins(u,v,w,c);ins(v,u,0,-c);}
void init()
{
    m=read();n=read();
    for (int i=1;i<=n;i++) c[i]=read();
    for (int i=1;i<=m;i++)
        for (int j=1;j<=n;j++)
            a[i][j]=read();
    for(int i=1;i<=m;i++)
    {
        s[i]=read();
        for(int j=1;j<=s[i];j++)
            t[i][j]=read();
        t[i][s[i]+1]=inf;
        for(int j=1;j<=s[i]+1;j++)
            w[i][j]=read();
    }
}
void build()
{
    for(int i=1;i<=n;i++)
        insert(0,i,c[i],0);
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n;j++)
            if(a[i][j])insert(j,n+i,inf,0);
    for(int i=1;i<=m;i++)
       for(int j=1;j<=s[i]+1;j++)
          insert(n+i,T,t[i][j]-t[i][j-1],w[i][j]);
}
bool spfa()
{
    for (int i=0;i<=n+m;i++) dis[i]=inf;dis[T]=inf;
    queue<int> q;
    q.push(0);dis[0]=0;inq[0]=1;
    while (!q.empty())
    {
        int now=q.front();q.pop();
        for (int i=head[now];i;i=e[i].next)
        {
            if (e[i].v&&dis[e[i].to]>dis[now]+e[i].c)
            {
                dis[e[i].to]=dis[now]+e[i].c;
                from[e[i].to]=i;
                if (!inq[e[i].to])
                {
                    inq[e[i].to]=1;
                    q.push(e[i].to);
                }
            }
        }
        inq[now]=0;
    }
    if (dis[T]==inf) return 0;return 1;
}
void mcf()
{
    int i=from[T],x=inf;
    while (i)
    {
        x=min(e[i].v,x);
        i=from[e[i].from];
    }
    i=from[T];
    while (i)
    {
        e[i].v-=x;e[i^1].v+=x;
        ans+=x*e[i].c;
        i=from[e[i].from];
    }
}
int main()
{
    init();
    build();
    while (spfa()) mcf();
    printf("%lld",ans);
    return 0;   
}