AtCoder keyence2021 C エラーが出る理由を教えてください

Question

AtCoder keyence2021 C (https://atcoder.jp/contests/keyence2021) において以下のコードを書きました。方針は解説にあるものと同じです。入力例のうち、１と２では正解を出力するのですが、3では"abort() has been called." とエラーが出ます。これについて少し調べてみたのですが、いまいちよく分かりません。 このエラーについての解説と、できれば私のコードのどの部分がエラーを起こしているのかの指摘を頂きたいです。よろしくお願いします。

コード

#include <bits/stdc++.h>
#define _GLIBCXX_DEBUG
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define rep2(i, s, n) for (int i = (s); i < (int)(n); i++)
#define all(v) v.begin(), v.end()
#define pb push_back
#define sz(x) ((int)(x).size())
typedef long long ll;
using namespace std;
using Graph = vector<vector<int>>;
template<class T>bool chmax(T& a, const T& b) { if (a < b) { a = b; return 1; } return 0; }
template<class T>bool chmin(T& a, const T& b) { if (b < a) { a = b; return 1; } return 0; }
const ll INF = 1LL << 60;

int H, W, K; vector<vector<char>> grid; vector<vector<vector<ll>>> dp;

vector<ll> rec(int h, int w) {
    if (dp[h][w][0] != -1) return dp[h][w];

    vector<ll> res(H * W - K + 1, 0);

    if (h == 0 && w == 0) {
        res[0] = 1;
        return res;
    }

    if (h != 0) {
        if (grid[h - 1][w] == 'X' || grid[h - 1][w] == 'D') rep(i, H * W - K + 1) res[i] += rec(h - 1, w)[i] % 998244353;
        if (grid[h - 1][w] == 'O') rep(i, H * W - K) res[i + 1] += rec(h - 1, w)[i] % 998244353;
    }
    if (w != 0) {
        if (grid[h][w - 1] == 'X' || grid[h][w - 1] == 'R') rep(i, H * W - K + 1) res[i] += rec(h, w - 1)[i] % 998244353;
        if (grid[h][w - 1] == 'O') rep(i, H * W - K) res[i + 1] += rec(h, w - 1)[i] % 998244353;
    }
    
    return dp[h][w] = res;
}

int main() {
    cin.tie(0);
    ios::sync_with_stdio(false);

    cin >> H >> W >> K; grid = vector<vector<char>>(H, vector<char>(W, 'O'));
    rep(i, K) {
        int h, w; char c; cin >> h >> w >> c;

        grid[h - 1][w - 1] = c;
    }
    dp = vector<vector<vector<ll>>>(H, vector<vector<ll>>(W, vector<ll>(H * W - K + 1, -1)));

    vector<ll> fres = rec(H - 1, W - 1);

    ll ans = 0;
    rep(i, H* W - K + 1) {
        ans += fres[i] * pow(3, H * W - K - i) * pow(2, i);
        ans %= 998244353;
    }

    cout << ans << "\n";

    return 0;
}

入力例3
5000 5000 10
585 1323 R
2633 3788 X
1222 4989 D
1456 4841 X
2115 3191 R
2120 4450 X
4325 2864 X
222 3205 D
2134 2388 X
2262 3565 R

Answer 1

    dp = vector<vector<vector<ll>>>(H, vector<vector<ll>>(W, vector<ll>(H * W - K + 1, -1)));

llが8バイト、H=5000,W=5000,K=10で計算すると、この式で必要なメモリは5ペタバイト弱になります。
さすがにこのサイズを確保するのは無理だと思われます。