Problem p02555 (C++) — buggy submission, verdict: Runtime Error

// include
//------------------------------------------
#include <bits/stdc++.h>
using namespace std;

// typedef
//------------------------------------------
using LL = int64_t;
using VL = vector<LL>;
using VVL = vector<VL>;
using VS = vector<string>;
using PLL = pair<LL, LL>;

// container util
//------------------------------------------
#define whole(f, x, ...)                                                       \\
  ([&](decltype((x)) whole) {                                                  \\
    return (f)(begin(whole), end(whole), ##__VA_ARGS__);                       \\
  })(x)
#define rwhole(f, x, ...)                                                      \\
  ([&](decltype((x)) whole) {                                                  \\
    return (f)(rbegin(whole), rend(whole), ##__VA_ARGS__);                     \\
  })(x)
#define EACH(i, c)                                                             \\
  for (decltype((c).begin()) i = (c).begin(); i != (c).end(); ++i)
#define EXIST(s, e) ((s).find(e) != (s).end())
#define ALL(x) ::std::begin(x), ::std::end(x)
#define RALL(x) ::std::rbegin(x), ::std::rend(x)
#define SORT(c) whole(sort, c)
#define RSORT(c) rwhole(sort, c)

// constant
//--------------------------------------------
constexpr double EPS = 1e-10;
constexpr double PI = 3.14159265358979323846;
constexpr int MOD = 1000000007;

// grid
//--------------------------------------------
VL dx = {0, 1, 0, -1};
VL dy = {1, 0, -1, 0};
VL dx2 = {-1, 0, 1, -1, 1, -1, 0, 1};
VL dy2 = {-1, -1, -1, 0, 0, 1, 1, 1};

// debug
//--------------------------------------------
#define dump(x) cerr << #x << " = " << (x) << endl;
#define debug(x)                                                               \\
  cerr << #x << " = " << (x) << " (L" << __LINE__ << ")"                       \\
       << " " << __FILE__ << endl;

// IO accelerate
//--------------------------------------------
struct InitIO {
  InitIO() {
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);
    cout << fixed << setprecision(30);
  }
} init_io;

// template
//--------------------------------------------
//  declaretion
template <typename T> istream &operator>>(istream &is, vector<T> &vec);
template <typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p);
template <typename T> ostream &operator<<(ostream &os, const vector<T> &vec);
template <typename T>
ostream &operator<<(ostream &os, const vector<vector<T>> &vv);
template <typename T> vector<T> make_v(size_t a);
template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts);
template <typename T, typename V>
typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v);
template <typename T, typename V>
typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v);

// implementation
template <typename T> istream &operator>>(istream &is, vector<T> &vec) {
  for (T &x : vec)
    is >> x;
  return is;
}
template <typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p) {
  os << p.first << "," << p.second;
  return os;
}
template <typename T> ostream &operator<<(ostream &os, const vector<T> &vec) {
  for (int i = 0; i < vec.size(); i++) {
    os << vec[i] << (i + 1 == vec.size() ? "" : "\\t");
  }
  return os;
}
template <typename T>
ostream &operator<<(ostream &s, const vector<vector<T>> &vv) {
  for (int i = 0; i < vv.size(); ++i) {
    s << vv[i] << endl;
  }
  return s;
}

// 多重vector
// auto dp=make_v<int>(4,h,w) みたいに使える
template <typename T> vector<T> make_v(size_t a) { return vector<T>(a); }

template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts) {
  return vector<decltype(make_v<T>(ts...))>(a, make_v<T>(ts...));
}

// 多重vectorのためのfill
// fill_v(dp,0) みたいに使える
template <typename T, typename V>
typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v) {
  t = v;
}

template <typename T, typename V>
typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v) {
  for (auto &e : t)
    fill_v(e, v);
}

template <class T, class U> void chmax(T &a, U b) {
  if (a < b)
    a = b;
}
template <class T, class U> void chmin(T &a, U b) {
  if (b < a)
    a = b;
}
template <typename T> T gcd(T a, T b) { return b ? gcd(b, a % b) : a; }
template <typename T> T lcm(T a, T b) { return a / gcd(a, b) * b; }

// main code
//  auto mod int
//
//  https://youtu.be/L8grWxBlIZ4?t=9858
//  https://youtu.be/ERZuLAxZffQ?t=4807 : optimize
//  https://youtu.be/8uowVvQ_-Mo?t=1329 : division
//  const int MOD = 1000000007;
struct mint {
  LL x; // typedef long long LL;
  constexpr mint(LL x = 0) : x((x % MOD + MOD) % MOD) {}
  constexpr mint operator-() const { return mint(-x); }
  constexpr mint &operator+=(const mint &a) {
    if ((this->x += a.x) >= MOD)
      this->x -= MOD;
    return *this;
  }
  constexpr mint &operator-=(const mint &a) {
    if ((this->x += MOD - a.x) >= MOD)
      this->x -= MOD;
    return *this;
  }
  constexpr mint &operator*=(const mint &a) {
    (this->x *= a.x) %= MOD;
    return *this;
  }
  constexpr mint operator+(const mint &a) const {
    mint res(*this);
    return res += a;
  }
  constexpr mint operator-(const mint &a) const {
    mint res(*this);
    return res -= a;
  }
  constexpr mint operator*(const mint &a) const {
    mint res(*this);
    return res *= a;
  }
  constexpr mint pow(LL t) const {
    if (!t)
      return 1;
    mint a = pow(t >> 1);
    a *= a;
    if (t & 1)
      a *= *this;
    return a;
  }

