2245: [SDOI2011]工作安排

2245: [SDOI2011]工作安排(传送门点我QAQ)

---

**

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

---

**

Solution

**

把源点向所有货物连一条边 容量是需求量 cost为0

把货物向可以工作的人连一条边 容量为inf cost为0

把所有工人向汇点连好几条边……cost按照题目需求分段就好……

恩 就建图麻烦点QAQ

然而好久不打了……调了好久QAQ

---

**

Code

**

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <queue>
using namespace std;
const int maxn=255,inf=(1<<30);
int a[maxn][maxn],dis[maxn+maxn],from[maxn+maxn],inq[maxn+maxn],head[maxn+maxn],cnt=1,n,m,S,T;
int fw[10];
long long ans;
struct data{int to,from,next,v,c;}e[maxn*maxn];
void ins(int u,int v,int w,int c)
{
    cnt++;
    e[cnt].to=v;e[cnt].from=u;
    e[cnt].next=head[u];head[u]=cnt;
    e[cnt].v=w;e[cnt].c=c;
}
void insert(int u,int v,int w,int c)
{
    ins(u,v,w,c);ins(v,u,0,-c);
}
bool spfa()
{
    queue<int> q;
    for(int i=1;i<=T;i++)dis[i]=inf;
    inq[S]=1;dis[S]=0;q.push(S);
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=head[x];i;i=e[i].next)
            if(e[i].v>0&&dis[x]+e[i].c<dis[e[i].to])
            {
                dis[e[i].to]=dis[x]+e[i].c;
                from[e[i].to]=i;
                if(!inq[e[i].to])
                {
                    inq[e[i].to]=1;q.push(e[i].to);
                }
            }
        inq[x]=0;
    }
    return dis[T]!=inf;
}
void mcf()
{
    int x=inf;
    for(int i=from[T];i;i=from[e[i].from])
        x=min(x,e[i].v);
    for(int i=from[T];i;i=from[e[i].from])
    {
        e[i].v-=x;e[i^1].v+=x;
        ans+=e[i].c*x;
    }
}



int main()
{
    scanf("%d%d",&m,&n);S=n+m+1;T=S+1;
    for(int i=1,w;i<=n;i++)scanf("%d",&w),insert(S,i,w,0);
    for(int i=1,w;i<=m;i++)
        for(int j=1;j<=n;j++)
        {   
            scanf("%d",&w);
            if(w)insert(j,i+n,inf,0);
        }
    for(int i=1;i<=m;i++)
    {
        int s,w;
        scanf("%d",&s);
        for(int j=1;j<=s;j++)scanf("%d",&fw[j]);fw[++s]=inf;
        for(int j=1;j<=s;j++)scanf("%d",&w),insert(n+i,T,fw[j]-fw[j-1],w);
    }
    while(spfa())
        mcf();
    printf("%lld",ans);
    return 0;
}
//QQ:834026555 

---

---

倒要试南墙！