Problem p02839 (C++)

Buggy submission (status: Runtime Error):
#include <bits/stdc++.h>
using namespace std;

// Define
using ll = long long;
using ull = unsigned long long;
using ld = long double;
template <class T> using pvector = vector<pair<T, T>>;
template <class T>
using rpriority_queue = priority_queue<T, vector<T>, greater<T>>;
constexpr const ll dx[4] = {1, 0, -1, 0};
constexpr const ll dy[4] = {0, 1, 0, -1};
constexpr const ll MOD = 1e9 + 7;
constexpr const ll mod = 998244353;
constexpr const ll INF = 1LL << 60;
constexpr const ll inf = 1 << 30;
constexpr const char rt = '\n';
constexpr const char sp = ' ';
#define mp make_pair
#define mt make_tuple
#define pb push_back
#define eb emplace_back
#define elif else if
#define all(a, v, ...)                                                         \
  ([&](decltype((v)) w) { return (a)(begin(w), end(w), ##__VA_ARGS__); })(v)
#define fi first
#define se second

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 (a > b) {
    a = b;
    return 1;
  }
  return 0;
}

// Debug
#define debug(...)                                                             \
  {                                                                            \
    cerr << __LINE__ << ": " << #__VA_ARGS__ << " = ";                         \
    for (auto &&X : {__VA_ARGS__})                                             \
      cerr << "[" << X << "] ";                                                \
    cerr << rt;                                                                \
  }

#define dump(a, h, w)                                                          \
  {                                                                            \
    cerr << __LINE__ << ": " << #a << " = [" << rt;                            \
    rep(i, h) {                                                                \
      rep(j, w) cerr << a[i][j] << sp;                                         \
      cerr << rt;                                                              \
    }                                                                          \
    cerr << "]" << rt;                                                         \
  }

#define vdump(a, n)                                                            \
  {                                                                            \
    cerr << __LINE__ << ": " << #a << " = [";                                  \
    rep(i, n) if (i) cerr << sp << a[i];                                       \
    else cerr << a[i];                                                         \
    cerr << "]" << rt;                                                         \
  }

// Loop
#define inc(i, a, n) for (ll i = (a), _##i = (n); i <= _##i; ++i)
#define dec(i, a, n) for (ll i = (a), _##i = (n); i >= _##i; --i)
#define rep(i, n) for (ll i = 0, _##i = (n); i < _##i; ++i)
#define each(i, a) for (auto &&i : a)

// Stream
#define fout(n) cout << fixed << setprecision(n)
struct io {
  io() { cin.tie(nullptr), ios::sync_with_stdio(false); }
} io;

// Speed
#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

// Math
inline constexpr ll gcd(const ll a, const ll b) {
  return b ? gcd(b, a % b) : a;
}
inline constexpr ll lcm(const ll a, const ll b) { return a / gcd(a, b) * b; }

inline constexpr ll modulo(const ll n, const ll m = MOD) {
  ll k = n % m;
  return k + m * (k < 0);
}
inline constexpr ll chmod(ll &n, const ll m = MOD) {
  n %= m;
  return n += m * (n < 0);
}
inline constexpr ll mpow(ll a, ll n, const ll m = MOD) {
  ll r = 1;
  rep(i, 64) {
    if (n & (1LL << i))
      r *= a;
    chmod(r, m);
    a *= a;
    chmod(a, m);
  }
  return r;
}
inline ll inv(const ll n, const ll m = MOD) {
  ll a = n, b = m, x = 1, y = 0;
  while (b) {
    ll t = a / b;
    a -= t * b;
    swap(a, b);
    x -= t * y;
    swap(x, y);
  }
  return modulo(x, m);
}

bool DP[80][80][20000 + 1];
signed main() {
  ll H, W;
  cin >> H >> W;
  ll A[H][W], B[H][W];
  rep(i, H) rep(j, W) cin >> A[i][j];
  rep(i, H) rep(j, W) cin >> B[i][j];
  rep(i, H + 1) rep(j, W + 1) rep(k, 20000 + 1) DP[i][j][k] = 0;
  DP[0][0][10000] = 1;
  rep(i, H) rep(j, W) rep(k, 20000 + 1) {
    if (i + 1 <= H - 1) {
      if (0 <= k + A[i][j] - B[i][j] && k + A[i][j] - B[i][j] <= 20000) {
        DP[i + 1][j][k + A[i][j] - B[i][j]] |= DP[i][j][k];
      }
      if (0 <= k + B[i][j] - A[i][j] && k + B[i][j] - A[i][j] <= 20000) {
        DP[i + 1][j][k + B[i][j] - A[i][j]] |= DP[i][j][k];
      }
    }
    if (j + 1 <= W - 1) {
      if (0 <= k + A[i][j] - B[i][j] && k + A[i][j] - B[i][j] <= 20000) {
        DP[i][j + 1][k + A[i][j] - B[i][j]] |= DP[i][j][k];
      }
      if (0 <= k + B[i][j] - A[i][j] && k + B[i][j] - A[i][j] <= 20000) {
        DP[i][j + 1][k + B[i][j] - A[i][j]] |= DP[i][j][k];
      }
    }
  }
  bool get[20001] = {};
  rep(i, 20000 + 1) {
    if (abs(i + A[H - 1][W - 1] - B[H - 1][W - 1] - 10000) <= 20000) {
      get[abs(i + A[H - 1][W - 1] - B[H - 1][W - 1] - 10000)] |=
          DP[H - 1][W - 1][i];
    }
    if (abs(i + B[H - 1][W - 1] - A[H - 1][W - 1] - 10000) <= 20000) {
      get[abs(i + B[H - 1][W - 1] - A[H - 1][W - 1] - 10000)] |=
          DP[H - 1][W - 1][i];
    }
  }
  ll res = INF;
  rep(i, 20000) if (get[i]) chmin(res, i);
  cout << res << rt;
}

// -g -D_GLIBCXX_DEBUG -fsanitize=undefined


Fixed submission:
#include <bits/stdc++.h>
using namespace std;

// Define
using ll = long long;
using ull = unsigned long long;
using ld = long double;
template <class T> using pvector = vector<pair<T, T>>;
template <class T>
using rpriority_queue = priority_queue<T, vector<T>, greater<T>>;
constexpr const ll dx[4] = {1, 0, -1, 0};
constexpr const ll dy[4] = {0, 1, 0, -1};
constexpr const ll MOD = 1e9 + 7;
constexpr const ll mod = 998244353;
constexpr const ll INF = 1LL << 60;
constexpr const ll inf = 1 << 30;
constexpr const char rt = '\n';
constexpr const char sp = ' ';
#define mp make_pair
#define mt make_tuple
#define pb push_back
#define eb emplace_back
#define elif else if
#define all(a, v, ...)                                                         \
  ([&](decltype((v)) w) { return (a)(begin(w), end(w), ##__VA_ARGS__); })(v)
#define fi first
#define se second

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 (a > b) {
    a = b;
    return 1;
  }
  return 0;
}

// Debug
#define debug(...)                                                             \
  {                                                                            \
    cerr << __LINE__ << ": " << #__VA_ARGS__ << " = ";                         \
    for (auto &&X : {__VA_ARGS__})                                             \
      cerr << "[" << X << "] ";                                                \
    cerr << rt;                                                                \
  }

#define dump(a, h, w)                                                          \
  {                                                                            \
    cerr << __LINE__ << ": " << #a << " = [" << rt;                            \
    rep(i, h) {                                                                \
      rep(j, w) cerr << a[i][j] << sp;                                         \
      cerr << rt;                                                              \
    }                                                                          \
    cerr << "]" << rt;                                                         \
  }

#define vdump(a, n)                                                            \
  {                                                                            \
    cerr << __LINE__ << ": " << #a << " = [";                                  \
    rep(i, n) if (i) cerr << sp << a[i];                                       \
    else cerr << a[i];                                                         \
    cerr << "]" << rt;                                                         \
  }

// Loop
#define inc(i, a, n) for (ll i = (a), _##i = (n); i <= _##i; ++i)
#define dec(i, a, n) for (ll i = (a), _##i = (n); i >= _##i; --i)
#define rep(i, n) for (ll i = 0, _##i = (n); i < _##i; ++i)
#define each(i, a) for (auto &&i : a)

// Stream
#define fout(n) cout << fixed << setprecision(n)
struct io {
  io() { cin.tie(nullptr), ios::sync_with_stdio(false); }
} io;

// Speed
#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

// Math
inline constexpr ll gcd(const ll a, const ll b) {
  return b ? gcd(b, a % b) : a;
}
inline constexpr ll lcm(const ll a, const ll b) { return a / gcd(a, b) * b; }

inline constexpr ll modulo(const ll n, const ll m = MOD) {
  ll k = n % m;
  return k + m * (k < 0);
}
inline constexpr ll chmod(ll &n, const ll m = MOD) {
  n %= m;
  return n += m * (n < 0);
}
inline constexpr ll mpow(ll a, ll n, const ll m = MOD) {
  ll r = 1;
  rep(i, 64) {
    if (n & (1LL << i))
      r *= a;
    chmod(r, m);
    a *= a;
    chmod(a, m);
  }
  return r;
}
inline ll inv(const ll n, const ll m = MOD) {
  ll a = n, b = m, x = 1, y = 0;
  while (b) {
    ll t = a / b;
    a -= t * b;
    swap(a, b);
    x -= t * y;
    swap(x, y);
  }
  return modulo(x, m);
}

bool DP[80][80][20000 + 1];
signed main() {
  ll H, W;
  cin >> H >> W;
  ll A[H][W], B[H][W];
  rep(i, H) rep(j, W) cin >> A[i][j];
  rep(i, H) rep(j, W) cin >> B[i][j];
  rep(i, H) rep(j, W) rep(k, 20000 + 1) DP[i][j][k] = 0;
  DP[0][0][10000] = 1;
  rep(i, H) rep(j, W) rep(k, 20000 + 1) {
    if (i + 1 <= H - 1) {
      if (0 <= k + A[i][j] - B[i][j] && k + A[i][j] - B[i][j] <= 20000) {
        DP[i + 1][j][k + A[i][j] - B[i][j]] |= DP[i][j][k];
      }
      if (0 <= k + B[i][j] - A[i][j] && k + B[i][j] - A[i][j] <= 20000) {
        DP[i + 1][j][k + B[i][j] - A[i][j]] |= DP[i][j][k];
      }
    }
    if (j + 1 <= W - 1) {
      if (0 <= k + A[i][j] - B[i][j] && k + A[i][j] - B[i][j] <= 20000) {
        DP[i][j + 1][k + A[i][j] - B[i][j]] |= DP[i][j][k];
      }
      if (0 <= k + B[i][j] - A[i][j] && k + B[i][j] - A[i][j] <= 20000) {
        DP[i][j + 1][k + B[i][j] - A[i][j]] |= DP[i][j][k];
      }
    }
  }
  bool get[20001] = {};
  rep(i, 20000 + 1) {
    if (abs(i + A[H - 1][W - 1] - B[H - 1][W - 1] - 10000) <= 20000) {
      get[abs(i + A[H - 1][W - 1] - B[H - 1][W - 1] - 10000)] |=
          DP[H - 1][W - 1][i];
    }
    if (abs(i + B[H - 1][W - 1] - A[H - 1][W - 1] - 10000) <= 20000) {
      get[abs(i + B[H - 1][W - 1] - A[H - 1][W - 1] - 10000)] |=
          DP[H - 1][W - 1][i];
    }
  }
  ll res = INF;
  rep(i, 20000) if (get[i]) chmin(res, i);
  cout << res << rt;
}

// -g -D_GLIBCXX_DEBUG -fsanitize=undefined