Problem p02782 (C++)

Buggy submission (status: Runtime Error):
#include <algorithm>
#include <cmath>
#include <cstring>
#include <functional>
#include <iomanip>
#include <iostream>
#include <limits>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>

using namespace std;

#ifdef LOCAL
#define eprintf(...) fprintf(stderr, __VA_ARGS__)
#else
#define eprintf(...) 42
#endif
using ll = long long;
using vll = std::vector<ll>;
using vvll = std::vector<vll>;
using vvvll = std::vector<vvll>;
#define reps(i, S, E) for (ll i = (S); i <= (E); i++)
#define rep(i, N) reps(i, 0, N - 1)
#define deps(i, E, S) for (ll i = (E); i >= (S); i--)
#define dep(i, N) deps(i, N - 1, 0)
const ll INF = 1LL << 60;
const int INF_INT = 1 << 30;

template <class T> inline bool chmax(T &a, T b) {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class T> inline bool chmin(T &a, T b) {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
struct mll {
  static ll MOD;
  ll val;
  mll(ll v = 0) : val(v % MOD) {
    if (val < 0)
      val += MOD;
  }
  mll operator-() const { return -val; }
  mll operator+(const mll &b) const { return val + b.val; }
  mll operator-(const mll &b) const { return val - b.val; }
  mll operator*(const mll &b) const { return val * b.val; }
  mll operator/(const mll &b) const { return mll(*this) /= b; }
  mll operator+(ll b) const { return *this + mll(b); }
  mll operator-(ll b) const { return *this - mll(b); }
  mll operator*(ll b) const { return *this * mll(b); }
  friend mll operator+(ll a, const mll &b) { return b + a; }
  friend mll operator-(ll a, const mll &b) { return -b + a; }
  friend mll operator*(ll a, const mll &b) { return b * a; }
  mll &operator+=(const mll &b) {
    val = (val + b.val) % MOD;
    return *this;
  }
  mll &operator-=(const mll &b) {
    val = (val + MOD - b.val) % MOD;
    return *this;
  }
  mll &operator*=(const mll &b) {
    val = (val * b.val) % MOD;
    return *this;
  }
  mll &operator/=(const mll &b) {
    ll c = b.val, d = MOD, u = 1, v = 0;
    while (d) {
      ll t = c / d;
      c -= t * d;
      swap(c, d);
      u -= t * v;
      swap(u, v);
    }
    val = val * u % MOD;
    if (val < 0)
      val += MOD;
    return *this;
  }
  mll &operator+=(ll b) { return *this += mll(b); }
  mll &operator-=(ll b) { return *this -= mll(b); }
  mll &operator*=(ll b) { return *this *= mll(b); }
  mll &operator/=(ll b) { return *this /= mll(b); }
  bool operator==(const mll &b) { return val == b.val; }
  bool operator!=(const mll &b) { return val != b.val; }
  bool operator==(ll b) { return *this == mll(b); }
  bool operator!=(ll b) { return *this != mll(b); }
  friend bool operator==(ll a, const mll &b) { return mll(a) == b.val; }
  friend bool operator!=(ll a, const mll &b) { return mll(a) != b.val; }
  friend ostream &operator<<(ostream &os, const mll &a) { return os << a.val; }
  friend istream &operator>>(istream &is, mll &a) { return is >> a.val; }
  static mll Combination(ll a, ll b) {
    chmin(b, a - b);
    if (b < 0)
      return mll(0);
    mll c = 1;
    rep(i, b) c *= a - i;
    rep(i, b) c /= i + 1;
    return c;
  }
};
using vmll = std::vector<mll>;
using vvmll = std::vector<vmll>;
using vvvmll = std::vector<vvmll>;
using vvvvmll = std::vector<vvvmll>;

struct Fast {
  Fast() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(std::numeric_limits<double>::max_digits10);
  }
} fast; // cin,cout高速化のおまじない＋桁数指定

ll mll::MOD = (ll)(1e9 + 7); // 998244353ll;

//--debug--
#ifdef LOCAL
void _debug_out_(string s = "") {
  if (s != "") {
    eprintf(s.c_str());
  } else {
    eprintf("\n");
  }
}
void _debug_out_(ll x, string s = "") {
  if (s != "") {
    eprintf(s.c_str(), x);
  } else {
    eprintf("%lld\n", x);
  }
}
void _debug_out_(vvll x, string s = "") {
  string t = s;
  if (s == "") {
    ll maxLen = 0;
    rep(i, (ll)x.size()) rep(j, (ll)x[i].size()) {
      ll n = x[i][j] / 10, len = 1;
      while (n > 0) {
        len++;
        n /= 10;
      }
      chmax(maxLen, len);
    }
    t = "%" + to_string(maxLen) + "lld ";
  }
  rep(i, (ll)x.size()) {
    rep(j, (ll)x[i].size()) { eprintf(t.c_str(), x[i][j]); }
    eprintf("\n");
  }
  eprintf("\n");
}
void _debug_out_(vll x, string s = "") {
  rep(i, (ll)x.size()) {
    if (s != "") {
      eprintf(s.c_str(), x[i]);
    } else {
      eprintf("%lld ", x[i]);
    }
  }
  eprintf("\n");
}
void _debug_out_(vector<bool> x) {
  rep(i, (ll)x.size()) { eprintf("%s", x[i] ? "o" : "x"); }
  eprintf("\n");
}
void _debug_out_(vector<vector<bool>> x) {
  rep(i, (ll)x.size()) {
    rep(j, (ll)x[i].size()) { eprintf("%s", x[i][j] ? "o" : "x"); }
    eprintf("\n");
  }
  eprintf("\n");
}
void _debug_out_(map<ll, ll> x) {
  auto begin = x.begin(), end = x.end();
  for (auto iter = begin; iter != end; iter++) {
    eprintf("%lld:%lld ", iter->first, iter->second);
  }
  eprintf("\n");
}
void _debug_out_(map<char, ll> x) {
  auto begin = x.begin(), end = x.end();
  for (auto iter = begin; iter != end; iter++) {
    eprintf("%c:%lld ", iter->first, iter->second);
  }
  eprintf("\n");
}
#define debug_out(...) _debug_out_(__VA_ARGS__)
#else
#define debug_out(...) 42
#endif
//--debug--

const ll MAX = 1000010;
const ll MOD = 1000000007;
ll fac[MAX], finv[MAX], inv[MAX];

// テーブルを作る前処理
void COMinit() {
  fac[0] = fac[1] = 1;
  finv[0] = finv[1] = 1;
  inv[1] = 1;
  for (ll i = 2; i < MAX; i++) {
    fac[i] = fac[i - 1] * i % MOD;
    inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD;
    finv[i] = finv[i - 1] * inv[i] % MOD;
  }
}

// 二項係数計算
mll COM(ll n, ll k) {
  if (n < k)
    return 0;
  if (n < 0 || k < 0)
    return 0;
  return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;
}

mll f(ll r, ll c) { //(r,c)へのパス数
  return COM(r + c, r);
}

int main() {
  ll r1;
  cin >> r1;
  ll c1;
  cin >> c1;
  ll r2;
  cin >> r2;
  ll c2;
  cin >> c2;

  COMinit();

  mll result = f(r2 + 1, c2 + 1);
  result -= f(r1 - 1 + 1, c2 + 1);
  result -= f(r2 + 1, c1 - 1 + 1);
  result += f(r1 - 1 + 1, c1 - 1 + 1);

  cout << result << endl;

  return 0;
}

Fixed submission:
#include <algorithm>
#include <cmath>
#include <cstring>
#include <functional>
#include <iomanip>
#include <iostream>
#include <limits>
#include <map>
#include <queue>
#include <set>
#include <stack>
#include <string>
#include <vector>

using namespace std;

#ifdef LOCAL
#define eprintf(...) fprintf(stderr, __VA_ARGS__)
#else
#define eprintf(...) 42
#endif
using ll = long long;
using vll = std::vector<ll>;
using vvll = std::vector<vll>;
using vvvll = std::vector<vvll>;
#define reps(i, S, E) for (ll i = (S); i <= (E); i++)
#define rep(i, N) reps(i, 0, N - 1)
#define deps(i, E, S) for (ll i = (E); i >= (S); i--)
#define dep(i, N) deps(i, N - 1, 0)
const ll INF = 1LL << 60;
const int INF_INT = 1 << 30;

template <class T> inline bool chmax(T &a, T b) {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class T> inline bool chmin(T &a, T b) {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
struct mll {
  static ll MOD;
  ll val;
  mll(ll v = 0) : val(v % MOD) {
    if (val < 0)
      val += MOD;
  }
  mll operator-() const { return -val; }
  mll operator+(const mll &b) const { return val + b.val; }
  mll operator-(const mll &b) const { return val - b.val; }
  mll operator*(const mll &b) const { return val * b.val; }
  mll operator/(const mll &b) const { return mll(*this) /= b; }
  mll operator+(ll b) const { return *this + mll(b); }
  mll operator-(ll b) const { return *this - mll(b); }
  mll operator*(ll b) const { return *this * mll(b); }
  friend mll operator+(ll a, const mll &b) { return b + a; }
  friend mll operator-(ll a, const mll &b) { return -b + a; }
  friend mll operator*(ll a, const mll &b) { return b * a; }
  mll &operator+=(const mll &b) {
    val = (val + b.val) % MOD;
    return *this;
  }
  mll &operator-=(const mll &b) {
    val = (val + MOD - b.val) % MOD;
    return *this;
  }
  mll &operator*=(const mll &b) {
    val = (val * b.val) % MOD;
    return *this;
  }
  mll &operator/=(const mll &b) {
    ll c = b.val, d = MOD, u = 1, v = 0;
    while (d) {
      ll t = c / d;
      c -= t * d;
      swap(c, d);
      u -= t * v;
      swap(u, v);
    }
    val = val * u % MOD;
    if (val < 0)
      val += MOD;
    return *this;
  }
  mll &operator+=(ll b) { return *this += mll(b); }
  mll &operator-=(ll b) { return *this -= mll(b); }
  mll &operator*=(ll b) { return *this *= mll(b); }
  mll &operator/=(ll b) { return *this /= mll(b); }
  bool operator==(const mll &b) { return val == b.val; }
  bool operator!=(const mll &b) { return val != b.val; }
  bool operator==(ll b) { return *this == mll(b); }
  bool operator!=(ll b) { return *this != mll(b); }
  friend bool operator==(ll a, const mll &b) { return mll(a) == b.val; }
  friend bool operator!=(ll a, const mll &b) { return mll(a) != b.val; }
  friend ostream &operator<<(ostream &os, const mll &a) { return os << a.val; }
  friend istream &operator>>(istream &is, mll &a) { return is >> a.val; }
  static mll Combination(ll a, ll b) {
    chmin(b, a - b);
    if (b < 0)
      return mll(0);
    mll c = 1;
    rep(i, b) c *= a - i;
    rep(i, b) c /= i + 1;
    return c;
  }
};
using vmll = std::vector<mll>;
using vvmll = std::vector<vmll>;
using vvvmll = std::vector<vvmll>;
using vvvvmll = std::vector<vvvmll>;

struct Fast {
  Fast() {
    cin.tie(0);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(std::numeric_limits<double>::max_digits10);
  }
} fast; // cin,cout高速化のおまじない＋桁数指定

ll mll::MOD = (ll)(1e9 + 7); // 998244353ll;

//--debug--
#ifdef LOCAL
void _debug_out_(string s = "") {
  if (s != "") {
    eprintf(s.c_str());
  } else {
    eprintf("\n");
  }
}
void _debug_out_(ll x, string s = "") {
  if (s != "") {
    eprintf(s.c_str(), x);
  } else {
    eprintf("%lld\n", x);
  }
}
void _debug_out_(vvll x, string s = "") {
  string t = s;
  if (s == "") {
    ll maxLen = 0;
    rep(i, (ll)x.size()) rep(j, (ll)x[i].size()) {
      ll n = x[i][j] / 10, len = 1;
      while (n > 0) {
        len++;
        n /= 10;
      }
      chmax(maxLen, len);
    }
    t = "%" + to_string(maxLen) + "lld ";
  }
  rep(i, (ll)x.size()) {
    rep(j, (ll)x[i].size()) { eprintf(t.c_str(), x[i][j]); }
    eprintf("\n");
  }
  eprintf("\n");
}
void _debug_out_(vll x, string s = "") {
  rep(i, (ll)x.size()) {
    if (s != "") {
      eprintf(s.c_str(), x[i]);
    } else {
      eprintf("%lld ", x[i]);
    }
  }
  eprintf("\n");
}
void _debug_out_(vector<bool> x) {
  rep(i, (ll)x.size()) { eprintf("%s", x[i] ? "o" : "x"); }
  eprintf("\n");
}
void _debug_out_(vector<vector<bool>> x) {
  rep(i, (ll)x.size()) {
    rep(j, (ll)x[i].size()) { eprintf("%s", x[i][j] ? "o" : "x"); }
    eprintf("\n");
  }
  eprintf("\n");
}
void _debug_out_(map<ll, ll> x) {
  auto begin = x.begin(), end = x.end();
  for (auto iter = begin; iter != end; iter++) {
    eprintf("%lld:%lld ", iter->first, iter->second);
  }
  eprintf("\n");
}
void _debug_out_(map<char, ll> x) {
  auto begin = x.begin(), end = x.end();
  for (auto iter = begin; iter != end; iter++) {
    eprintf("%c:%lld ", iter->first, iter->second);
  }
  eprintf("\n");
}
#define debug_out(...) _debug_out_(__VA_ARGS__)
#else
#define debug_out(...) 42
#endif
//--debug--

const ll MAX = 2000010;
const ll MOD = 1000000007;
ll fac[MAX], finv[MAX], inv[MAX];

// テーブルを作る前処理
void COMinit() {
  fac[0] = fac[1] = 1;
  finv[0] = finv[1] = 1;
  inv[1] = 1;
  for (ll i = 2; i < MAX; i++) {
    fac[i] = fac[i - 1] * i % MOD;
    inv[i] = MOD - inv[MOD % i] * (MOD / i) % MOD;
    finv[i] = finv[i - 1] * inv[i] % MOD;
  }
}

// 二項係数計算
mll COM(ll n, ll k) {
  if (n < k)
    return 0;
  if (n < 0 || k < 0)
    return 0;
  return fac[n] * (finv[k] * finv[n - k] % MOD) % MOD;
}

mll f(ll r, ll c) { //(r,c)へのパス数
  return COM(r + c, r);
}

int main() {
  ll r1;
  cin >> r1;
  ll c1;
  cin >> c1;
  ll r2;
  cin >> r2;
  ll c2;
  cin >> c2;

  COMinit();

  mll result = f(r2 + 1, c2 + 1);
  result -= f(r1 - 1 + 1, c2 + 1);
  result -= f(r2 + 1, c1 - 1 + 1);
  result += f(r1 - 1 + 1, c1 - 1 + 1);

  cout << result << endl;

  return 0;
}