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

#include <bits/stdc++.h>
// #include <boost/multiprecision/cpp_int.hpp>
using namespace std;
// using Graph = vector<vector<edge>>;
using dou = long double;
string yes = "yes";
string Yes = "Yes";
string YES = "YES";
string no = "no";
string No = "No";
string NO = "NO";

template <class T> inline bool chmax(T &a, T b) {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class T> inline bool chmin(T &a, T b) {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
typedef long long ll;
typedef pair<int, int> P;
typedef pair<ll, ll> PL;
// const ll mod = 1000000007;
const ll mod = 1000000007ll;
struct mint {
  ll x; // typedef long long ll;
  mint(ll x = 0) : x((x % mod + mod) % mod) {}
  mint operator-() const { return mint(-x); }
  mint &operator+=(const mint a) {
    if ((x += a.x) >= mod)
      x -= mod;
    return *this;
  }
  mint &operator-=(const mint a) {
    if ((x += mod - a.x) >= mod)
      x -= mod;
    return *this;
  }
  mint &operator*=(const mint a) {
    (x *= a.x) %= mod;
    return *this;
  }
  mint operator+(const mint a) const { return mint(*this) += a; }
  mint operator-(const mint a) const { return mint(*this) -= a; }
  mint operator*(const mint a) const { return mint(*this) *= a; }
  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
  mint inv() const { return pow(mod - 2); }
  mint &operator/=(const mint a) { return *this *= a.inv(); }
  mint operator/(const mint a) const { return mint(*this) /= a; }
};
istream &operator>>(istream &is, const mint &a) { return is >> a.x; }
ostream &operator<<(ostream &os, const mint &a) { return os << a.x; }

#define rep(i, n) for (ll i = 0; i < (ll)(n); i++)
#define brep(n) for (int bit = 0; bit < (1 << n); bit++)
#define bbrep(n) for (int bbit = 0; bbit < (1 << n); bbit++)
#define erep(i, container) for (auto &i : container)
#define itrep(i, container) for (auto i : container)
#define irep(i, n) for (ll i = n - 1; i >= (ll)0ll; i--)
#define rrep(i, m, n) for (ll i = m; i < (ll)(n); i++)
#define reprep(i, j, h, w) rep(i, h) rep(j, w)
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define VEC(type, name, n)                                                     \\
  std::vector<type> name(n);                                                   \\
  rep(i, n) std::cin >> name[i];
#define pb push_back
#define pf push_front
#define query                                                                  \\
  int qq;                                                                      \\
  std::cin >> qq;                                                              \\
  rep(qqq, qq)
#define lb lower_bound
#define ub upper_bound
#define fi first
#define se second
#define itn int
#define mp make_pair
#define sum(a) accumulate(all(a), 0ll)
#define keta fixed << setprecision
#define vout(a)                                                                \\
  erep(qxqxqx, a) std::cout << qxqxqx << ' ';                                  \\
  std::cout << std::endl;
#define vvector(name, typ, m, n, a)                                            \\
  vector<vector<typ>> name(m, vector<typ>(n, a))
// #define vvector(name,typ,m,n)vector<vector<typ> > name(m,vector<typ> (n))
#define vvvector(name, t, l, m, n, a)                                          \\
  vector<vector<vector<t>>> name(l, vector<vector<t>>(m, vector<t>(n, a)));
#define vvvvector(name, t, k, l, m, n, a)                                      \\
  vector<vector<vector<vector<t>>>> name(                                      \\
      k, vector<vector<vector<t>>>(l, vector<vector<t>>(m, vector<t>(n, a))));
#define case std::cout << "Case #" << qqq + 1 << ": "
#define res resize
#define as assign
#define ffor for (;;)
#define ppri(a, b) std::cout << a << " " << b << std::endl
#define pppri(a, b, c) std::cout << a << " " << b << " " << c << std::endl
#define aall(x, n) (x).begin(), (x).begin() + (n)
#define ssum(a) accumulate(a, 0ll)
#define stirng string
// #define grid_input(a,type) int h,w;std::cin >>
// h>>w;vvector(a,type,h,w,0);reprep(i,j,h,w)std::cin >> a[i][j];

// typedef long long T;
ll ceil(ll a, ll b) { return ((a + b - 1) / b); }
const int INF = 2000000000;
// const ll INF64 =3223372036854775807ll;
// const ll INF64 = 9223372036854775807ll;
const ll INF64 = 9223372036854775ll;
const ll MOD = 1000000007ll;
// const ll MOD = 1000003ll;
const ll OD = 1000000000000007ll;
const dou pi = 3.141592653589793;
long long modpow(long long a, long long n, ll MOD) { // 累乗の余剰
  long long res = 1;
  while (n > 0) {
    if (n & 1)
      res = res * a % MOD;
    a = a * a % MOD;
    n >>= 1;
  }
  return res;
}
// メモ
// tuple<hoge,huga,foo> tのget関数はget<i>(t)
// resizeはメンバ関数
// 三項演算子を覚えよう
//  a?b:cと書くと aならばbの処理をする
// i本目にbitをたてる　-> (1<<i)

int main() {
  /*
  rep(i,10){
      std::cout << modpow(2,i,8) << std::endl;
  }
  */
  int n;
  std::cin >> n;
  string s;
  std::cin >> s;
  // mint d=0;
  reverse(all(s));
  int ccnt = 0;
  rep(i, n) {
    if (s[i] == '1') {
      // d+=modpow(2,i);
      ccnt++;
    }
  }
  bool toku = 0;
  if (ccnt == 1)
    toku = 1;
  ll dp = 0, dm = 0;
  rep(i, n) {
    if (s[i] == '0')
      continue;
    dp += modpow(2, i, ccnt + 1);
    if (ccnt != 1) {
      dm += modpow(2, i, ccnt - 1);
      dm %= (ccnt - 1);
    }
    dp %= (ccnt + 1);
  }
  // ppri(dp,dm);
  reverse(all(s));
  rep(i, n) {
    if (toku) {
      if (s[i] == '1') {
        std::cout << 0 << std::endl;
        continue;
      }
    }
    int cnt = ccnt;
    ll d;
    if (s[i] == '1') {
      cnt--;
      d = dm;
      d -= modpow(2, n - 1 - i, cnt);
      d += cnt * 10000;
      d %= cnt;
    } else {
      cnt++;
      d = dp;
      d += modpow(2, n - 1 - i, cnt);
      //    std::cout << d << std::endl;
      // d+=cnt*10000;
      d %= cnt;
    }
    // std::cout << d << std::endl;
    //    ppri(d,cnt);
    if (d == 0) {
      std::cout << 1 << std::endl;
      continue;
    }
    int ans = 1;
    //  std::cout << 123 << std::endl;
    while (d != 0) {
      //    ppri(ans,d);
      ans++;
      int cnt = 0;
      rep(j, 20) {
        if (d & (1 << j))
          cnt++;
      }
      d %= cnt;
    }
    std::cout << ans << std::endl;
  }
}

Error: 0

Fixed submission:

#include <bits/stdc++.h>
// #include <boost/multiprecision/cpp_int.hpp>
using namespace std;
// using Graph = vector<vector<edge>>;
using dou = long double;
string yes = "yes";
string Yes = "Yes";
string YES = "YES";
string no = "no";
string No = "No";
string NO = "NO";

template <class T> inline bool chmax(T &a, T b) {
  if (a < b) {
    a = b;
    return true;
  }
  return false;
}
template <class T> inline bool chmin(T &a, T b) {
  if (a > b) {
    a = b;
    return true;
  }
  return false;
}
typedef long long ll;
typedef pair<int, int> P;
typedef pair<ll, ll> PL;
// const ll mod = 1000000007;
const ll mod = 1000000007ll;
struct mint {
  ll x; // typedef long long ll;
  mint(ll x = 0) : x((x % mod + mod) % mod) {}
  mint operator-() const { return mint(-x); }
  mint &operator+=(const mint a) {
    if ((x += a.x) >= mod)
      x -= mod;
    return *this;
  }
  mint &operator-=(const mint a) {
    if ((x += mod - a.x) >= mod)
      x -= mod;
    return *this;
  }
  mint &operator*=(const mint a) {
    (x *= a.x) %= mod;
    return *this;
  }
  mint operator+(const mint a) const { return mint(*this) += a; }
  mint operator-(const mint a) const { return mint(*this) -= a; }
  mint operator*(const mint a) const { return mint(*this) *= a; }
  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
  mint inv() const { return pow(mod - 2); }
  mint &operator/=(const mint a) { return *this *= a.inv(); }
  mint operator/(const mint a) const { return mint(*this) /= a; }
};
istream &operator>>(istream &is, const mint &a) { return is >> a.x; }
ostream &operator<<(ostream &os, const mint &a) { return os << a.x; }

#define rep(i, n) for (ll i = 0; i < (ll)(n); i++)
#define brep(n) for (int bit = 0; bit < (1 << n); bit++)
#define bbrep(n) for (int bbit = 0; bbit < (1 << n); bbit++)
#define erep(i, container) for (auto &i : container)
#define itrep(i, container) for (auto i : container)
#define irep(i, n) for (ll i = n - 1; i >= (ll)0ll; i--)
#define rrep(i, m, n) for (ll i = m; i < (ll)(n); i++)
#define reprep(i, j, h, w) rep(i, h) rep(j, w)
#define all(x) (x).begin(), (x).end()
#define rall(x) (x).rbegin(), (x).rend()
#define VEC(type, name, n)                                                     \\
  std::vector<type> name(n);                                                   \\
  rep(i, n) std::cin >> name[i];
#define pb push_back
#define pf push_front
#define query                                                                  \\
  int qq;                                                                      \\
  std::cin >> qq;                                                              \\
  rep(qqq, qq)
#define lb lower_bound
#define ub upper_bound
#define fi first
#define se second
#define itn int
#define mp make_pair
#define sum(a) accumulate(all(a), 0ll)
#define keta fixed << setprecision
#define vout(a)                                                                \\
  erep(qxqxqx, a) std::cout << qxqxqx << ' ';                                  \\
  std::cout << std::endl;
#define vvector(name, typ, m, n, a)                                            \\
  vector<vector<typ>> name(m, vector<typ>(n, a))
// #define vvector(name,typ,m,n)vector<vector<typ> > name(m,vector<typ> (n))
#define vvvector(name, t, l, m, n, a)                                          \\
  vector<vector<vector<t>>> name(l, vector<vector<t>>(m, vector<t>(n, a)));
#define vvvvector(name, t, k, l, m, n, a)                                      \\
  vector<vector<vector<vector<t>>>> name(                                      \\
      k, vector<vector<vector<t>>>(l, vector<vector<t>>(m, vector<t>(n, a))));
#define case std::cout << "Case #" << qqq + 1 << ": "
#define res resize
#define as assign
#define ffor for (;;)
#define ppri(a, b) std::cout << a << " " << b << std::endl
#define pppri(a, b, c) std::cout << a << " " << b << " " << c << std::endl
#define aall(x, n) (x).begin(), (x).begin() + (n)
#define ssum(a) accumulate(a, 0ll)
#define stirng string
// #define grid_input(a,type) int h,w;std::cin >>
// h>>w;vvector(a,type,h,w,0);reprep(i,j,h,w)std::cin >> a[i][j];

// typedef long long T;
ll ceil(ll a, ll b) { return ((a + b - 1) / b); }
const int INF = 2000000000;
// const ll INF64 =3223372036854775807ll;
// const ll INF64 = 9223372036854775807ll;
const ll INF64 = 9223372036854775ll;
const ll MOD = 1000000007ll;
// const ll MOD = 1000003ll;
const ll OD = 1000000000000007ll;
const dou pi = 3.141592653589793;
long long modpow(long long a, long long n) { // 累乗の余剰
  long long res = 1;
  while (n > 0) {
    if (n & 1)
      res = res * a % MOD;
    a = a * a % MOD;
    n >>= 1;
  }
  return res;
}
// メモ
// tuple<hoge,huga,foo> tのget関数はget<i>(t)
// resizeはメンバ関数
// 三項演算子を覚えよう
//  a?b:cと書くと aならばbの処理をする
// i本目にbitをたてる　-> (1<<i)

int main() {
  int n;
  std::cin >> n;
  string s;
  std::cin >> s;
  mint d = 0;
  reverse(all(s));
  int ccnt = 0;
  rep(i, n) {
    if (s[i] == '1') {
      d += modpow(2, i);
      ccnt++;
    }
  }
  // std::cout << d << std::endl;
  reverse(all(s));
  rep(i, n) {
    mint dd = d;
    int cnt = ccnt;
    //   std::cout << 1 << std::endl;
    if (s[i] == '0') {
      cnt++;
      dd += modpow(2, n - i - 1);
    } else {
      cnt--;
      dd -= modpow(2, n - i - 1);
    }
    // std::cout << dd << std::endl;
    // std::cout << 2 << std::endl;
    // ppri(2,cnt);
    // std::cout << cnt << std::endl;
    if (dd.x % cnt == 0) {
      std::cout << 1 << std::endl;
      continue;
    }
    int ans = 1;
    int dum = dd.x;
    dum %= cnt;
    // std::cout << 3 << std::endl;
    while (dum % cnt != 0) {
      ans++;
      cnt = 0;
      rep(i, 20) {
        if (dum & (1 << i))
          cnt++;
      }
      dum %= cnt;
    }
    std::cout << ans << std::endl;
    // std::cout << 4 << std::endl;
  }
}