Largest Common Submatrix

刚学的悬线法dp，热乎的。

悬线法不懂的话可以去了解一下，非常简单易懂。

这题的解法主要在判断怎样是公共上。
相信学习dp的朋友门应该做过一道题，是求公共子序列的就是用的映射方法，传送门

将第一个矩阵映射到一个有序矩阵上，再用映射数组去将第二个矩阵映射一遍，那么第二个矩阵中有序的子矩阵就是公共的。因为无论之前的矩阵是怎样的，这个映射数组将两个矩阵映射后总有一块是有序的且公共的。

//
// Created by acer on 2021/2/28.
//

#include <bits/stdc++.h>

using namespace std;
int lft[1100][1100];
int rgt[1100][1100], up[1100][1100];
int n, m;
int field[1100][1100];    //第一个矩阵
int cp[1100][1100];      //第二个矩阵
int mp[5000000];           //映射数组
int last[1100][1100];          //映射之后的第二个数组
signed main() {
    cin >> n >> m;
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            cin >> field[i][j];
        }
    }
    int cnt = 1;
    for (int i = 1; i <= n; ++i)
        for (int j = 1; j <= m; ++j)
            cp[i][j] = cnt++;
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            mp[field[i][j]] = cp[i][j];
        }
    }

    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            cin >> last[i][j];
            last[i][j] = mp[last[i][j]];
        }
    }

    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            up[i][j] = 1;
            rgt[i][j] = lft[i][j] = j;
        }
    }
    for (int i = 1; i <= n; ++i) {
        for (int j = 2; j <= m; ++j) {
            if (last[i][j-1]+1 == last[i][j])
            {

                lft[i][j] = lft[i][j-1];
            }
        }
    }
    for (int i = 1; i <= n; ++i) {
        for (int j = m; j > 1; --j) {
            if (last[i][j-1]+1 == last[i][j])
            {

                rgt[i][j-1] = rgt[i][j];
            }
        }
    }
    int ans = 0;
    for (int i = 1; i <= n; ++i) {
        for (int j = 1; j <= m; ++j) {
            if ( i>1 && last[i-1][j]+m == last[i][j])
            {
                up[i][j] = up[i-1][j]+1;
                lft[i][j] = max(lft[i-1][j],lft[i][j]);    //木桶效应
                rgt[i][j] = min(rgt[i-1][j],rgt[i][j]);
            }
            int a = rgt[i][j] - lft[i][j] + 1;
            ans = max(ans,a*up[i][j]);
        }
    }
    cout << ans << endl;
}