Problem p02758 (C++)

Buggy submission (status: Runtime Error):
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int MOD = 998244353;
typedef pair<int, int> P;
int dp[123456]; // dp[i]:v[i+1]に干渉しv[i+2]以降に干渉しない通り
int a[123456];  // a[i]:v[i]を起動した時に干渉する右端
template <typename T> void chmax(T &a, T b) {
  if (a < b)
    a = b;
}
template <typename T> struct SegmentTree {
  using F = function<T(T, T)>;
  int n;         // セグ木の幅
  F f;           // モノイドの作用
  T ti;          // 初期値
  vector<T> dat; // セグ木で使う配列 一番上に使われるのがdat[1]
                 // 一番下の段の左端がdat[n]
  SegmentTree(){};
  SegmentTree(F f, T ti) : f(f), ti(ti) {} // fがモノイドとしての作用、tiが単位元

  SegmentTree(F f, T ti, int n_) : f(f), ti(ti) {
    n = 1;
    while (n < n_)
      n <<= 1;
    dat.assign(n << 1, ti);
  }

  SegmentTree(F f, T ti, const vector<T> &v) : f(f), ti(ti) {
    int n_ = v.size();
    n = 1;
    while (n < n_)
      n <<= 1;
    dat.assign(n << 1, ti);
    for (int i = 0; i < n_; i++)
      dat[n + i] = v[i];
    for (int i = n - 1; i; i--)
      dat[i] = f(dat[(i << 1) | 0], dat[(i << 1) | 1]);
  }

  T get(int x) { return dat[x + n]; }

  void init(int n_) { // ただの初期化
    n = 1;
    while (n < n_)
      n <<= 1;              // nが実際にとる２べきの幅
    dat.assign(n << 1, ti); // dat配列は幅の２倍あってその全てにtiを代入
  }

  void build(const vector<T> &v) { // vector vに合わせた初期化
    int n_ = v.size();
    init(n_);
    for (int i = 0; i < n_; i++)
      dat[n + i] = v[i]; // v[i]を一番子の所に代入
    for (int i = n - 1; i; i--)
      dat[i] = f(dat[(i << 1) | 0], dat[(i << 1) | 1]); // iの子供は2*i、2*i+1
  }

  void set_val(int k, T x) { // （0-indexedで）kにxを代入
    dat[k += n] = x;
    while (k >>= 1) // kがどんどん自分の親になっていく、0になったらおしまい
      dat[k] = f(dat[(k << 1) | 0], dat[(k << 1) | 1]);
  }

  T query(int a, int b) { //[a,b)の範囲のモノイドの作用を返す
    T vl = ti, vr = ti;
    for (int l = a + n, r = b + n; l < r; l >>= 1, r >>= 1) {
      if (l & 1)
        vl = f(vl, dat[l++]); // 左端を使いたくなるのはブロックが奇数の時
      if (r & 1)
        vr = f(dat[--r], vr); // 右端〃偶数の時（開区間なのでr:奇数で反応）
    }
    return f(vl, vr);
  }

  template <typename C>
  int find(int st, C &check, T &acc, int k, int l,
           int r) { // dat[k]:[l,r)の間にいるenを返す
    // 下のfind用なので先にそっち読むべき
    // accはst~lとかで確定している値
    if (l + 1 == r) { // 答えがちゃんと存在するならこのlが答え
      acc = f(acc, dat[k]);
      return check(acc) ? k - n : -1; // 入るならk-n(==l)を、入らないなら-1を返す
    }
    int m = (l + r) >> 1; // にぶたんのそれ
    if (m <= st)
      return find(st, check, acc, (k << 1) | 1, m,
                  r); // mがstより後ろなら明らかにl=m;
    if (st <= l &&
        !check(
            f(acc,
              dat[k]))) { // 今見てるk全部でもダメ(mだと足りないのでl=mにする）
      acc = f(acc, dat[k]); // kは全部accに加えた上で-1を返す
      return -1;
    }
    int vl = find(st, check, acc, (k << 1) | 0, l,
                  m); // まず[l,m)の中にenがいるか探す
    if (~vl)
      return vl; // いたならそれを返す(if(~vl)は-1以外にのみ反応する)
    return find(st, check, acc, (k << 1) | 1, m,
                r); // いなかったら[m,r)の中にenがいるのでそれを返す
    // acc=f(acc,dat[k]);がここで効いてる
  }

  template <typename C> int find(int st, C &check) {
    // check(f(dat[st]~dat[en-1]))=false,check(f(dat[st]~dat[en]))＝trueとなるenを返す
    T acc = ti;
    return find(st, check, acc, 1, 0, n);
  }
};