  // for prime MOD
  constexpr mint inv() const { return pow(MOD - 2); }
  constexpr mint &operator/=(const mint &a) { return (*this) *= a.inv(); }
  constexpr mint operator/(const mint &a) const {
    mint res(*this);
    return res /= a;
  }

  constexpr bool operator==(const mint &a) const { return this->x == a.x; }
  constexpr bool operator!=(const mint &a) const { return this->x != a.x; }
  constexpr bool operator<(const mint &a) const { return this->x < a.x; }

  friend istream &operator>>(istream &is, mint &m) { return is >> m.x; }
  friend ostream &operator<<(ostream &os, const mint &m) { return os << m.x; }
};
typedef vector<mint> VM;
typedef vector<VM> VVM;

const int MAX_C = 2000;
long long Com[MAX_C][MAX_C];

void calc_com() {
  // memset(Com, 1, sizeof(Com));
  for (int i = 0; i < MAX_C; i++) {
    for (int j = 0; j < MAX_C; j++) {
      Com[i][j] = 0;
    }
  }
  Com[0][0] = 1;
  for (int i = 1; i < MAX_C; ++i) {
    Com[i][0] = 1;
    for (int j = 1; j < MAX_C; ++j) {
      Com[i][j] = (Com[i - 1][j - 1] + Com[i - 1][j]) % MOD;
    }
  }
}

int main(int argc, char *argv[]) {
  LL s;
  cin >> s;
  if (s <= 2) {
    cout << 0 << endl;
    return 0;
  }
  mint ans = 0;
  LL d = s / 3;
  LL m = s % 3;
  calc_com();
  for (int i = 0; i < d; i++) {
    LL ball = m + i * 3;
    LL box = d - i;
    if (box <= 0) {
      ans += 1;
    } else {
      ans += Com[box + ball - 1][ball];
    }
  }
  cout << ans << endl;
  return 0;
}

Error: -11

Fixed submission:

// include
//------------------------------------------
#include <bits/stdc++.h>
using namespace std;

// typedef
//------------------------------------------
using LL = int64_t;
using VL = vector<LL>;
using VVL = vector<VL>;
using VS = vector<string>;
using PLL = pair<LL, LL>;

// container util
//------------------------------------------
#define whole(f, x, ...)                                                       \\
  ([&](decltype((x)) whole) {                                                  \\
    return (f)(begin(whole), end(whole), ##__VA_ARGS__);                       \\
  })(x)
#define rwhole(f, x, ...)                                                      \\
  ([&](decltype((x)) whole) {                                                  \\
    return (f)(rbegin(whole), rend(whole), ##__VA_ARGS__);                     \\
  })(x)
#define EACH(i, c)                                                             \\
  for (decltype((c).begin()) i = (c).begin(); i != (c).end(); ++i)
#define EXIST(s, e) ((s).find(e) != (s).end())
#define ALL(x) ::std::begin(x), ::std::end(x)
#define RALL(x) ::std::rbegin(x), ::std::rend(x)
#define SORT(c) whole(sort, c)
#define RSORT(c) rwhole(sort, c)

// constant
//--------------------------------------------
constexpr double EPS = 1e-10;
constexpr double PI = 3.14159265358979323846;
constexpr int MOD = 1000000007;

// grid
//--------------------------------------------
VL dx = {0, 1, 0, -1};
VL dy = {1, 0, -1, 0};
VL dx2 = {-1, 0, 1, -1, 1, -1, 0, 1};
VL dy2 = {-1, -1, -1, 0, 0, 1, 1, 1};

// debug
//--------------------------------------------
#define dump(x) cerr << #x << " = " << (x) << endl;
#define debug(x)                                                               \\
  cerr << #x << " = " << (x) << " (L" << __LINE__ << ")"                       \\
       << " " << __FILE__ << endl;

// IO accelerate
//--------------------------------------------
struct InitIO {
  InitIO() {
    cin.tie(nullptr);
    ios_base::sync_with_stdio(false);
    cout << fixed << setprecision(30);
  }
} init_io;

// template
//--------------------------------------------
//  declaretion
template <typename T> istream &operator>>(istream &is, vector<T> &vec);
template <typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p);
template <typename T> ostream &operator<<(ostream &os, const vector<T> &vec);
template <typename T>
ostream &operator<<(ostream &os, const vector<vector<T>> &vv);
template <typename T> vector<T> make_v(size_t a);
template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts);
template <typename T, typename V>
typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v);
template <typename T, typename V>
typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v);

