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

/*    特徴    */

/*    教訓    */

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

using ll = long long;
using ld = long double;
using pint = pair<int, int>;
using pll = pair<ll, ll>;
using pld = pair<ld, ld>;
const int INF = 1e9 + 7;
const ll LINF = 9223372036854775807;
const ll MOD = 1e9 + 7;
const ld PI = acos(-1);
const ld EPS = 1e-9; // 微調整用（EPSより小さいと0と判定など）

int ii() {
  int x;
  if (scanf("%d", &x) == 1)
    return x;
  else
    return 0;
}
long long il() {
  long long x;
  if (scanf("%lld", &x) == 1)
    return x;
  else
    return 0;
}
string is() {
  string x;
  cin >> x;
  return x;
}
char ic() {
  char x;
  cin >> x;
  return x;
}
void oi(int x) { printf("%d ", x); }
void ol(long long x) { printf("%lld ", x); }
void od_nosp(double x) { printf("%.15f", x); } // 古い問題用
void od(double x) { printf("%.15f ", x); }
// long doubleで受け取り、fをLfなどに変えて出力すると、変な数値が出る
// それをなんとかするには独自の出力を作らなければならなそう
void os(const string &s) { printf("%s ", s.c_str()); }
void oc(const char &c) { printf("%c ", c); }
#define o_map(v)                                                               \\
  {                                                                            \\
    cerr << #v << endl;                                                        \\
    for (const auto &xxx : v) {                                                \\
      cout << xxx.first << " " << xxx.second << "\\n";                          \\
    }                                                                          \\
  } // 動作未確認
void br() { putchar('\\n'); }

#define gcd __gcd // llは受け取らない C++17~のgcdと違うので注意
int lcm(int a, int b) { return a / gcd(a, b) * b; }
#define begin_end(a) a.begin(), a.end() // sort(begin_end(vec));
#define REP(i, m, n) for (ll i = (ll)(m); i < (ll)(n); i++)
#define rep(i, n) REP(i, 0, n)
#define m_p(a, b) make_pair(a, b)
#define SORT_UNIQUE(c)                                                         \\
  (sort(c.begin(), c.end()),                                                   \\
   c.resize(distance(c.begin(), unique(c.begin(), c.end()))))
#define GET_POS(c, x) (lower_bound(c.begin(), c.end(), x) - c.begin())
#define p_b push_back
#define SZ(x) ((int)(x).size) // size()がunsignedなのでエラー避けに

// 4近傍（上下左右） rep(i, 2) にすると右・下だけに進む
vector<int> dx_4 = {1, 0, -1, 0};
vector<int> dy_4 = {0, 1, 0, -1};

// coutによるpairの出力（空白区切り）
template <typename T1, typename T2>
ostream &operator<<(ostream &s, const pair<T1, T2> &p) {
  return s << "(" << p.first << " " << p.second << ")";
}
// coutによるvectorの出力（空白区切り）
template <typename T> ostream &operator<<(ostream &s, const vector<T> &v) {
  int len = v.size();
  for (int i = 0; i < len; ++i) {
    s << v[i];
    if (i < len - 1)
      s << " "; //"\\t"に変えるとTabで見やすく区切る
  }
  return s;
}
// coutによる多次元vectorの出力（空白区切り）
template <typename T>
ostream &operator<<(ostream &s, const vector<vector<T>> &vv) {
  int len = vv.size();
  for (int i = 0; i < len; ++i) {
    s << vv[i] << endl;
  }
  return s;
}

// 最大値、最小値の更新。更新したor等しければtrueを返す
template <typename T> bool chmax(T &a, T b) { return (a = max(a, b)) == b; }
template <typename T> bool chmin(T &a, T b) { return (a = min(a, b)) == b; }

// -------- template end - //

// - library ------------- //

// prefix : ceil でスニペットにしてある
// ゼロ除算エラーの対策をしていないので注意
// 関数の引数の評価順序が未規定のため、ceil_divide(ii(), ii())はNG
ll ceil_divide(ll a, ll b) { return (a + b - 1) / b; }

// --------- library end - //

// - main() -------------- //

int main() {

  ll N = il();
  vector<ll> As;
  rep(i, N + 1) As.p_b(il());

  if (N == 0) {
    if (As[0] == 1)
      cout << 1 << endl;
    else
      cout << -1 << endl;
    return 0;
  }

  // ここから、深さ1以上の場合
  if (As[0] != 0) {
    cout << -1 << endl;
    return 0;
  }

  // 最大の頂点数を下から算出
  ll nowmax = 0;
  vector<ll> max_v(N + 1, 0);
  for (ll i = N; i >= 0; --i) {
    nowmax += As[i];
    max_v[i] = nowmax;
  }

  // 算出
  ll nownode = 1;
  vector<ll> ans_v(N + 1, 1);
  REP(i, 1, N + 1) {
    nownode = (nownode - As[i - 1]) * 2;
    if (i < N && nownode <= 0) {
      cout << -1 << endl;
      return 0;
    }
    if (i == N && (nownode - As[N]) < 0) {
      cout << -1 << endl;
      return 0;
    }
    chmin(nownode, max_v[i]);
    ans_v[i] = nownode;
  }
  // 矛盾が無いか見る
  for (ll i = N; i >= 1; --i) {
    assert(ceil_divide(ans_v[i], 2) + As[i - 1] <= ans_v[i - 1]);
    assert(ans_v[i] + As[i - 1] >= ans_v[i - 1]);
  }

  // 合計
  ll ans = 0;
  for (ll eachv : ans_v)
    ans += eachv;

  cout << ans << endl;

  //

  // ---------- main() end - //
}


