初见安~这里是传送门:洛谷P4012 [网络流24题] 深海机器人问题
题解
题意很简单,多个点出发多个机器人,每条路径上有价值,但是只能收集一次而路径可以被多个机器人走,求最大可收集价值。
网络流建图很明显就如题意建网格图。但是对于每条边我们要整两条,一条容量为1,费用为价值,一条容量为INF,费用为0。换言之,一条拿来收集标本,一条拿来过路。
源点连向每个机器人出发的点,终点连向汇点。
还是比较简单,事后跑最大费用流就可以了。
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#define maxn 1005
#define maxm 100005
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
int read() {
int x = 0, f = 1, ch = getchar();
while(!isdigit(ch)) {if(ch == '-') f = -1; ch = getchar();}
while(isdigit(ch)) x = (x << 1) + (x << 3) + ch - '0', ch = getchar();
return x * f;
}
struct edge {int to, w, f, nxt;} e[maxm];
int head[maxn], k = 0;
void add(int u, int v, int w, int f) {
e[k] = {v, w, f, head[u]}; head[u] = k++;
e[k] = {u, 0, -f, head[v]}; head[v] = k++;
}
int n, m, a, b, S, T;
int id(int x, int y) {return x * m + y;};
int dis[maxn], pre[maxn], incf[maxn], ans = 0;
bool vis[maxn];
bool spfa() {
memset(dis, 0x3f, sizeof dis), dis[S] = 0;
memset(vis, 0, sizeof vis); vis[S] = true;
queue<int> q; incf[S] = INF, q.push(S);
while(q.size()) {
register int u = q.front(), v; q.pop(); vis[u] = false;
for(int i = head[u], v; ~i; i = e[i].nxt) {
v = e[i].to;
if(e[i].w && dis[u] + e[i].f < dis[v]) {
dis[v] = dis[u] + e[i].f;
incf[v] = min(incf[u], e[i].w); pre[v] = i;
if(!vis[v]) q.push(v), vis[v] = true;
}
}
}
if(dis[T] == INF) return false; return true;
}
void update() {
register int u = T, i;
while(u != S) {
i = pre[u];
e[i].w -= incf[T], e[i ^ 1].w += incf[T];
u = e[i ^ 1].to;
}
ans += 1ll * incf[T] * dis[T];
}
signed main() {
memset(head, -1, sizeof head);
a = read(), b = read(), n = read() + 1, m = read() + 1;
S = n * m + 5, T = n * m + 10;
register int w, x, y;//最大费用流=用负边权跑最小费用流。因为有费用为0的边。
for(int i = 0; i < n; i++) for(int j = 1; j < m; j++) w = read(), add(id(i, j - 1), id(i, j), 1, -w), add(id(i, j - 1), id(i, j), INF, 0);
for(int j = 0; j < m; j++) for(int i = 1; i < n; i++) w = read(), add(id(i - 1, j), id(i, j), 1, -w), add(id(i - 1, j), id(i, j), INF, 0);
for(int i = 1; i <= a; i++) w = read(), x = read(), y = read(), add(S, id(x, y), w, 0);
for(int i = 1; i <= b; i++) w = read(), x = read(), y = read(), add(id(x, y), T, w, 0);
while(spfa()) update();//最大费用流
printf("%lld\n", -ans);
return 0;
}
连续水了几篇最近做的题……
迎评:)
——End——