Problem p02782 (C++)

Buggy submission (status: Runtime Error):
#define _CRT_SECURE_NO_WARNINGS
#define _USE_MATH_DEFINES
#include "bits/stdc++.h"
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define FOR(i, m, n) for (int i = (m); i < (n); ++i)
#define rrep(i, n) for (int i = (n)-1; i >= 0; --i)
#define rfor(i, m, n) for (int i = (m); i >= (n); --i)
#define sz(x) ((int)(x).size())
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define mp make_pair
#define pb push_back
#define eb emplace_back

using namespace std;
using LL = long long;
using ld = long double;
using VB = vector<bool>;
using VVB = vector<VB>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<LL>;
using VVL = vector<VL>;
using VS = vector<string>;
using VD = vector<double>;
using PII = pair<int, int>;
using VP = vector<PII>;
using PLL = pair<LL, LL>;
using VPL = vector<PLL>;
template <class T> using Grid = vector<vector<T>>;
template <class T> using PQ = priority_queue<T>;
template <class T> using PQS = priority_queue<T, vector<T>, greater<T>>;
constexpr int inf = (int)1e9;
constexpr LL inf_ll = (LL)1e18, MOD = 1000000007;
constexpr ld PI = M_PI, EPS = 1e-12;

template <class T> inline void Sort(T &a) noexcept { sort(all(a)); }
template <class T> inline void RSort(T &a) noexcept { sort(rall(a)); }
template <class T> inline void Reverse(T &a) noexcept { reverse(all(a)); }
template <class T> inline void Unique(T &a) noexcept {
  a.erase(unique(all(a)), a.end());
}
template <class T> inline T Sorted(T a) noexcept {
  Sort(a);
  return a;
}
template <class T> inline T RSorted(T a) noexcept {
  RSort(a);
  return a;
}
template <class T> inline T Reversed(T a) noexcept {
  Reverse(a);
  return a;
}
template <class T> inline T Uniqued(T a) noexcept {
  Unique(a);
  return a;
}
template <class T> inline auto Max(const T &a) noexcept {
  return *max_element(all(a));
}
template <class T> inline auto Min(const T &a) noexcept {
  return *min_element(all(a));
}
template <class T> inline int MaxPos(const T &a) noexcept {
  return max_element(all(a)) - a.begin();
}
template <class T> inline int MinPos(const T &a) noexcept {
  return min_element(all(a)) - a.begin();
}
template <class T, class U> inline int Count(const T &a, const U &v) noexcept {
  return count(all(a), v);
}
template <class T, class U> inline int Find(const T &a, const U &v) noexcept {
  auto pos = find(all(a), v);
  return pos == a.end() ? -1 : pos - a.begin();
}
template <class T, class U> inline U Sum(const T &a, const U &v) noexcept {
  return accumulate(all(a), v);
}
template <class T, class U> inline int Lower(const T &a, const U &v) noexcept {
  return lower_bound(all(a), v) - a.begin();
}
template <class T, class U> inline int Upper(const T &a, const U &v) noexcept {
  return upper_bound(all(a), v) - a.begin();
}
template <class T, class P> inline void RemoveIf(T &a, P f) noexcept {
  a.erase(remove_if(all(a), f), a.end());
}
template <class T> inline T Age(T n, T m) noexcept { return (n + m - 1) / m; }
template <class T> inline T Age2(T n, T m) noexcept { return Age(n, m) * m; }
template <class T> inline T Tri(T n) noexcept {
  return (n & 1) ? (n + 1) / 2 * n : n / 2 * (n + 1);
}
template <class T = long long> inline T BIT(int b) noexcept {
  return T{1} << b;
}
template <class T> inline T Gcd(T n, T m) noexcept {
  return m ? Gcd(m, n % m) : n;
}
template <class T> inline T Lcm(T n, T m) noexcept { return n / Gcd(n, m) * m; }
template <class T> inline T Pow(T a, T n) noexcept {
  T r = 1;
  while (n > 0) {
    if (n & 1)
      r *= a;
    a *= a;
    n /= 2;
  }
  return r;
}
template <class T> inline T Powmod(T a, T n, T m = MOD) noexcept {
  T r = 1;
  while (n > 0) {
    if (n & 1)
      r = r * a % m, n--;
    else
      a = a * a % m, n /= 2;
  }
  return r;
}
template <class T> inline bool chmax(T &a, const T &b) noexcept {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class T> inline bool chmin(T &a, const T &b) noexcept {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
template <class T>
inline bool inRange(const T &v, const T &min, const T &max) noexcept {
  return min <= v && v <= max;
}
inline string operator*(string s, int n) noexcept {
  string ret;
  rep(i, n) ret += s;
  return ret;
}

// --- input --- //
#if defined(_WIN32) || defined(ONLINE_JUDGE)
#define getchar_unlocked _getchar_nolock
#define putchar_unlocked _putchar_nolock
#endif
inline int gc() noexcept { return getchar_unlocked(); }
template <class T> inline void InputF(T &v) noexcept { cin >> v; }
inline void InputF(char &v) noexcept {
  while (isspace(v = gc()))
    ;
}
inline void InputF(string &v) noexcept {
  char c;
  for (InputF(c); !isspace(c); c = gc())
    v += c;
}
inline void InputF(int &v) noexcept {
  bool neg = false;
  v = 0;
  char c;
  InputF(c);
  if (c == '-') {
    neg = true;
    c = gc();
  }
  for (; isdigit(c); c = gc())
    v = v * 10 + (c - '0');
  if (neg)
    v = -v;
}
inline void InputF(long long &v) noexcept {
  bool neg = false;
  v = 0;
  char c;
  InputF(c);
  if (c == '-') {
    neg = true;
    c = gc();
  }
  for (; isdigit(c); c = gc())
    v = v * 10 + (LL)(c - '0');
  if (neg)
    v = -v;
}
inline void InputF(double &v) noexcept {
  double dp = 1;
  bool neg = false, adp = false;
  v = 0;
  char c;
  InputF(c);
  if (c == '-') {
    neg = true;
    c = gc();
  }
  for (; isdigit(c) || c == '.'; c = gc()) {
    if (c == '.')
      adp = true;
    else if (adp)
      v += (c - '0') * (dp *= 0.1);
    else
      v = v * 10 + (c - '0');
  }
  if (neg)
    v = -v;
}
template <class T, class U> inline void InputF(pair<T, U> &v) noexcept {
  InputF(v.first);
  InputF(v.second);
}
template <class T> inline void InputF(vector<T> &v) noexcept {
  for (auto &e : v)
    InputF(e);
}
template <class T> inline T input() {
  T v;
  InputF(v);
  return v;
}
struct InputV {
  int n, m;
  InputV(int N) : n(N), m(0) {}
  InputV(pair<int, int> N) : n(N.first), m(N.second) {}
  template <class T> operator vector<T>() noexcept {
    vector<T> v(n);
    InputF(v);
    return v;
  }
  template <class T> operator vector<vector<T>>() noexcept {
    vector<vector<T>> v(n, vector<T>(m));
    InputF(v);
    return v;
  }
};
struct Input {
  template <class T> operator T() noexcept { return input<T>(); }
  int operator--(int) {
    int v;
    InputF(v);
    v--;
    return v;
  }
  InputV operator[](int n) noexcept { return InputV(n); }
  InputV operator[](pair<int, int> n) noexcept { return InputV(n); }
  void operator()() {}
  template <class H, class... T> void operator()(H &&h, T &&...t) {
    InputF(h);
    operator()(forward<T>(t)...);
  }
  template <class T, size_t W> array<vector<T>, W> get(int H) {
    array<vector<T>, W> ret;
    rep(i, H) rep(j, W) {
      T x = *this;
      ret[j].push_back(x);
    }
    return ret;
  }
} in;

// --- output --- //
struct BoolStr {
  const char *t, *f;
  BoolStr(const char *_t, const char *_f) : t(_t), f(_f) {}
} Yes("Yes", "No"), yes("yes", "no"), YES("YES", "NO"), Int("1", "0");
struct DivStr {
  const char *d, *l;
  DivStr(const char *_d, const char *_l) : d(_d), l(_l) {}
} spc(" ", "\n"), no_spc("", "\n"), end_line("\n", "\n"), comma(",", "\n"),
    no_endl(" ", "");
class Output {
  BoolStr B{Yes};
  DivStr D{spc};
  bool isPrint = true;
  void p(double v) { printf("%.20f", v); }
  void p(long double v) { printf("%.20Lf", v); }
  void p(int v) { printf("%d", v); }
  void p(LL v) { printf("%lld", v); }
  void p(char v) { putchar(v); }
  void p(bool v) { printf("%s", v ? B.t : B.f); }
  template <class T> void p(const T &v) { cout << v; }
  template <class T, class U> void p(const pair<T, U> &v) {
    p(v.first);
    printf("%s", D.d);
    p(v.second);
  }
  template <class T> void p(const vector<T> &v) {
    rep(i, sz(v)) {
      if (i)
        printf("%s", D.d);
      p(v[i]);
    }
  }
  template <class T> void p(const vector<vector<T>> &v) {
    rep(i, sz(v)) {
      if (i)
        printf("%s", D.l);
      p(v[i]);
    }
  }
  void p(const BoolStr &v) {
    B = v;
    isPrint = false;
  }
  void p(const DivStr &v) {
    D = v;
    isPrint = false;
  }

public:
  void operator()() { printf("%s", D.l); }
  template <class H> void operator()(H &&h) {
    p(h);
    if (isPrint)
      printf("%s", D.l);
    isPrint = true;
    B = Yes;
    D = spc;
  }
  template <class H, class... T> void operator()(H &&h, T &&...t) {
    p(h);
    if (isPrint)
      printf("%s", D.d);
    isPrint = true;
    operator()(forward<T>(t)...);
  }
  template <class... T> void exit(T &&...t) {
    operator()(forward<T>(t)...);
    std::exit(EXIT_SUCCESS);
  }
} out;

// --- dump --- //
#if __has_include("dump.hpp")
#include "dump.hpp"
#else
#define dump(...) (void(0))
#endif

// ---------------------------------------------------------------- //

template <class T> class Combi {
  int Max;
  T Mod;
  vector<T> fac, finv, inv;

public:
  Combi(int _Max = 300009, T _Mod = 1000000007)
      : Max(_Max), Mod(_Mod), fac(_Max), finv(_Max), inv(_Max) {
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (int 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;
    }
  }
  T P(int n, int r) { // n P r
    return (n < r || n < 0 || r < 0) ? 0 : fac[n] * finv[n - r] % MOD;
  }
  T C(int n, int r) { // n C r
    return (n < r || n < 0 || r < 0)
               ? 0
               : fac[n] * (finv[r] * finv[n - r] % Mod) % Mod;
  }
  T H(int n, int r) { // n H r
    return (n < 0 || r < 0) ? 0 : r == 0 ? 1 : C(n + r - 1, r);
  }
};

template <int MOD> struct modint {
  using T = long long;
  T n;
  constexpr modint(const T x = 0) noexcept : n(x % MOD) {
    if (n < 0)
      n += MOD;
  }
  constexpr int get_mod() const noexcept { return MOD; }
  constexpr modint operator+() const noexcept { return *this; }
  constexpr modint operator-() const noexcept { return n ? MOD - n : 0; }
  constexpr modint &operator++() noexcept {
    if (MOD <= ++n)
      n = 0;
    return *this;
  }
  constexpr modint &operator--() noexcept {
    if (n <= 0)
      n = MOD;
    n--;
    return *this;
  }
  constexpr modint &operator++(int) noexcept {
    modint t = *this;
    ++*this;
    return t;
  }
  constexpr modint &operator--(int) noexcept {
    modint t = *this;
    --*this;
    return t;
  }
  constexpr modint next() const noexcept {
    modint t = *this;
    ++t;
    return t;
  }
  constexpr modint pred() const noexcept {
    modint t = *this;
    --t;
    return t;
  }
  constexpr modint operator+(const modint &m) const noexcept {
    return modint(*this) += m;
  }
  constexpr modint operator-(const modint &m) const noexcept {
    return modint(*this) -= m;
  }
  constexpr modint operator*(const modint &m) const noexcept {
    return modint(*this) *= m;
  }
  constexpr modint operator/(const modint &m) const noexcept {
    return modint(*this) /= m;
  }
  constexpr modint &operator+=(const modint &m) noexcept {
    n += m.n;
    if (n >= MOD)
      n -= MOD;
    return *this;
  }
  constexpr modint &operator-=(const modint &m) noexcept {
    n -= m.n;
    if (n < 0)
      n += MOD;
    return *this;
  }
  constexpr modint &operator*=(const modint &m) noexcept {
    n = n * m.n % MOD;
    return *this;
  }
  constexpr modint &operator/=(const modint &m) noexcept {
    assert(m.n);
    T a = m.n, b = MOD, u = 1, v = 0;
    while (b) {
      T t = a / b;
      a -= t * b;
      swap(a, b);
      u -= t * v;
      swap(u, v);
    }
    n = n * u % MOD;
    if (n < 0)
      n += MOD;
    return *this;
  }
  constexpr bool operator==(const modint &m) const noexcept { return n == m.n; }
  constexpr bool operator!=(const modint &m) const noexcept { return n != m.n; }
  constexpr modint pow(modint m) const noexcept {
    modint t = n, res = 1;
    while (m.n > 0) {
      if (m.n & 1)
        res *= t;
      t *= t;
      m.n >>= 1;
    }
    return res;
  }
  constexpr modint operator^(modint m) const noexcept { return pow(m); }
};
using mint = modint<1000000007>;
using VM = vector<mint>;
ostream &operator<<(ostream &os, const modint<MOD> &m) noexcept {
  return os << m.n;
}
istream &operator>>(istream &is, modint<MOD> &m) noexcept { return is >> m.n; }
mint operator""_m(unsigned long long n) { return n; }

int main() {
  int x1, y1, x2, y2;
  in(x1, y1, x2, y2);
  Combi<LL> comb(1000007);

  auto calc = [&](int x1, int y1, int x2, int y2) {
    mint ans;
    FOR(x, x1, x2 + 1) {
      mint val1 = mint(comb.C(x + y2, y2)) * 2;
      mint val2 = mint(comb.C(x + y1 - 1, y1 - 1)) * 2;
      ans += val1 - val2;
    }
    return ans;
  };

  out(calc(x1, y1, x2, y2));
}

Fixed submission:
#define _CRT_SECURE_NO_WARNINGS
#define _USE_MATH_DEFINES
#include "bits/stdc++.h"
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define FOR(i, m, n) for (int i = (m); i < (n); ++i)
#define rrep(i, n) for (int i = (n)-1; i >= 0; --i)
#define rfor(i, m, n) for (int i = (m); i >= (n); --i)
#define sz(x) ((int)(x).size())
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define mp make_pair
#define pb push_back
#define eb emplace_back

using namespace std;
using LL = long long;
using ld = long double;
using VB = vector<bool>;
using VVB = vector<VB>;
using VI = vector<int>;
using VVI = vector<VI>;
using VL = vector<LL>;
using VVL = vector<VL>;
using VS = vector<string>;
using VD = vector<double>;
using PII = pair<int, int>;
using VP = vector<PII>;
using PLL = pair<LL, LL>;
using VPL = vector<PLL>;
template <class T> using Grid = vector<vector<T>>;
template <class T> using PQ = priority_queue<T>;
template <class T> using PQS = priority_queue<T, vector<T>, greater<T>>;
constexpr int inf = (int)1e9;
constexpr LL inf_ll = (LL)1e18, MOD = 1000000007;
constexpr ld PI = M_PI, EPS = 1e-12;

template <class T> inline void Sort(T &a) noexcept { sort(all(a)); }
template <class T> inline void RSort(T &a) noexcept { sort(rall(a)); }
template <class T> inline void Reverse(T &a) noexcept { reverse(all(a)); }
template <class T> inline void Unique(T &a) noexcept {
  a.erase(unique(all(a)), a.end());
}
template <class T> inline T Sorted(T a) noexcept {
  Sort(a);
  return a;
}
template <class T> inline T RSorted(T a) noexcept {
  RSort(a);
  return a;
}
template <class T> inline T Reversed(T a) noexcept {
  Reverse(a);
  return a;
}
template <class T> inline T Uniqued(T a) noexcept {
  Unique(a);
  return a;
}
template <class T> inline auto Max(const T &a) noexcept {
  return *max_element(all(a));
}
template <class T> inline auto Min(const T &a) noexcept {
  return *min_element(all(a));
}
template <class T> inline int MaxPos(const T &a) noexcept {
  return max_element(all(a)) - a.begin();
}
template <class T> inline int MinPos(const T &a) noexcept {
  return min_element(all(a)) - a.begin();
}
template <class T, class U> inline int Count(const T &a, const U &v) noexcept {
  return count(all(a), v);
}
template <class T, class U> inline int Find(const T &a, const U &v) noexcept {
  auto pos = find(all(a), v);
  return pos == a.end() ? -1 : pos - a.begin();
}
template <class T, class U> inline U Sum(const T &a, const U &v) noexcept {
  return accumulate(all(a), v);
}
template <class T, class U> inline int Lower(const T &a, const U &v) noexcept {
  return lower_bound(all(a), v) - a.begin();
}
template <class T, class U> inline int Upper(const T &a, const U &v) noexcept {
  return upper_bound(all(a), v) - a.begin();
}
template <class T, class P> inline void RemoveIf(T &a, P f) noexcept {
  a.erase(remove_if(all(a), f), a.end());
}
template <class T> inline T Age(T n, T m) noexcept { return (n + m - 1) / m; }
template <class T> inline T Age2(T n, T m) noexcept { return Age(n, m) * m; }
template <class T> inline T Tri(T n) noexcept {
  return (n & 1) ? (n + 1) / 2 * n : n / 2 * (n + 1);
}
template <class T = long long> inline T BIT(int b) noexcept {
  return T{1} << b;
}
template <class T> inline T Gcd(T n, T m) noexcept {
  return m ? Gcd(m, n % m) : n;
}
template <class T> inline T Lcm(T n, T m) noexcept { return n / Gcd(n, m) * m; }
template <class T> inline T Pow(T a, T n) noexcept {
  T r = 1;
  while (n > 0) {
    if (n & 1)
      r *= a;
    a *= a;
    n /= 2;
  }
  return r;
}
template <class T> inline T Powmod(T a, T n, T m = MOD) noexcept {
  T r = 1;
  while (n > 0) {
    if (n & 1)
      r = r * a % m, n--;
    else
      a = a * a % m, n /= 2;
  }
  return r;
}
template <class T> inline bool chmax(T &a, const T &b) noexcept {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class T> inline bool chmin(T &a, const T &b) noexcept {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
template <class T>
inline bool inRange(const T &v, const T &min, const T &max) noexcept {
  return min <= v && v <= max;
}
inline string operator*(string s, int n) noexcept {
  string ret;
  rep(i, n) ret += s;
  return ret;
}

// --- input --- //
#if defined(_WIN32) || defined(ONLINE_JUDGE)
#define getchar_unlocked _getchar_nolock
#define putchar_unlocked _putchar_nolock
#endif
inline int gc() noexcept { return getchar_unlocked(); }
template <class T> inline void InputF(T &v) noexcept { cin >> v; }
inline void InputF(char &v) noexcept {
  while (isspace(v = gc()))
    ;
}
inline void InputF(string &v) noexcept {
  char c;
  for (InputF(c); !isspace(c); c = gc())
    v += c;
}
inline void InputF(int &v) noexcept {
  bool neg = false;
  v = 0;
  char c;
  InputF(c);
  if (c == '-') {
    neg = true;
    c = gc();
  }
  for (; isdigit(c); c = gc())
    v = v * 10 + (c - '0');
  if (neg)
    v = -v;
}
inline void InputF(long long &v) noexcept {
  bool neg = false;
  v = 0;
  char c;
  InputF(c);
  if (c == '-') {
    neg = true;
    c = gc();
  }
  for (; isdigit(c); c = gc())
    v = v * 10 + (LL)(c - '0');
  if (neg)
    v = -v;
}
inline void InputF(double &v) noexcept {
  double dp = 1;
  bool neg = false, adp = false;
  v = 0;
  char c;
  InputF(c);
  if (c == '-') {
    neg = true;
    c = gc();
  }
  for (; isdigit(c) || c == '.'; c = gc()) {
    if (c == '.')
      adp = true;
    else if (adp)
      v += (c - '0') * (dp *= 0.1);
    else
      v = v * 10 + (c - '0');
  }
  if (neg)
    v = -v;
}
template <class T, class U> inline void InputF(pair<T, U> &v) noexcept {
  InputF(v.first);
  InputF(v.second);
}
template <class T> inline void InputF(vector<T> &v) noexcept {
  for (auto &e : v)
    InputF(e);
}
template <class T> inline T input() {
  T v;
  InputF(v);
  return v;
}
struct InputV {
  int n, m;
  InputV(int N) : n(N), m(0) {}
  InputV(pair<int, int> N) : n(N.first), m(N.second) {}
  template <class T> operator vector<T>() noexcept {
    vector<T> v(n);
    InputF(v);
    return v;
  }
  template <class T> operator vector<vector<T>>() noexcept {
    vector<vector<T>> v(n, vector<T>(m));
    InputF(v);
    return v;
  }
};
struct Input {
  template <class T> operator T() noexcept { return input<T>(); }
  int operator--(int) {
    int v;
    InputF(v);
    v--;
    return v;
  }
  InputV operator[](int n) noexcept { return InputV(n); }
  InputV operator[](pair<int, int> n) noexcept { return InputV(n); }
  void operator()() {}
  template <class H, class... T> void operator()(H &&h, T &&...t) {
    InputF(h);
    operator()(forward<T>(t)...);
  }
  template <class T, size_t W> array<vector<T>, W> get(int H) {
    array<vector<T>, W> ret;
    rep(i, H) rep(j, W) {
      T x = *this;
      ret[j].push_back(x);
    }
    return ret;
  }
} in;

// --- output --- //
struct BoolStr {
  const char *t, *f;
  BoolStr(const char *_t, const char *_f) : t(_t), f(_f) {}
} Yes("Yes", "No"), yes("yes", "no"), YES("YES", "NO"), Int("1", "0");
struct DivStr {
  const char *d, *l;
  DivStr(const char *_d, const char *_l) : d(_d), l(_l) {}
} spc(" ", "\n"), no_spc("", "\n"), end_line("\n", "\n"), comma(",", "\n"),
    no_endl(" ", "");
class Output {
  BoolStr B{Yes};
  DivStr D{spc};
  bool isPrint = true;
  void p(double v) { printf("%.20f", v); }
  void p(long double v) { printf("%.20Lf", v); }
  void p(int v) { printf("%d", v); }
  void p(LL v) { printf("%lld", v); }
  void p(char v) { putchar(v); }
  void p(bool v) { printf("%s", v ? B.t : B.f); }
  template <class T> void p(const T &v) { cout << v; }
  template <class T, class U> void p(const pair<T, U> &v) {
    p(v.first);
    printf("%s", D.d);
    p(v.second);
  }
  template <class T> void p(const vector<T> &v) {
    rep(i, sz(v)) {
      if (i)
        printf("%s", D.d);
      p(v[i]);
    }
  }
  template <class T> void p(const vector<vector<T>> &v) {
    rep(i, sz(v)) {
      if (i)
        printf("%s", D.l);
      p(v[i]);
    }
  }
  void p(const BoolStr &v) {
    B = v;
    isPrint = false;
  }
  void p(const DivStr &v) {
    D = v;
    isPrint = false;
  }

public:
  void operator()() { printf("%s", D.l); }
  template <class H> void operator()(H &&h) {
    p(h);
    if (isPrint)
      printf("%s", D.l);
    isPrint = true;
    B = Yes;
    D = spc;
  }
  template <class H, class... T> void operator()(H &&h, T &&...t) {
    p(h);
    if (isPrint)
      printf("%s", D.d);
    isPrint = true;
    operator()(forward<T>(t)...);
  }
  template <class... T> void exit(T &&...t) {
    operator()(forward<T>(t)...);
    std::exit(EXIT_SUCCESS);
  }
} out;

// --- dump --- //
#if __has_include("dump.hpp")
#include "dump.hpp"
#else
#define dump(...) (void(0))
#endif

// ---------------------------------------------------------------- //

template <class T> class Combi {
  int Max;
  T Mod;
  vector<T> fac, finv, inv;

public:
  Combi(int _Max = 300009, T _Mod = 1000000007)
      : Max(_Max), Mod(_Mod), fac(_Max), finv(_Max), inv(_Max) {
    fac[0] = fac[1] = 1;
    finv[0] = finv[1] = 1;
    inv[1] = 1;
    for (int 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;
    }
  }
  T P(int n, int r) { // n P r
    return (n < r || n < 0 || r < 0) ? 0 : fac[n] * finv[n - r] % MOD;
  }
  T C(int n, int r) { // n C r
    return (n < r || n < 0 || r < 0)
               ? 0
               : fac[n] * (finv[r] * finv[n - r] % Mod) % Mod;
  }
  T H(int n, int r) { // n H r
    return (n < 0 || r < 0) ? 0 : r == 0 ? 1 : C(n + r - 1, r);
  }
};

template <int MOD> struct modint {
  using T = long long;
  T n;
  constexpr modint(const T x = 0) noexcept : n(x % MOD) {
    if (n < 0)
      n += MOD;
  }
  constexpr int get_mod() const noexcept { return MOD; }
  constexpr modint operator+() const noexcept { return *this; }
  constexpr modint operator-() const noexcept { return n ? MOD - n : 0; }
  constexpr modint &operator++() noexcept {
    if (MOD <= ++n)
      n = 0;
    return *this;
  }
  constexpr modint &operator--() noexcept {
    if (n <= 0)
      n = MOD;
    n--;
    return *this;
  }
  constexpr modint &operator++(int) noexcept {
    modint t = *this;
    ++*this;
    return t;
  }
  constexpr modint &operator--(int) noexcept {
    modint t = *this;
    --*this;
    return t;
  }
  constexpr modint next() const noexcept {
    modint t = *this;
    ++t;
    return t;
  }
  constexpr modint pred() const noexcept {
    modint t = *this;
    --t;
    return t;
  }
  constexpr modint operator+(const modint &m) const noexcept {
    return modint(*this) += m;
  }
  constexpr modint operator-(const modint &m) const noexcept {
    return modint(*this) -= m;
  }
  constexpr modint operator*(const modint &m) const noexcept {
    return modint(*this) *= m;
  }
  constexpr modint operator/(const modint &m) const noexcept {
    return modint(*this) /= m;
  }
  constexpr modint &operator+=(const modint &m) noexcept {
    n += m.n;
    if (n >= MOD)
      n -= MOD;
    return *this;
  }
  constexpr modint &operator-=(const modint &m) noexcept {
    n -= m.n;
    if (n < 0)
      n += MOD;
    return *this;
  }
  constexpr modint &operator*=(const modint &m) noexcept {
    n = n * m.n % MOD;
    return *this;
  }
  constexpr modint &operator/=(const modint &m) noexcept {
    assert(m.n);
    T a = m.n, b = MOD, u = 1, v = 0;
    while (b) {
      T t = a / b;
      a -= t * b;
      swap(a, b);
      u -= t * v;
      swap(u, v);
    }
    n = n * u % MOD;
    if (n < 0)
      n += MOD;
    return *this;
  }
  constexpr bool operator==(const modint &m) const noexcept { return n == m.n; }
  constexpr bool operator!=(const modint &m) const noexcept { return n != m.n; }
  constexpr modint pow(modint m) const noexcept {
    modint t = n, res = 1;
    while (m.n > 0) {
      if (m.n & 1)
        res *= t;
      t *= t;
      m.n >>= 1;
    }
    return res;
  }
  constexpr modint operator^(modint m) const noexcept { return pow(m); }
};
using mint = modint<1000000007>;
using VM = vector<mint>;
ostream &operator<<(ostream &os, const modint<MOD> &m) noexcept {
  return os << m.n;
}
istream &operator>>(istream &is, modint<MOD> &m) noexcept { return is >> m.n; }
mint operator""_m(unsigned long long n) { return n; }

int main() {
  int x1, y1, x2, y2;
  in(x1, y1, x2, y2);
  Combi<LL> comb(3000007);

  mint ans;
  FOR(x, x1, x2 + 1) { ans += comb.C(x + y2 + 1, y2) - comb.C(x + y1, y1 - 1); }
  out(ans);
}