Vijos -- 1143 三取方格数

设有N*N的方格图，我们将其中的某些方格填入正整数，
而其他的方格中放入0。
某人从图得左上角出发，可以向下走，也可以向右走，直到到达右下角。
在走过的路上，他取走了方格中的数。（取走后方格中数字变为0）

此人从左上角到右下角共走3次，试找出3条路径，使得取得的数总和最大。

https://vijos.org/p/1143

好久没做题了QAQ。。。

一道多进程dp。从（1,1）走到（n,n）一共需要走2*n-2步，并且如果可以确定这一步所在的横坐标，那么就可以求出唯一确定的纵坐标。

用一个数组：dp[step][i][j][k] 表示在第step步，第一次走到了横坐标为i的点，第二次走到了横坐标为j的点，第三次走到了横坐标为k的点。所取的数的总和最大。

从第step-1步走到第step步共有8中走法。取它们的最大值加上这一步取的数的总和，然后再去重，得到的结果就是dp[step][i][j][k] 。

最终结果就是dp[2*n-2][n][n][n] + mat[1][1]（因为之前没有算（1,1）的值）

#include <map>
#include <queue>
#include <stack>
#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;

#define pb push_back
#define MP make_pair

typedef long long ll;
typedef unsigned long long ull;
typedef pair<int, int> PII;

const int inf = 0x3f3f3f3f;
const int mod = 1000000007;
const int maxn = 22;

int max8(int a, int b, int c, int d, int e, int f, int g, int h){
    int ret = 0;
    if(a > ret) ret = a;
    if(b > ret) ret = b;
    if(c > ret) ret = c;
    if(d > ret) ret = d;
    if(e > ret) ret = e;
    if(f > ret) ret = f;
    if(g > ret) ret = g;
    if(h > ret) ret = h;
    return ret;
}

int dp[maxn * 2][maxn][maxn][maxn], mat[maxn][maxn];
int n;
int main(){
    //freopen("out.txt", "w", stdout);
    scanf("%d", &n);
    memset(dp, 0, sizeof(dp));
    for(int i=1; i<=n; i++)
        for(int j=1; j<=n; j++)
        scanf("%d", &mat[i][j]);
    for(int step=1; step <= 2*n-2; step++){
        for(int i=1; i<=min(step, n)+1; i++){
            for(int j=1; j<=min(step, n)+1; j++){
                for(int k=1; k<=min(step, n)+1; k++){
                    dp[step][i][j][k] = max8(dp[step - 1][i][j][k], dp[step - 1][i][j][k - 1],
                                             dp[step - 1][i][j - 1][k], dp[step - 1][i][j - 1][k - 1],
                                             dp[step - 1][i - 1][j][k], dp[step - 1][i - 1][j][k - 1],
                                             dp[step - 1][i - 1][j - 1][k], dp[step - 1][i - 1][j - 1][k - 1]);
                    dp[step][i][j][k] += mat[i][min(step - i + 2, n)] + mat[j][min(step - j + 2, n)] + mat[k][min(step - k + 2, n)];
                    if(i == j) dp[step][i][j][k] -= mat[i][min(step - i + 2, n)];
                    if(i == k) dp[step][i][j][k] -= mat[i][min(step - i + 2, n)];
                    if(j == k) dp[step][i][j][k] -= mat[j][min(step - j + 2, n)];
                    if(i == j && i == k) dp[step][i][j][k] += mat[i][min(step - i + 2, n)];
//                    cout << "dp ";
                }
            }
        }
    }
    printf("%d\n", dp[2 * n - 2][n][n][n] + mat[1][1]);
    return 0;
}