signed main() {
  dp[0] = 1;
  int n;
  cin >> n;
  vector<P> v(n);
  for (int i = 0; i < n; i++) {
    int x, d;
    cin >> x >> d;
    v[i] = P(x, d);
  }
  sort(v.begin(), v.end());
  int ans = 0;
  /*
  for(int i=0;i<n;i++){
    int now=i,right=v[i].first+v[i].second;
    while(now<n&&right>v[now+1].first){
      now++;
      chmax(right,v[now].first+v[now].second);
    }
    a[i]=now;
  }
  */
  auto f = [](int a, int b) { return max(a, b); };
  SegmentTree<int> seg(f, -1);
  seg.init(n + 5);
  for (int i = 0; i < n; i++)
    seg.set_val(i, i);
  for (int i = n - 2; i >= 0; i--) {
    int l = i, r = n;
    while (r - l > 1) {
      int mid = (r + l) >> 1;
      if (v[mid].first >= v[i].first + v[i].second)
        r = mid;
      else
        l = mid;
    }
    seg.set_val(i, seg.query(i, r));
  }
  for (int i = 0; i < n; i++) {
    dp[seg.get(i) + 1] += dp[i];
    dp[seg.get(i) + 1] %= MOD;
    dp[i + 1] += dp[i];
    dp[i + 1] %= MOD;
  }
  cout << (dp[n] + dp[n + 1]) % MOD << endl;
}


