Lưu ý: Các code mẫu dưới đây chỉ mang tính tham khảo và có thể không AC được bài tập này

Code mẫu của ladpro98

#include <bits/stdc++.h>

using namespace std;

const int N = 404;
const int LIMIT = 1e6 + 6;

short pos[N * N][N]; // optimize cache
short dp[N][N];
int a[N][N];
vector<pair<int, int> > V[LIMIT];
int n, m;

int main() {
    ios::sync_with_stdio(0);
    cin >> m >> n;
    for (int i = 1; i <= m; ++i) for (int j = 1; j <= n; ++j) cin >> a[i][j];
    for (int i = 1; i <= m; ++i) for (int j = 1; j <= n; ++j) V[a[i][j]].push_back(make_pair(i, j));
    int cnt = 0;
    for (int i = 1; i < LIMIT; ++i) if (!V[i].empty()) {
        for (int j = 0; j < V[i].size(); ++j) a[V[i][j].first][V[i][j].second] = cnt;
        V[i].clear();
        ++cnt;
    }
    int ans = 0;
    for (int k = 1; k <= m; ++k) {
        for (int i = n; i >= 1; --i) for (int j = i; j <= n; ++j) {
            dp[i][j] = max(dp[i][j], max(max(pos[a[k][i]][i], pos[a[k][j]][j]), max(dp[i + 1][j], dp[i][j - 1])));
            if (i < j) {
                dp[i][j] = max(dp[i][j], max(pos[a[k][i]][j], pos[a[k][j]][i]));
                if (a[k][i] == a[k][j]) dp[i][j] = k;
            }
            ans = max(ans, (j - i + 1) * (k - dp[i][j]));
        }
        for (int i = 1; i <= n; ++i) pos[a[k][i]][i] = k;
    }
    cout << ans << endl;
    return 0;
}