// implementation
template <typename T> istream &operator>>(istream &is, vector<T> &vec) {
  for (T &x : vec)
    is >> x;
  return is;
}
template <typename T1, typename T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &p) {
  os << p.first << "," << p.second;
  return os;
}
template <typename T> ostream &operator<<(ostream &os, const vector<T> &vec) {
  for (int i = 0; i < vec.size(); i++) {
    os << vec[i] << (i + 1 == vec.size() ? "" : "\\t");
  }
  return os;
}
template <typename T>
ostream &operator<<(ostream &s, const vector<vector<T>> &vv) {
  for (int i = 0; i < vv.size(); ++i) {
    s << vv[i] << endl;
  }
  return s;
}

// 多重vector
// auto dp=make_v<int>(4,h,w) みたいに使える
template <typename T> vector<T> make_v(size_t a) { return vector<T>(a); }

template <typename T, typename... Ts> auto make_v(size_t a, Ts... ts) {
  return vector<decltype(make_v<T>(ts...))>(a, make_v<T>(ts...));
}

// 多重vectorのためのfill
// fill_v(dp,0) みたいに使える
template <typename T, typename V>
typename enable_if<is_class<T>::value == 0>::type fill_v(T &t, const V &v) {
  t = v;
}

template <typename T, typename V>
typename enable_if<is_class<T>::value != 0>::type fill_v(T &t, const V &v) {
  for (auto &e : t)
    fill_v(e, v);
}

template <class T, class U> void chmax(T &a, U b) {
  if (a < b)
    a = b;
}
template <class T, class U> void chmin(T &a, U b) {
  if (b < a)
    a = b;
}
template <typename T> T gcd(T a, T b) { return b ? gcd(b, a % b) : a; }
template <typename T> T lcm(T a, T b) { return a / gcd(a, b) * b; }

// main code
//  auto mod int
//
//  https://youtu.be/L8grWxBlIZ4?t=9858
//  https://youtu.be/ERZuLAxZffQ?t=4807 : optimize
//  https://youtu.be/8uowVvQ_-Mo?t=1329 : division
//  const int MOD = 1000000007;
struct mint {
  LL x; // typedef long long LL;
  constexpr mint(LL x = 0) : x((x % MOD + MOD) % MOD) {}
  constexpr mint operator-() const { return mint(-x); }
  constexpr mint &operator+=(const mint &a) {
    if ((this->x += a.x) >= MOD)
      this->x -= MOD;
    return *this;
  }
  constexpr mint &operator-=(const mint &a) {
    if ((this->x += MOD - a.x) >= MOD)
      this->x -= MOD;
    return *this;
  }
  constexpr mint &operator*=(const mint &a) {
    (this->x *= a.x) %= MOD;
    return *this;
  }
  constexpr mint operator+(const mint &a) const {
    mint res(*this);
    return res += a;
  }
  constexpr mint operator-(const mint &a) const {
    mint res(*this);
    return res -= a;
  }
  constexpr mint operator*(const mint &a) const {
    mint res(*this);
    return res *= a;
  }
  constexpr mint pow(LL t) const {
    if (!t)
      return 1;
    mint a = pow(t >> 1);
    a *= a;
    if (t & 1)
      a *= *this;
    return a;
  }

  // for prime MOD
  constexpr mint inv() const { return pow(MOD - 2); }
  constexpr mint &operator/=(const mint &a) { return (*this) *= a.inv(); }
  constexpr mint operator/(const mint &a) const {
    mint res(*this);
    return res /= a;
  }

  constexpr bool operator==(const mint &a) const { return this->x == a.x; }
  constexpr bool operator!=(const mint &a) const { return this->x != a.x; }
  constexpr bool operator<(const mint &a) const { return this->x < a.x; }

  friend istream &operator>>(istream &is, mint &m) { return is >> m.x; }
  friend ostream &operator<<(ostream &os, const mint &m) { return os << m.x; }
};
typedef vector<mint> VM;
typedef vector<VM> VVM;

const int MAX_C = 4000;
long long Com[MAX_C][MAX_C];

void calc_com() {
  memset(Com, 0, sizeof(Com));
  Com[0][0] = 1;
  for (int i = 1; i < MAX_C; ++i) {
    Com[i][0] = 1;
    for (int j = 1; j < MAX_C; ++j) {
      Com[i][j] = (Com[i - 1][j - 1] + Com[i - 1][j]) % MOD;
    }
  }
}

int main(int argc, char *argv[]) {
  LL s;
  cin >> s;
  if (s < 2) {
    cout << 0 << endl;
    return 0;
  }
  mint ans = 0;
  LL d = s / 3;
  LL m = s % 3;
  calc_com();
  for (int i = 0; i < d; i++) {
    LL ball = m + i * 3;
    LL box = d + m - i - 1;
    ans += Com[box + ball - 1][ball];
  }
  cout << ans << endl;
  return 0;
}