洛谷 P3980 [NOI2008]志愿者招募（网络流-费用流）

题目链接

题意：中文题不说了。

题解：这个题建边非常巧妙，1连源点，n+1连汇点。然后每点之间都连一条流量为inf-a[i]，费用为0的边。然后雇佣的人起点为u，终点为v，费用为val，就连一条u到v+1，流量inf（因为可以雇佣任意多的人），费用为val的边。

现在说说为啥这样建边是对的。现在我们从源点流到汇点一条流表示要雇佣一个人，所以至少要雇工作人数最多那天的人数 ， 但是现在有些天并不用雇那么多人，所以在这里有一些人是不用雇的，但是我们又必须跑那么多条流，怎么办？我们就在要不了那么多条人的地方建免费边，表示在这里有可以跑max-a[i]人条免费的边。 但是现在我们把每条边都建了inf-a[i]（为什么不建max-a[i]免费边？，其实可以的，只是要单独建反向边，代码里面解释）条免费边。 这样建边跑第一条流的流量肯定是min（inf-a[i]），且费用为0的边，然后剩下的费用为0的边的就是每天能跑免费边的数目了，这样就处理出每天能跑免费边的条数了，后面能跑免费边肯定直接跑免费边的流。

#pragma comment(linker, "/STACK:102400000,102400000")

#include <iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<map>
#include<queue>
#include<set>
#include<cmath>
#include<stack>
#include<string>

const int mod = 998244353;
const int maxn = 1e5 + 5;
const int inf = 1e9;
const long long onf = 1e18;
#define me(a, b) memset(a,b,sizeof(a))
#define lowbit(x) x&(-x)
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
#define PI 3.14159265358979323846
typedef long long ll;
typedef unsigned long long ull;
using namespace std;
int head[maxn], tot;
int minflow[maxn], pre[maxn];
int dis[maxn];
bool vis[maxn];
int n, m, st, ed;
int mincost;
int Max = 0, a[maxn];
struct node {
    int v, next, flow, val;
} maps[maxn << 2];

inline void add_edge(int u, int v, int flow, int val) {
    maps[tot].v = v;
    maps[tot].flow = flow;
    maps[tot].val = val;
    maps[tot].next = head[u];
    head[u] = tot++;
}

inline void ins(int u, int v, int flow, int val) {
    add_edge(u, v, flow, val);
    add_edge(v, u, 0, -val);
}

inline bool bfs() {
    fill(dis, dis + maxn, inf);
    dis[st] = 0, pre[st] = -1, minflow[st] = inf;
    queue<int> q;
    q.push(st);
    while (!q.empty()) {
        int u = q.front();
        q.pop(), vis[u] = 0;
        for (int i = head[u]; i != -1; i = maps[i].next) {
            int v = maps[i].v, flow = maps[i].flow, val = maps[i].val;
            if (flow && dis[v] > dis[u] + val) {
                dis[v] = dis[u] + val;
                pre[v] = i;
                minflow[v] = min(minflow[u], flow);
                if (!vis[v])
                    q.push(v), vis[v] = 1;
            }
        }
    }
    return dis[ed] < inf;
}

inline void MCMF() {
    while (bfs()) {
        for (int i = pre[ed]; i != -1; i = pre[maps[i ^ 1].v]) {
            maps[i].flow -= minflow[ed];
            maps[i ^ 1].flow += minflow[ed];
        }
        //     cout << minflow[ed] << " " << dis[ed] << endl;
        mincost += minflow[ed] * dis[ed];
    }
}

inline void init() {
    me(head, -1);
    tot = mincost = 0;
}

void build_inf() {///建inf-a[i]的免费边
    for (int i = 1; i <= n; i++) {
        ins(i, i + 1, inf - a[i], 0);
    }
    for (int i = 1; i <= m; i++) {
        int u, v, val;
        scanf("%d%d%d", &u, &v, &val);
        ins(u, v + 1, inf, val);
    }
    ins(st, 1, inf, 0), ins(n + 1, ed, inf, 0);
}

void build_max() {///建max-a[i]的免费边
    for (int i = 1; i <= n; i++) {
        add_edge(i, i + 1, Max - a[i], 0);
        add_edge(i + 1, i, inf - a[i], 0);///因为可能出现不跑这条免费边情况更优的情况，所以要给足反悔的余地
    }
    for (int i = 1; i <= m; i++) {
        int u, v, val;
        scanf("%d%d%d", &u, &v, &val);
        ins(u, v + 1, inf, val);
    }
    ins(st, 1, Max, 0), ins(n + 1, ed, Max, 0);
}

int main() {
    init();
    scanf("%d%d", &n, &m);
    st = 0, ed = n + 2;
    for (int i = 1; i <= n; i++) {
        scanf("%d", &a[i]);
        Max = max(Max, a[i]);
    }
    build_inf();
    ///build_max();
    MCMF();
    printf("%d\n", mincost);
    return 0;
}