[BZOJ1061][Noi2008]志愿者招募 线性规划+费用流

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=1061

根据题意列方程，然后用网络流解线性规划。

题解直接贴ByVoid的吧，太神了：https://www.byvoid.com/blog/noi-2008-employee/

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
const int INF=1<<30;
int inline readint(){
    int Num;char ch;
    while((ch=getchar())<'0'||ch>'9');Num=ch-'0';
    while((ch=getchar())>='0'&&ch<='9') Num=Num*10+ch-'0';
    return Num;
}
int n,m;
int A[1010];
int S,T;
int to[30000],ne[30000],w[30000],c[30000],fir[1010],cnt=0;
void add(int a,int b,int x,int y=INF){
    to[cnt]=b,w[cnt]=x,c[cnt]=y,ne[cnt]=fir[a],fir[a]=cnt++;
    to[cnt]=a,w[cnt]=-x,c[cnt]=0,ne[cnt]=fir[b],fir[b]=cnt++;
}
int dis[1010],pre[1010];
bool in[1010];
queue <int> q;
bool Spfa(){
    memset(dis,127/3,sizeof(dis));
    q.push(S);
    in[S]=true;
    dis[S]=0;
    int u;
    while(!q.empty()){
        u=q.front();
        q.pop();
        in[u]=false;
        for(int i=fir[u];i!=-1;i=ne[i]){
            int v=to[i];
            if(dis[v]>dis[u]+w[i]&&c[i]){
                dis[v]=dis[u]+w[i];
                pre[v]=i;
                if(!in[v]){
                    q.push(v);
                    in[v]=true;
                }
            }
        }
    }
    return dis[T]!=dis[0];
}
void Mcmf(){
    int cost=0;
    pre[S]=-1;
    while(Spfa()){
        int f=INF;
        for(int i=pre[T];i!=-1;i=pre[to[i^1]]) f=min(f,c[i]);
        for(int i=pre[T];i!=-1;i=pre[to[i^1]]){
            c[i]-=f;
            c[i^1]+=f;
        }
        cost+=dis[T]*f;
    }
    printf("%d
",cost);
}
int main(){
    n=readint();
    m=readint();
    for(int i=1;i<=n;i++) A[i]=readint();
    memset(fir,-1,sizeof(fir));
    S=n+2;
    T=n+3;
    for(int i=1;i<=n+1;i++){
        if(A[i]>=A[i-1]) add(S,i,0,A[i]-A[i-1]);
        else add(i,T,0,A[i-1]-A[i]);
        if(i>1) add(i,i-1,0);
    }
    for(int i=1;i<=m;i++){
        int a=readint(),
            b=readint(),
            c=readint();
        add(a,b+1,c);
    }
    Mcmf();
    return 0;
}