bzoj3265 志愿者招募加强版 线性规划与网络流（费用流）

题目分析

设$S_i$为第$i$天可以工作的志愿者的种类集合，$x_j$为第$j$种志愿者招募的个数，那么我们可以得到若干形如：

$$\sum_{j \in S_i} x_j \geq a_i$$
的式子。我们添加变量$y_i$使得原式变成：

$$(\sum_{j \in S_i} x_j) -y_i= a_i$$
设该式为 $P_i$。添加式子 $P_0:y_0=0$和 $P_{n+1}:y_{n+1}=0$后，式子之间两两相减可得若干形如 $P_{i-1}-P_i$的式子，记 $K_i=P_{i-1}-P_i$。并且将 $K_i$种等号右边的常数全部移到左边。
可以发现，假如第 $j$种志愿者工作的某一段时间是$[s_j,t_j]$，那么在 $K_{s_j}$中， $x_j$取负号。在第 $K_{t_j+1}$中， $x_j$取正号。所有这些式子中，变量 $x$和$y$一定拥有相等数量的取正号的式子和取负号的式子。
将一个 $K$式子看做点，一个取负号的变量或是常数看做流出，一个取正号的变量或常数看做流入，那么我们的目标是流量平衡，也就是等号右边为0。由于每种志愿者有不同的费用，所以使用费用流解决该问题。
对于一个$y_i$取负号的式子，要有一条容量为inf费用为0的边指向 $y_i$取正号的式子。也就是每个 $i$到$i+1$之间有一条容量为inf的边。
对于常数，一种方法是若正则从源点流过来，若负则流到汇点去。当然你也可以从点 $i$到点$i+1$之间流一条容量为 $-a_i$费用为0的边（也就是将那条容量为 $inf$的边的容量减少）
对于一个 $x_i$取负号的式子( $K_{s_j}$)，向一个 $x_i$取正号的式子( $K_{t_j+1}$)，连一条容量为inf费用为 $c_i$的边。

代码

#include<bits/stdc++.h>
using namespace std;
#define RI register int
int read() {
    int q=0;char ch=' ';
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9') q=q*10+ch-'0',ch=getchar();
    return q;
}
typedef long long LL;
const int N=1005,M=1000005,inf=0x3f3f3f3f;
int n,m,tot=1,S,T,ans;
int a[N],h[N],ne[M],to[M],flow[M],w[M],l[N],r[N];
int dis[N],pre[N],liu[N],inq[N];
void add(int x,int y,int z,int c) {
    to[++tot]=y,ne[tot]=h[x],h[x]=tot,flow[tot]=z,w[tot]=c;
    to[++tot]=x,ne[tot]=h[y],h[y]=tot,flow[tot]=0,w[tot]=-c;
}
int bfs() {
    for(RI i=S;i<=T;++i) dis[i]=0x3f3f3f3f,pre[i]=inq[i]=0;
    queue<int> q;
    dis[S]=0,liu[S]=inf,q.push(S);
    while(!q.empty()) {
        int x=q.front();q.pop(),inq[x]=0;
        for(RI i=h[x];i;i=ne[i])
            if(flow[i]>0&&dis[x]+w[i]<dis[to[i]]) {
                dis[to[i]]=dis[x]+w[i],pre[to[i]]=i;
                liu[to[i]]=min(flow[i],liu[x]);
                if(!inq[to[i]]) inq[to[i]]=1,q.push(to[i]);
            }
    }
    if(!pre[T]) return 0;
    ans+=liu[T]*dis[T];int x=T;
    while(x) {
        int kl=pre[x];
        flow[kl]-=liu[T],flow[kl^1]+=liu[T],x=to[kl^1];
    }
    return 1;
}
int main()
{
    int x,num;
    n=read(),m=read();
    S=0,T=n+1;
    add(S,1,inf,0);
    for(RI i=1;i<=n;++i) add(i,i+1,inf-read(),0);
    for(RI i=1;i<=m;++i) {
        num=read();
        for(RI j=1;j<=num;++j) l[j]=read(),r[j]=read();
        x=read();
        for(RI j=1;j<=num;++j) add(l[j],r[j]+1,inf,x);
    }
    while(bfs());
    printf("%d\n",ans);
    return 0;
}