Problem p03176 (C++)

Buggy submission (status: Runtime Error):
#include <algorithm>
#include <cstdlib>
#include <iostream>
#include <map>
#include <math.h>
#include <numeric>
#include <queue>
#include <set>
#include <unordered_map>

using namespace std;

#define fs first
#define sc second

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

const ll mod = 1e9 + 7;
ll fact[200200];
ll invfact[200200];

inline ll take_mod(ll a) { return (a % mod + mod) % mod; }

inline ll add(ll a, ll b) { return take_mod(a + b); }

inline ll sub(ll a, ll b) { return take_mod(a - b); }

inline ll mul(ll a, ll b) { return take_mod(a * b); }

inline ll pow(ll x, ll n) {
  ll res = 1LL;
  while (n > 0) {
    if (n & 1)
      res = mul(res, x);
    x = mul(x, x);
    n >>= 1;
  }
  return res;
}

ll mod_inv(ll x) { return pow(x, mod - 2); }

// nは上限
void make_fact(ll n) {
  fact[0] = 1;
  ll res = 1;
  for (int i = 1; i <= n; i++) {
    fact[i] = res;
    res = mul(res, i + 1);
  }
}

// nは上限
void make_invfact(ll n) {
  invfact[0] = 1;
  invfact[n] = mod_inv(fact[n]);
  for (int i = n - 1; i >= 1; i--) {
    invfact[i] = mul(invfact[i + 1], i + 1);
  }
}

ll perm(ll n, ll k) { return mul(fact[n], invfact[n - k]); }

ll comb(ll n, ll k) { return mul(mul(fact[n], invfact[n - k]), invfact[k]); }

#define MAX_N 210000
template <typename T> class SegTree {
public:
  // セグメント木を持つグローバル変数
  int n;
  T dat[2 * MAX_N - 1];

  void init(int n_, T init_num) {
    n = 1;
    while (n < n_)
      n *= 2;
    for (int i = 0; i < 2 * n - 1; i++) {
      dat[i] = init_num;
    }
  }

  // k番目の値(0-indexed)をaに変更
  void update(int k, T a) {
    k += n - 1;
    dat[k] = a;
    while (k > 0) {
      k = (k - 1) / 2;
      dat[k] = max(dat[k * 2 + 1], dat[k * 2 + 2]);
    }
  }

  // [a, b)の最小値を求める。
  // kは節点の番号、l, rはその節点が[l, r)に対応づいていることを表す。
  // query(a, b, 0, 0, n)で外からは呼び出す。
private:
  T query(int a, int b, int k, int l, int r) {
    if (r <= a || b <= l) {
      return 0;
    }

    if (a <= l && r <= b) {
      return dat[k];
    } else {
      ll vl = query(a, b, k * 2 + 1, l, (l + r) / 2);
      ll vr = query(a, b, k * 2 + 2, (l + r) / 2, r);
      return max(vl, vr);
    }
  }
  // [a, b)の最小値を求める。
public:
  T query(int a, int b) { return query(a, b, 0, 0, this->n); }
};

int main() {
  int N, h[210000];
  ll a[210000];
  cin >> N;
  h[0] = 0;
  a[0] = 0;
  for (int i = 1; i <= N; i++) {
    cin >> h[i];
  }
  for (int i = 1; i <= N; i++) {
    cin >> a[i];
  }

  SegTree<ll> seg;
  seg.init(N + 1, 0);

  ll res = 0;
  for (int i = 1; i <= N; i++) {
    ll tmp = seg.query(0, h[i]) + a[i];
    seg.update(h[i], tmp);
    res = max(res, tmp);
  }

  cout << res << endl;

  return 0;
}

Fixed submission:
#include <algorithm>
#include <cstdlib>
#include <iostream>
#include <map>
#include <math.h>
#include <numeric>
#include <queue>
#include <set>
#include <unordered_map>

using namespace std;

#define fs first
#define sc second

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

const ll mod = 1e9 + 7;
ll fact[200200];
ll invfact[200200];

inline ll take_mod(ll a) { return (a % mod + mod) % mod; }

inline ll add(ll a, ll b) { return take_mod(a + b); }

inline ll sub(ll a, ll b) { return take_mod(a - b); }

inline ll mul(ll a, ll b) { return take_mod(a * b); }

inline ll pow(ll x, ll n) {
  ll res = 1LL;
  while (n > 0) {
    if (n & 1)
      res = mul(res, x);
    x = mul(x, x);
    n >>= 1;
  }
  return res;
}

ll mod_inv(ll x) { return pow(x, mod - 2); }

// nは上限
void make_fact(ll n) {
  fact[0] = 1;
  ll res = 1;
  for (int i = 1; i <= n; i++) {
    fact[i] = res;
    res = mul(res, i + 1);
  }
}

// nは上限
void make_invfact(ll n) {
  invfact[0] = 1;
  invfact[n] = mod_inv(fact[n]);
  for (int i = n - 1; i >= 1; i--) {
    invfact[i] = mul(invfact[i + 1], i + 1);
  }
}

ll perm(ll n, ll k) { return mul(fact[n], invfact[n - k]); }

ll comb(ll n, ll k) { return mul(mul(fact[n], invfact[n - k]), invfact[k]); }

#define MAX_N 410000
template <typename T> class SegTree {
public:
  // セグメント木を持つグローバル変数
  int n;
  T dat[2 * MAX_N - 1];

  void init(int n_, T init_num) {
    n = 1;
    while (n < n_)
      n *= 2;
    for (int i = 0; i < 2 * n - 1; i++) {
      dat[i] = init_num;
    }
  }

  // k番目の値(0-indexed)をaに変更
  void update(int k, T a) {
    k += n - 1;
    dat[k] = a;
    while (k > 0) {
      k = (k - 1) / 2;
      dat[k] = max(dat[k * 2 + 1], dat[k * 2 + 2]);
    }
  }

  // [a, b)の最小値を求める。
  // kは節点の番号、l, rはその節点が[l, r)に対応づいていることを表す。
  // query(a, b, 0, 0, n)で外からは呼び出す。
private:
  T query(int a, int b, int k, int l, int r) {
    if (r <= a || b <= l) {
      return 0;
    }

    if (a <= l && r <= b) {
      return dat[k];
    } else {
      ll vl = query(a, b, k * 2 + 1, l, (l + r) / 2);
      ll vr = query(a, b, k * 2 + 2, (l + r) / 2, r);
      return max(vl, vr);
    }
  }
  // [a, b)の最小値を求める。
public:
  T query(int a, int b) { return query(a, b, 0, 0, this->n); }
};

int main() {
  int N, h[210000];
  ll a[210000];
  cin >> N;
  h[0] = 0;
  a[0] = 0;
  for (int i = 1; i <= N; i++) {
    cin >> h[i];
  }
  for (int i = 1; i <= N; i++) {
    cin >> a[i];
  }

  SegTree<ll> seg;
  seg.init(N + 1, 0);

  ll res = 0;
  for (int i = 1; i <= N; i++) {
    ll tmp = seg.query(0, h[i]) + a[i];
    seg.update(h[i], tmp);
    res = max(res, tmp);
  }

  cout << res << endl;

  return 0;
}