给定一个由 nnn 行数字组成的数字梯形如下图所示。
梯形的第一行有 mmm 个数字。从梯形的顶部的 mmm 个数字开始,在每个数字处可以沿左下或右下方向移动,形成一条从梯形的顶至底的路径。
分别遵守以下规则:
输入格式:
第 111 行中有 222 个正整数 mmm 和 nnn ,分别表示数字梯形的第一行有 mmm 个数字,共有 nnn 行。接下来的 nnn 行是数字梯形中各行的数字。
第 111 行有 mmm 个数字,第 222 行有 m+1m+1m+1 个数字,以此类推。
输出格式:
将按照规则 111 ,规则 222 ,和规则 333 计算出的最大数字总和并输出,每行一个最大总和。
输入样例#1: 复制
2 5
2 3
3 4 5
9 10 9 1
1 1 10 1 1
1 1 10 12 1 1
输出样例#1: 复制
66
75
77
1≤m,n≤201\leq m,n \leq 201≤m,n≤20
我只会第一问qwq。。
因为有数量的限制,考虑拆点建图,把每个点拆为$a_1$和$b_1$,两点之间连流量为$1$,费用为权值的边
从$b_i$向下方和右下的$a_1$连一条流量为$1$,费用为$0$边
从$S$向第一层的$a_1$连流量为$1$,费用为$0$的边,从$b_N$到$T$连流量为$1$,费用为$0$的边
对于第二问,因为没有点的个数的限制,那么就不用拆点了,直接向能到达的点连流量为$1$,费用为点权的边
对于第三问,直接把第二问中的所有边为流量设为$INF$(除了从$S$出发的)
// luogu-judger-enable-o2
/*
因为有数量的限制,考虑拆点建图,把每个点拆为$a_1$和$b_1$,两点之间连流量为$1$,费用为权值的边
从$b_i$向下方和右下的$a_1$连一条流量为$1$,费用为$0$边
从$S$向第一层的$a_1$连流量为$1$,费用为$0$的边,从$b_N$到$T$连流量为$1$,费用为$0$的边
对于第二问,因为没有点的个数的限制,那么就不用拆点了,直接向能到达的点连流量为$1$,费用为点权的边
对于第三问,直接把第二问中的所有边为流量设为$INF$(除了从$S$出发的)
*/
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
using namespace std;
const int MAXN = 1e5 + 10, INF = 1e9 + 10;
inline int read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = 1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N, M, S = 0, T = 1e5 + 1;
int a[21][21];
struct Edge {
int u, v, w, f, nxt;
}E[MAXN];
int head[MAXN << 1], num = 0;
inline void add_edge(int x, int y, int w, int f) {
E[num] = (Edge){x, y, w, f, head[x]};
head[x] = num++;
}
inline void AddEdge(int x, int y, int w, int f) {
add_edge(x, y, w, f);
add_edge(y, x, -w, 0);
}
int anscost, dis[MAXN], vis[MAXN], Pre[MAXN];
bool SPFA() {
memset(dis, -0x3f, sizeof(dis));
memset(vis, 0, sizeof(vis));
queue<int> q; q.push(S); dis[S] = 0;
while(!q.empty()) {
int p = q.front(); q.pop(); vis[p] = 0;
for(int i = head[p]; i !=- 1; i = E[i].nxt) {
int to = E[i].v;
if((dis[to] < dis[p] + E[i].w) && E[i].f > 0) {
dis[to] = dis[p] + E[i].w;
Pre[to] = i;
if(!vis[to]) q.push(to), vis[to] = 1;
}
}
}
return dis[T] > 0;
}
int F() {
int nowflow = INF;
for(int i = T; i != S; i = E[Pre[i]].u) nowflow = min(nowflow, E[Pre[i]].f);
for(int i = T; i != S; i = E[Pre[i]].u) E[Pre[i]].f -= nowflow, E[Pre[i] ^ 1].f += nowflow;
anscost += dis[T] * nowflow;
}
int MCMF() {
anscost = 0;
while(SPFA())
F();
return anscost;
}
int be[21][21], tot = 0, X;
void Solve1() {
memset(head, -1, sizeof(head)); num = 0;
for(int i = 1; i < N; i++) {
for(int j = 1; j <= M + i - 1; j++) {
AddEdge(be[i][j], be[i][j] + X, a[i][j], 1);
AddEdge(be[i][j] + X, be[i + 1][j], 0, 1);
AddEdge(be[i][j] + X, be[i + 1][j + 1], 0, 1);
}
}
for(int i = 1; i <= M; i++) AddEdge(S, be[1][i], 0, 1);
for(int i = 1; i <= N + M - 1; i++)
AddEdge(be[N][i], be[N][i] + X, a[N][i], 1),
AddEdge(be[N][i] + X, T, 0, 1);
printf("%d\n", MCMF());
}
void Solve2() {
memset(head, -1, sizeof(head)); num = 0;
for(int i = 1; i < N; i++) {
for(int j = 1; j <= M + i - 1; j++) {
AddEdge(be[i][j], be[i + 1][j + 1], a[i][j], 1);
AddEdge(be[i][j], be[i + 1][j], a[i][j], 1);
}
}
for(int i = 1; i <= M; i++) AddEdge(S, be[1][i], 0, 1);
for(int i = 1; i <= N + M - 1; i++) AddEdge(be[N][i], T, a[N][i], INF);
printf("%d\n", MCMF());
}
void Solve3() {
memset(head, -1, sizeof(head)); num = 0;
for(int i = 1; i < N; i++)
for(int j = 1; j <= M + i - 1; j++) {
AddEdge(be[i][j], be[i + 1][j + 1], a[i][j], INF);
AddEdge(be[i][j], be[i + 1][j], a[i][j], INF);
}
for(int i = 1; i <= M; i++) AddEdge(S, be[1][i], 0, 1);
for(int i = 1; i <= N + M - 1; i++) AddEdge(be[N][i], T, a[N][i], INF);
printf("%d\n", MCMF());
}
int main() {
memset(head, -1, sizeof(head));
M = read(); N = read(); X = (N + M - 1) * N;
for(int i = 1; i <= N; i++)
for(int j = 1; j <= M + i - 1; j++)
a[i][j] = read(), be[i][j] = ++tot;
Solve1();
Solve2();
Solve3();
return 0;
}
/*
*/