Problem p02782 (C++)

Buggy submission (status: Runtime Error):
#include <bits/stdc++.h>
#include <stdlib.h>
#include <sys/time.h>
#include <unistd.h>

using namespace std;

typedef long long ll;
typedef pair<ll, ll> P;
typedef pair<P, ll> T;

long long int INF = 2e18;
const ll fact_table = 1200008;
double Pi = 3.1415926535897932384626;

priority_queue<ll> pql;
priority_queue<P> pqp;
priority_queue<P> bag;
// big priority queue
priority_queue<ll, vector<ll>, greater<ll>> pqls;
priority_queue<P, vector<P>, greater<P>> pqps;
// small priority queue
// top pop

ll dx[8] = {1, 0, -1, 0, 1, 1, -1, -1};
ll dy[8] = {0, 1, 0, -1, 1, -1, -1, 1};
char dir[] = "DRUL";
// ↓,→,↑,←

#define endl "\n"
#ifdef ENJAPMA
#undef endl
#endif
#define p(x) cout << x << endl;
#define el cout << endl;
#define pe(x) cout << x << " ";
#define ps(x) cout << fixed << setprecision(25) << x << endl;
#define pu(x) cout << (x);
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define CLEAR(a) a = decltype(a)();
#define pc(x) cout << x << ",";

// const ll mod = 998244353ll;
const ll mod = 1000000007ll;

ll mypow(ll number1, ll number2, ll mod) {
  if (number2 == 0) {
    return 1ll;
  } else {
    ll number3 = mypow(number1, number2 / 2, mod);
    number3 *= number3;
    number3 %= mod;
    if (number2 % 2 == 1) {
      number3 *= number1;
      number3 %= mod;
    }
    return number3;
  }
}
void YES(bool condition) {
  if (condition) {
    p("YES");
  } else {
    p("NO");
  }
  return;
}
void Yes(bool condition) {
  if (condition) {
    p("Yes");
  } else {
    p("No");
  }
  return;
}

ll fact[fact_table + 5], rfact[fact_table + 5];

void c3_init() {
  fact[0] = rfact[0] = 1;
  for (ll i = 1; i <= fact_table; i++) {
    fact[i] = (fact[i - 1] * i) % mod;
  }
  rfact[fact_table] = mypow(fact[fact_table], mod - 2, mod);
  for (ll i = fact_table; i >= 1; i--) {
    rfact[i - 1] = rfact[i] * i;
    rfact[i - 1] %= mod;
  }
  return;
}
ll c3(ll n, ll r) {
  return (((fact[n] * rfact[r]) % mod) * rfact[n - r]) % mod;
}

bool multicase = false;

ll n, m, num, a, b, c, d, e, h, q;
ll w, k, ans;
ll x[500005], y[500005];
ll dp[15][500005];
string s;
ll query = 0;

ll calc(ll l, ll r) {
  ll res = c3(r + l + 2, r + 1);
  res += mod - 1;
  res %= mod;
  return res;
}

void solve() {
  c3_init();
  ll r1, r2, c1, c2;
  cin >> r1 >> c1 >> r2 >> c2;
  ll ans = 0;
  ll a = calc(r2, c2);
  ans =
      calc(r2, c2) - calc(r1 - 1, c2) - calc(r2, c1 - 1) + calc(r1 - 1, c1 - 1);

  ans += mod * 10;
  ans %= mod;
  p(ans);
  return;
}

int main() {
  // init();
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  ll q, testcase = 1;
  if (multicase) {
    cin >> q;
  } else {
    q = 1;
  }
  while (q--) {
    // pu("Case ");pu("#");pu(testcase);pu(": ");
    solve();
    testcase++;
  }
  // solve();
  return 0;
}


Fixed submission:
#include <bits/stdc++.h>
#include <stdlib.h>
#include <sys/time.h>
#include <unistd.h>

using namespace std;

typedef long long ll;
typedef pair<ll, ll> P;
typedef pair<P, ll> T;

long long int INF = 2e18;
const ll fact_table = 3200008;
double Pi = 3.1415926535897932384626;

priority_queue<ll> pql;
priority_queue<P> pqp;
priority_queue<P> bag;
// big priority queue
priority_queue<ll, vector<ll>, greater<ll>> pqls;
priority_queue<P, vector<P>, greater<P>> pqps;
// small priority queue
// top pop

ll dx[8] = {1, 0, -1, 0, 1, 1, -1, -1};
ll dy[8] = {0, 1, 0, -1, 1, -1, -1, 1};
char dir[] = "DRUL";
// ↓,→,↑,←

#define endl "\n"
#ifdef ENJAPMA
#undef endl
#endif
#define p(x) cout << x << endl;
#define el cout << endl;
#define pe(x) cout << x << " ";
#define ps(x) cout << fixed << setprecision(25) << x << endl;
#define pu(x) cout << (x);
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define CLEAR(a) a = decltype(a)();
#define pc(x) cout << x << ",";

// const ll mod = 998244353ll;
const ll mod = 1000000007ll;

ll mypow(ll number1, ll number2, ll mod) {
  if (number2 == 0) {
    return 1ll;
  } else {
    ll number3 = mypow(number1, number2 / 2, mod);
    number3 *= number3;
    number3 %= mod;
    if (number2 % 2 == 1) {
      number3 *= number1;
      number3 %= mod;
    }
    return number3;
  }
}
void YES(bool condition) {
  if (condition) {
    p("YES");
  } else {
    p("NO");
  }
  return;
}
void Yes(bool condition) {
  if (condition) {
    p("Yes");
  } else {
    p("No");
  }
  return;
}

ll fact[fact_table + 5], rfact[fact_table + 5];

void c3_init() {
  fact[0] = rfact[0] = 1;
  for (ll i = 1; i <= fact_table; i++) {
    fact[i] = (fact[i - 1] * i) % mod;
  }
  rfact[fact_table] = mypow(fact[fact_table], mod - 2, mod);
  for (ll i = fact_table; i >= 1; i--) {
    rfact[i - 1] = rfact[i] * i;
    rfact[i - 1] %= mod;
  }
  return;
}
ll c3(ll n, ll r) {
  return (((fact[n] * rfact[r]) % mod) * rfact[n - r]) % mod;
}

bool multicase = false;

ll n, m, num, a, b, c, d, e, h, q;
ll w, k, ans;
ll x[500005], y[500005];
ll dp[15][500005];
string s;
ll query = 0;

ll calc(ll l, ll r) {
  ll res = c3(r + l + 2, r + 1);
  res += mod - 1;
  res %= mod;
  return res;
}

void solve() {
  c3_init();
  ll r1, r2, c1, c2;
  cin >> r1 >> c1 >> r2 >> c2;
  ll ans = 0;
  ll a = calc(r2, c2);
  ans =
      calc(r2, c2) - calc(r1 - 1, c2) - calc(r2, c1 - 1) + calc(r1 - 1, c1 - 1);

  ans += mod * 10;
  ans %= mod;
  p(ans);
  return;
}

int main() {
  // init();
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  ll q, testcase = 1;
  if (multicase) {
    cin >> q;
  } else {
    q = 1;
  }
  while (q--) {
    // pu("Case ");pu("#");pu(testcase);pu(": ");
    solve();
    testcase++;
  }
  // solve();
  return 0;
}