https://www.lydsy.com/JudgeOnline/problem.php?id=1061

https://www.luogu.org/problemnew/show/P3980

申奧成功後，布布經過不懈努力，終於成為奧組委下屬公司人力資源部門的主管。布布剛上任就遇到了一個難題：為即將啟動的奧運新項目招募一批短期誌願者。經過估算，這個項目需要N 天才能完成，其中第i 天至少需要Ai 個人。 布布通過了解得知，一共有M 類誌願者可以招募。其中第i 類可以從第Si 天工作到第Ti 天，招募費用是每人Ci 元。新官上任三把火，為了出色地完成自己的工作，布布希望用盡量少的費用招募足夠的誌願者，但這並不是他的特長！於是布布找到了你，希望你幫他設計一種最優的招募方案。

（先吐槽這題說其他數據均不超過2^31-1，那麽這題如果c=2^31-1，a=2^31-1，怕不是要爆longlong（滑稽））

一個非常奇妙的建圖費用流。

開始想法是上下界費用流，然後對於每一天往匯點流一遍，然後對於能不能把人收回去這一個問題表示十分糾結於是棄療。

考慮最開始S結點有INF，我們希望到T仍然有INF的流量，然而在經過每天的時候（比如在s天），很不幸，這條邊權為INF-a[i]，怎麽辦呢？我們就需要付錢c將一些流量傳送到t+1（顯然s~t天就會少需要這些人從而變相相當於這些人在s~t天內幹活了）。

（當然如果最後得到的流量並非INF就說明無解了。）

……是的，這就是這道題的建圖思路，十分的妙啊。

仔細思考一下，它有效規避了下界，人員持續工作與進入撤出的問題（實在不行感性畫圖一下）。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cctype>
#include<cstdio>
#include<queue>
#include<cmath>
using namespace std;
typedef long long ll;
const int N=1010;
const int M=3e4+5
 
;
const int INF=1e9;
inline int read(){
    int X=0,w=0;char ch=0;
    while(!isdigit(ch)){w|=ch==‘-‘;ch=getchar();}
    while(isdigit(ch))X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return w?-X:X;
}
struct node{
    int nxt,to,w,b;
}edge[M];
int head[N],cnt=-1;
inline void add(int u,int v,int w,int b){
    edge[++cnt].to=v;edge[cnt].w=w;edge[cnt].b=b;
    edge[cnt].nxt=head[u];head[u]=cnt;
    edge[++cnt].to=u;edge[cnt].w=0;edge[cnt].b=-b;
    edge[cnt].nxt=head[v];head[v]=cnt;
}
int dis[N];
bool vis[N];
inline bool spfa(int s,int t,int n){
    deque<int>q;
    memset(vis,0,sizeof(vis));
    for(int i=1;i<=n;i++)dis[i]=INF;
    dis[t]=0;q.push_back(t);vis[t]=1;
    while(!q.empty()){
        int u=q.front();
        q.pop_front();vis[u]=0;
        for(int i=head[u];i!=-1;i=edge[i].nxt){
            int v=edge[i].to;
            int b=edge[i].b;
            if(edge[i^1].w&&dis[v]>dis[u]-b){
                dis[v]=dis[u]-b;
                if(!vis[v]){
                    vis[v]=1;
                    if(!q.empty()&&dis[v]<dis[q.front()]){
                        q.push_front(v);
                    }else{
                        q.push_back(v);
                    }
                }
            }
        }
    }
    return dis[s]<INF;
}
int ans,cur[N];
int dfs(int u,int flow,int m){
    if(u==m){
        vis[m]=1;
        return flow;
    }
    int res=0,delta;
    vis[u]=1;
    for(int &e=cur[u];e!=-1;e=edge[e].nxt){
        int v=edge[e].to;
        int b=edge[e].b;
        if(!vis[v]&&edge[e].w&&dis[u]-b==dis[v]){
            delta=dfs(v,min(edge[e].w,flow-res),m); 
            if(delta){
                edge[e].w-=delta;
                edge[e^1].w+=delta;
                res+=delta;
                   ans+=delta*b;
                if(res==flow)break;
            }
        }
    }
    return res;
}
inline int costflow(int S,int T,int n){
    int flow=0;
    while(spfa(S,T,n)){
        do{
            for(int i=1;i<=n;i++)cur[i]=head[i];
            memset(vis,0,sizeof(vis));
            flow+=dfs(S,INF,T);
        }while(vis[T]);
    }
    return ans;
}
int main(){
    memset(head,-1,sizeof(head));
    int n=read(),m=read(),S=n+2,T=S+1;
    add(S,1,INF,0);
    for(int i=1;i<=n;i++){
        add(i,i+1,INF-read(),0);
    }
    add(n+1,T,INF,0);
    for(int i=1;i<=m;i++){
        int s=read(),t=read(),c=read();
        add(s,t+1,INF,c);
    }
    printf("%d\n",costflow(S,T,T));
    return 0;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+歡迎訪問我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++

BZOJ1061：[NOI2008]誌願者招募——題解