Error: 0

Fixed submission:

/*    特徴    */

/*    教訓    */

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

using ll = long long;
using ld = long double;
using pint = pair<int, int>;
using pll = pair<ll, ll>;
using pld = pair<ld, ld>;
const int INF = 1e9 + 7;
const ll LINF = 9223372036854775807;
const ll MOD = 1e9 + 7;
const ld PI = acos(-1);
const ld EPS = 1e-9; // 微調整用（EPSより小さいと0と判定など）

int ii() {
  int x;
  if (scanf("%d", &x) == 1)
    return x;
  else
    return 0;
}
long long il() {
  long long x;
  if (scanf("%lld", &x) == 1)
    return x;
  else
    return 0;
}
string is() {
  string x;
  cin >> x;
  return x;
}
char ic() {
  char x;
  cin >> x;
  return x;
}
void oi(int x) { printf("%d ", x); }
void ol(long long x) { printf("%lld ", x); }
void od_nosp(double x) { printf("%.15f", x); } // 古い問題用
void od(double x) { printf("%.15f ", x); }
// long doubleで受け取り、fをLfなどに変えて出力すると、変な数値が出る
// それをなんとかするには独自の出力を作らなければならなそう
void os(const string &s) { printf("%s ", s.c_str()); }
void oc(const char &c) { printf("%c ", c); }
#define o_map(v)                                                               \\
  {                                                                            \\
    cerr << #v << endl;                                                        \\
    for (const auto &xxx : v) {                                                \\
      cout << xxx.first << " " << xxx.second << "\\n";                          \\
    }                                                                          \\
  } // 動作未確認
void br() { putchar('\\n'); }

#define gcd __gcd // llは受け取らない C++17~のgcdと違うので注意
int lcm(int a, int b) { return a / gcd(a, b) * b; }
#define begin_end(a) a.begin(), a.end() // sort(begin_end(vec));
#define REP(i, m, n) for (ll i = (ll)(m); i < (ll)(n); i++)
#define rep(i, n) REP(i, 0, n)
#define m_p(a, b) make_pair(a, b)
#define SORT_UNIQUE(c)                                                         \\
  (sort(c.begin(), c.end()),                                                   \\
   c.resize(distance(c.begin(), unique(c.begin(), c.end()))))
#define GET_POS(c, x) (lower_bound(c.begin(), c.end(), x) - c.begin())
#define p_b push_back
#define SZ(x) ((int)(x).size) // size()がunsignedなのでエラー避けに

// 4近傍（上下左右） rep(i, 2) にすると右・下だけに進む
vector<int> dx_4 = {1, 0, -1, 0};
vector<int> dy_4 = {0, 1, 0, -1};

// coutによるpairの出力（空白区切り）
template <typename T1, typename T2>
ostream &operator<<(ostream &s, const pair<T1, T2> &p) {
  return s << "(" << p.first << " " << p.second << ")";
}
// coutによるvectorの出力（空白区切り）
template <typename T> ostream &operator<<(ostream &s, const vector<T> &v) {
  int len = v.size();
  for (int i = 0; i < len; ++i) {
    s << v[i];
    if (i < len - 1)
      s << " "; //"\\t"に変えるとTabで見やすく区切る
  }
  return s;
}
// coutによる多次元vectorの出力（空白区切り）
template <typename T>
ostream &operator<<(ostream &s, const vector<vector<T>> &vv) {
  int len = vv.size();
  for (int i = 0; i < len; ++i) {
    s << vv[i] << endl;
  }
  return s;
}

// 最大値、最小値の更新。更新したor等しければtrueを返す
template <typename T> bool chmax(T &a, T b) { return (a = max(a, b)) == b; }
template <typename T> bool chmin(T &a, T b) { return (a = min(a, b)) == b; }

// -------- template end - //

// - library ------------- //

// prefix : ceil でスニペットにしてある
// ゼロ除算エラーの対策をしていないので注意
// 関数の引数の評価順序が未規定のため、ceil_divide(ii(), ii())はNG
ll ceil_divide(ll a, ll b) { return (a + b - 1) / b; }

// --------- library end - //

// - main() -------------- //

int main() {

  ll N = il();
  vector<ll> As;
  rep(i, N + 1) As.p_b(il());

  if (N == 0) {
    if (As[0] == 1)
      cout << 1 << endl;
    else
      cout << -1 << endl;
    return 0;
  }

  // ここから、深さ1以上の場合
  if (As[0] != 0) {
    cout << -1 << endl;
    return 0;
  }

  // 最大の頂点数を下から算出
  ll nowmax = 0;
  vector<ll> max_v(N + 1, 0);
  for (ll i = N; i >= 0; --i) {
    nowmax += As[i];
    max_v[i] = nowmax;
  }

  // 算出
  ll nownode = 1;
  vector<ll> ans_v(N + 1, 1);
  REP(i, 1, N + 2) {
    nownode = (nownode - As[i - 1]) * 2;
    if (i < N && nownode <= 0) {
      cout << -1 << endl;
      return 0;
    }
    if (i == N && nownode < 0) {
      cout << -1 << endl;
      return 0;
    }
    chmin(nownode, max_v[i]);
    ans_v[i] = nownode;
  }
  // 矛盾が無いか見る
  for (ll i = N; i >= 1; --i) {
    assert(ceil_divide(ans_v[i], 2) + As[i - 1] <= ans_v[i - 1]);
    assert(ans_v[i] + As[i - 1] >= ans_v[i - 1]);
  }

  // 合計
  ll ans = 0;
  for (ll eachv : ans_v)
    ans += eachv;

  cout << ans << endl;

  //

  // ---------- main() end - //
}