Fixed submission:
#include <bits/stdc++.h>
using namespace std;
#define int long long
const int MOD = 998244353;
typedef pair<int, int> P;
int dp[223456]; // dp[i]:v[i+1]に干渉しv[i+2]以降に干渉しない通り
int a[223456];  // a[i]:v[i]を起動した時に干渉する右端
template <typename T> void chmax(T &a, T b) {
  if (a < b)
    a = b;
}
template <typename T> struct SegmentTree {
  using F = function<T(T, T)>;
  int n;         // セグ木の幅
  F f;           // モノイドの作用
  T ti;          // 初期値
  vector<T> dat; // セグ木で使う配列 一番上に使われるのがdat[1]
                 // 一番下の段の左端がdat[n]
  SegmentTree(){};
  SegmentTree(F f, T ti) : f(f), ti(ti) {} // fがモノイドとしての作用、tiが単位元

  SegmentTree(F f, T ti, int n_) : f(f), ti(ti) {
    n = 1;
    while (n < n_)
      n <<= 1;
    dat.assign(n << 1, ti);
  }

  SegmentTree(F f, T ti, const vector<T> &v) : f(f), ti(ti) {
    int n_ = v.size();
    n = 1;
    while (n < n_)
      n <<= 1;
    dat.assign(n << 1, ti);
    for (int i = 0; i < n_; i++)
      dat[n + i] = v[i];
    for (int i = n - 1; i; i--)
      dat[i] = f(dat[(i << 1) | 0], dat[(i << 1) | 1]);
  }

  T get(int x) { return dat[x + n]; }

  void init(int n_) { // ただの初期化
    n = 1;
    while (n < n_)
      n <<= 1;              // nが実際にとる２べきの幅
    dat.assign(n << 1, ti); // dat配列は幅の２倍あってその全てにtiを代入
  }

  void build(const vector<T> &v) { // vector vに合わせた初期化
    int n_ = v.size();
    init(n_);
    for (int i = 0; i < n_; i++)
      dat[n + i] = v[i]; // v[i]を一番子の所に代入
    for (int i = n - 1; i; i--)
      dat[i] = f(dat[(i << 1) | 0], dat[(i << 1) | 1]); // iの子供は2*i、2*i+1
  }

  void set_val(int k, T x) { // （0-indexedで）kにxを代入
    dat[k += n] = x;
    while (k >>= 1) // kがどんどん自分の親になっていく、0になったらおしまい
      dat[k] = f(dat[(k << 1) | 0], dat[(k << 1) | 1]);
  }

  T query(int a, int b) { //[a,b)の範囲のモノイドの作用を返す
    T vl = ti, vr = ti;
    for (int l = a + n, r = b + n; l < r; l >>= 1, r >>= 1) {
      if (l & 1)
        vl = f(vl, dat[l++]); // 左端を使いたくなるのはブロックが奇数の時
      if (r & 1)
        vr = f(dat[--r], vr); // 右端〃偶数の時（開区間なのでr:奇数で反応）
    }
    return f(vl, vr);
  }

  template <typename C>
  int find(int st, C &check, T &acc, int k, int l,
           int r) { // dat[k]:[l,r)の間にいるenを返す
    // 下のfind用なので先にそっち読むべき
    // accはst~lとかで確定している値
    if (l + 1 == r) { // 答えがちゃんと存在するならこのlが答え
      acc = f(acc, dat[k]);
      return check(acc) ? k - n : -1; // 入るならk-n(==l)を、入らないなら-1を返す
    }
    int m = (l + r) >> 1; // にぶたんのそれ
    if (m <= st)
      return find(st, check, acc, (k << 1) | 1, m,
                  r); // mがstより後ろなら明らかにl=m;
    if (st <= l &&
        !check(
            f(acc,
              dat[k]))) { // 今見てるk全部でもダメ(mだと足りないのでl=mにする）
      acc = f(acc, dat[k]); // kは全部accに加えた上で-1を返す
      return -1;
    }
    int vl = find(st, check, acc, (k << 1) | 0, l,
                  m); // まず[l,m)の中にenがいるか探す
    if (~vl)
      return vl; // いたならそれを返す(if(~vl)は-1以外にのみ反応する)
    return find(st, check, acc, (k << 1) | 1, m,
                r); // いなかったら[m,r)の中にenがいるのでそれを返す
    // acc=f(acc,dat[k]);がここで効いてる
  }

  template <typename C> int find(int st, C &check) {
    // check(f(dat[st]~dat[en-1]))=false,check(f(dat[st]~dat[en]))＝trueとなるenを返す
    T acc = ti;
    return find(st, check, acc, 1, 0, n);
  }
};

signed main() {
  dp[0] = 1;
  int n;
  cin >> n;
  vector<P> v(n);
  for (int i = 0; i < n; i++) {
    int x, d;
    cin >> x >> d;
    v[i] = P(x, d);
  }
  sort(v.begin(), v.end());
  int ans = 0;
  /*
  for(int i=0;i<n;i++){
    int now=i,right=v[i].first+v[i].second;
    while(now<n&&right>v[now+1].first){
      now++;
      chmax(right,v[now].first+v[now].second);
    }
    a[i]=now;
  }
  */
  auto f = [](int a, int b) { return max(a, b); };
  SegmentTree<int> seg(f, -1);
  seg.init(n + 5);
  for (int i = 0; i < n; i++)
    seg.set_val(i, i);
  for (int i = n - 2; i >= 0; i--) {
    int l = i, r = n;
    while (r - l > 1) {
      int mid = (r + l) >> 1;
      if (v[mid].first >= v[i].first + v[i].second)
        r = mid;
      else
        l = mid;
    }
    seg.set_val(i, seg.query(i, r));
  }
  for (int i = 0; i < n; i++) {
    dp[seg.get(i) + 1] += dp[i];
    dp[seg.get(i) + 1] %= MOD;
    dp[i + 1] += dp[i];
    dp[i + 1] %= MOD;
  }
  cout << (dp[n] + dp[n + 1]) % MOD << endl;
}