bzoj2245 [SDOI2011]工作安排（费用流）

我可能还是有毒。本不需要把人拆成两个点的，直接源向人建边的时候建Si+1条边不就好了xxx。所以啦，源向人建Si+1条边，容量为$t_i-t_{i-1}$,费用为Wi。工作向汇建边，容量为Ci,费用为0，人和工作之间按矩阵建边，容量为inf，费用为0.跑最小费用最大流就是答案啦

#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define N 810
#define M 70000
#define inf 0x3f3f3f3f
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 m,n,T=800,h[N],num=1,a[N],dis[N],path[N];
ll ans=0;
bool inq[N];
struct edge{
    int to,next,w,c;
}data[M<<1];
inline void add(int x,int y,int w,int c){
    data[++num].to=y;data[num].next=h[x];h[x]=num;data[num].w=w;data[num].c=c;
    data[++num].to=x;data[num].next=h[y];h[y]=num;data[num].w=0;data[num].c=-c;
}
inline bool spfa(){
    queue<int>q;memset(dis,inf,sizeof(dis));memset(path,0,sizeof(path));
    q.push(0);inq[0]=1;dis[0]=0;
    while(!q.empty()){
        int x=q.front();q.pop();inq[x]=0;
        for(int i=h[x];i;i=data[i].next){
            int y=data[i].to;if(!data[i].w) continue;
            if(dis[x]+data[i].c<dis[y]){
                dis[y]=dis[x]+data[i].c;path[y]=i;
                if(!inq[y]) inq[y]=1,q.push(y);
            }
        }
    }return path[T];
}
int main(){
//  freopen("a.in","r",stdin);
    m=read();n=read();
    for(int i=1;i<=n;++i) add(m+m+i,T,read(),0);
    for(int i=1;i<=m;++i)
        for(int j=1;j<=n;++j){
            int op=read();if(!op) continue;
            add(m+i,m+m+j,inf,0);
        }
    for(int i=1;i<=m;++i){
        int si=read();for(int j=1;j<=si;++j) a[j]=read();
        for(int j=1;j<=si;++j) add(i,i+m,a[j]-a[j-1],read());
        add(i,i+m,inf,read());add(0,i,inf,0);
    }while(spfa()){
        int low=inf,now=T;
        while(path[now]) low=min(data[path[now]].w,low),now=data[path[now]^1].to;
        ans+=(ll)low*dis[T],now=T;
        while(path[now]) data[path[now]].w-=low,data[path[now]^1].w+=low,now=data[path[now]^1].to;
    }printf("%lld\n",ans);
    return 0;
}