P4012 深海机器人问题（线性规划网络优化，最大费用最大流）

题意：

深海资源考察探险队的潜艇将到达深海的海底进行科学考察。

潜艇内有多个深海机器人。潜艇到达深海海底后，深海机器人将离开潜艇向预定目标移动。

深海机器人在移动中还必须沿途采集海底生物标本。沿途生物标本由最先遇到它的深海机器人完成采集。

每条预定路径上的生物标本的价值是已知的，而且生物标本只能被采集一次。

本题限定深海机器人只能从其出发位置沿着向北或向东的方向移动，而且多个深海机器人可以在同一时间占据同一位置。

用一个P×Q 网格表示深海机器人的可移动位置。西南角的坐标为 (0,0)，东北角的坐标为 (Q,P)。

给定每个深海机器人的出发位置和目标位置，以及每条网格边上生物标本的价值。

计算深海机器人的最优移动方案， 使深海机器人到达目的地后，采集到的生物标本的总价值最高。

思路：

感觉这个描述不错，借来用下。

#include<bits/stdc++.h>
#define ll long long
#define INF 0x3f3f3f3f
using namespace std;
const int N = 1e5 + 10;
int n, m, s, t, h[N], cnt, dis[N], inq[N], pre[N], last[N], flow[N];
struct node {
    int v, w, c, nt;
} no[N];
void add(int u, int v, int w, int c) {
    no[cnt] = node{v, w, c, h[u]};
    h[u] = cnt++;
    no[cnt] = node{u, 0, -c, h[v]};
    h[v] = cnt++;
}
int spfa() {
    memset(inq, 0, sizeof inq);
    memset(dis, INF, sizeof dis);
    dis[s] = 0, inq[s] = 1, flow[s] = INF;
    queue<int> q;
    q.push(s);
    while(!q.empty()) {
        int u = q.front();
        q.pop(), inq[u] = 0;
        for(int i = h[u]; ~i; i = no[i].nt) {
            int v = no[i].v;
            if(no[i].w > 0 && dis[v] > dis[u] + no[i].c) {
                dis[v] = dis[u] + no[i].c;
                pre[v] = u, last[v] = i;
                flow[v] = min(flow[u], no[i].w);
                if(!inq[v])
                    inq[v] = 1, q.push(v);
            }
        }
    }
    return dis[t] != INF;
}
ll mcmf() {
    ll ans = 0;
    while(spfa()) {
        for(int i = t; i != s; i = pre[i]) {
            no[last[i]].w -= flow[t];
            no[last[i] ^ 1].w += flow[t];
        }
        ans += 1ll * flow[t] * dis[t];
    }
    return ans;
}
int a, b, P, Q;
int co(int x, int y) {
    return (x - 1) * Q + y;
}
int main() {
    memset(h, -1, sizeof h);
    scanf("%d%d%d%d", &a, &b, &P, &Q);
    P++, Q++;
    t = 1000;
    for(int x, i = 1; i <= P; i++)
        for(int j = 1; j < Q; j++) {
            scanf("%d", &x);
            add(co(i, j), co(i, j) + 1, 1, -x);
            add(co(i, j), co(i, j) + 1, INF, 0);
        }
    for(int x, j = 1; j <= Q; j++)
        for(int i = 1; i < P; i++) {
            scanf("%d", &x);
            add(co(i, j), co(i, j) + Q, 1, -x);
            add(co(i, j), co(i, j) + Q, INF, 0);
        }
    for(int k, x, y, i = 1; i <= a; i++) {
        scanf("%d%d%d", &k, &x, &y);
        x++, y++;
        add(s, co(x, y), k, 0);
    }
    for(int k, x, y, i = 1; i <= b; i++) {
        scanf("%d%d%d", &k, &x, &y);
        x++, y++;
        add(co(x, y), t, k, 0);
    }
    printf("%d\n", -mcmf());
    return 0;
}