bzoj2245 [SDOI2011]工作安排_2245: [sdoi2011]工作安排-CSDN博客

本文链接：https://blog.csdn.net/elijahqi/article/details/79091291

(http://www.elijahqi.win/2017/12/11/bzoj2245-sdoi2011%E5%B7%A5%E4%BD%9C%E5%AE%89%E6%8E%92/%20%E2%80%8E)
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 10

1 6

Sample Output

HINT
Source

第一轮day2

这是一份自己脑补的方法似乎有点慢在洛谷需要有o2加成才可过 bzoj15s过

http://blog.csdn.net/icefox_zhx/article/details/78768469 %%zhx大爷虐网络流题不眨眼啊 tjoi的省选题直接虐然而我是那个被题目虐的蒟蒻qwq

我的方法是针对每个人的每个愤怒值建一个点然后点与点之间的流量就是区间长度之差而费用呢就是下一段区间的愤怒值那么最后再把这些点去和商品连边最后商品连向汇点的时候费用为0 流量限制为商品数即可

#include<queue>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define inf 0x3f3f3f3f
#define N 2000
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if (T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0;char ch=gc();
    while(ch<'0'||ch>'9') ch=gc();
    while(ch<='9'&&ch>='0'){x=x*10+ch-'0';ch=gc();}
    return x;
}ll f[N];
int num=1,h[N],c[N],n,m,T,pre[N],path[N],s[300];bool flag[N],mp[300][300];
struct node{
    int y,next,z,c;
}data[N*N];
inline void insert1(int x,int y,int z,int c){
    data[++num].y=y;data[num].z=z;data[num].c=c;data[num].next=h[x];h[x]=num;
    data[++num].y=x;data[num].z=0;data[num].c=-c;data[num].next=h[y];h[y]=num;
}
inline bool spfa(){
    memset(f,127,sizeof(f));memset(flag,0,sizeof(flag));queue<int>q;q.push(0);flag[0]=1;f[0]=0;memset(pre,-1,sizeof(pre));
    while(!q.empty()){
        int x=q.front();flag[x]=0;q.pop();
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z,c=data[i].c;
            if(f[x]+c<f[y]&&z){
                f[y]=f[x]+c;pre[y]=x;path[y]=i;
                if(!flag[y]) flag[y]=1,q.push(y);
            }
        }
    }if(pre[T]==-1) return 0;else return 1;
}
int main(){
    freopen("bzoj2245.in","r",stdin);
    m=read();n=read();T=n+m+1;int cnt=T;
    for (int i=1;i<=n;++i) c[i]=read(),insert1(m+i,T,c[i],0);
    for (int i=1;i<=m;++i)
        for (int j=1;j<=n;++j) mp[i][j]=read();
    for (int i=1;i<=m;++i){
        int cc=read();insert1(0,i,inf,0);
        for (int j=1;j<=cc;++j) s[j]=read();s[cc+1]=inf;
        for (int j=1;j<=cc+1;++j) {
            int tmp=read();insert1(i,++cnt,s[j]-s[j-1],tmp);
            for (int z=1;z<=n;++z) if (mp[i][z]) insert1(cnt,m+z,inf,0);
        }
    }ll ans=0;
    while(spfa()){
        int minn=inf,now=T;
        while(now) minn=min(minn,data[path[now]].z),now=pre[now];now=T;
        while(now){ans+=(ll)data[path[now]].c*minn;data[path[now]].z-=minn;data[path[now]^1].z+=minn;now=pre[now];}
    }printf("%lld",ans);
    return 0;
}