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

#include <bits/stdc++.h>
#include <stdlib.h>
#include <sys/time.h>
#include <unistd.h>

using namespace std;

typedef long long ll;
typedef pair<ll, ll> P;
typedef pair<P, ll> T;
typedef pair<long double, ll> Ps;
typedef pair<ll, bool> Pb;

const ll INF = 1e18;
const ll fact_table = 3200008;
long double Pi = 3.1415926535897932384626;

priority_queue<ll> pql;
priority_queue<P> pqp;
priority_queue<P> bag;
// big priority queue
priority_queue<ll, vector<ll>, greater<ll>> pqls;
priority_queue<P, vector<P>, greater<P>> pqps;
// small priority queue
// top pop

ll dx[8] = {1, 0, -1, 0, 1, 1, -1, -1};
ll dy[8] = {0, 1, 0, -1, 1, -1, -1, 1};
char dir[] = "DRUL";
// ↓,→,↑,←

/*
#define endl "\\n"
#ifdef ENJAPMA
#undef endl
#endif
*/
#define p(x) cout << x << endl;
#define el cout << endl;
#define pe(x) cout << x << " ";
#define ps(x) cout << fixed << setprecision(25) << x << endl;
#define pu(x) cout << (x);
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define CLEAR(a) a = decltype(a)();
#define pc(x) cout << x << ",";
#define rep(i, n) for (ll i = 0; i < (n); i++)

// const ll mod = 998244353ll;
const ll mod = 1000000007ll;

ll mypow(ll a, ll b, ll mod) {
  ll x = 1;
  while (b) {
    while (!(b & 1)) {
      (a *= a) %= mod;
      b >>= 1;
    }
    (x *= a) %= mod;
    b--;
  }
  return x;
}
void YES(bool cond) {
  if (cond) {
    p("YES");
  } else {
    p("NO");
  }
  return;
}
void Yes(bool cond) {
  if (cond) {
    p("Yes");
  } else {
    p("No");
  }
  return;
}
void line() {
  p("--------------------");
  return;
}

/*
ll fact[fact_table + 5], rfact[fact_table + 5];

void c3_init() {
    fact[0] = rfact[0] = 1;
    for (ll i = 1; i <= fact_table; i++) {
        fact[i] = (fact[i - 1] * i) % mod;
    }
    rfact[fact_table] = mypow(fact[fact_table], mod - 2, mod);
    for (ll i = fact_table; i >= 1; i--) {
        rfact[i - 1] = rfact[i] * i;
        rfact[i - 1] %= mod;
    }
    return;
}
ll c3(ll n, ll r) {
    return (((fact[n] * rfact[r]) % mod ) * rfact[n - r]) % mod;
}
*/
struct Timer {
  int64_t start;
  const int64_t CYCLES_PER_SEC = 2800000000;
  Timer() { reset(); }
  void reset() { start = getCycle(); }
  inline double get() { return (double)(getCycle() - start) / CYCLES_PER_SEC; }
  inline int64_t getCycle() {
    uint32_t low, high;
    __asm__ volatile("rdtsc" : "=a"(low), "=d"(high));
    return ((int64_t)low) | ((int64_t)high << 32);
  }
};

Timer timer;

bool multicase = false;

ll n, m, k, w, num, a, b, c, d, e, h, q, ans, sum, l;

typedef vector<ll> vec;
typedef vector<vector<ll>> mat;

string s, t;
P table[40];
ll cost[50];
ll ans_table[50];

void solve() {
  cin >> n;
  for (int i = 0; i < n; i++) {
    cin >> table[i].first >> table[i].second >> cost[i];
  }
  for (int i = 0; i <= n; i++) {
    ans_table[i] = INF;
  }
  ll bit = 1;
  for (int i = 0; i < n; i++)
    bit *= 3;
  vec v;
  for (int i = 0; i < n; i++)
    v.pb(0);
  vec vx, vy;
  for (int i = 0; i < bit; i++) {
    ll cnt = 0;
    vx.clear();
    vy.clear();
    vx.pb(0);
    vy.pb(0);
    for (int j = 0; j < n; j++) {
      if (v[j] == 1) {
        vx.pb(table[j].first);
      } else if (v[j] == 2) {
        vy.pb(table[j].second);
      }
    }
    cnt = vx.size() + vy.size() - 2;
    ll allcost = 0;
    for (int j = 0; j < n; j++) {
      ll mcost = INF;
      if (v[j] != 0) {
        mcost = 0;
      } else {
        for (int k = 0; k < vx.size(); k++) {
          mcost = min(mcost, cost[j] * abs(vx[k] - table[j].first));
        }
        for (int k = 0; k < vy.size(); k++) {
          mcost = min(mcost, cost[j] * abs(vy[k] - table[j].second));
        }
      }
      allcost += mcost;
    }
    ans_table[cnt] = min(ans_table[cnt], allcost);

    // v をなんとかする処理
    ll back_cnt = 0;
    while (v.back() == 2) {
      v.pop_back();
      back_cnt++;
      if (v.size() == 0)
        break;
    }
    if (v.size() == 0)
      break;
    v.back() += 1;
    for (int j = 0; j < back_cnt; j++) {
      v.pb(0);
    }
  }

  for (int i = 0; i <= n; i++) {
    p(ans_table[i]);
  }

  return;
}

int main() {
  // init();
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  ll q, testcase = 1;
  if (multicase) {
    cin >> q;
  } else {
    q = 1;
  }
  while (q--) {
    // pu("Case ");pu("#");pu(testcase);pu(": ");
    solve();
    testcase++;
  }
  // solve();
  return 0;
}


Error: -6
munmap_chunk(): invalid pointer


Fixed submission:

#include <bits/stdc++.h>
#include <stdlib.h>
#include <sys/time.h>
#include <unistd.h>

using namespace std;

typedef long long ll;
typedef pair<ll, ll> P;
typedef pair<P, ll> T;
typedef pair<long double, ll> Ps;
typedef pair<ll, bool> Pb;

const ll INF = 1e18;
const ll fact_table = 3200008;
long double Pi = 3.1415926535897932384626;

priority_queue<ll> pql;
priority_queue<P> pqp;
priority_queue<P> bag;
// big priority queue
priority_queue<ll, vector<ll>, greater<ll>> pqls;
priority_queue<P, vector<P>, greater<P>> pqps;
// small priority queue
// top pop

ll dx[8] = {1, 0, -1, 0, 1, 1, -1, -1};
ll dy[8] = {0, 1, 0, -1, 1, -1, -1, 1};
char dir[] = "DRUL";
// ↓,→,↑,←

/*
#define endl "\\n"
#ifdef ENJAPMA
#undef endl
#endif
*/
#define p(x) cout << x << endl;
#define el cout << endl;
#define pe(x) cout << x << " ";
#define ps(x) cout << fixed << setprecision(25) << x << endl;
#define pu(x) cout << (x);
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define CLEAR(a) a = decltype(a)();
#define pc(x) cout << x << ",";
#define rep(i, n) for (ll i = 0; i < (n); i++)

// const ll mod = 998244353ll;
const ll mod = 1000000007ll;

ll mypow(ll a, ll b, ll mod) {
  ll x = 1;
  while (b) {
    while (!(b & 1)) {
      (a *= a) %= mod;
      b >>= 1;
    }
    (x *= a) %= mod;
    b--;
  }
  return x;
}
void YES(bool cond) {
  if (cond) {
    p("YES");
  } else {
    p("NO");
  }
  return;
}
void Yes(bool cond) {
  if (cond) {
    p("Yes");
  } else {
    p("No");
  }
  return;
}
void line() {
  p("--------------------");
  return;
}

/*
ll fact[fact_table + 5], rfact[fact_table + 5];

void c3_init() {
    fact[0] = rfact[0] = 1;
    for (ll i = 1; i <= fact_table; i++) {
        fact[i] = (fact[i - 1] * i) % mod;
    }
    rfact[fact_table] = mypow(fact[fact_table], mod - 2, mod);
    for (ll i = fact_table; i >= 1; i--) {
        rfact[i - 1] = rfact[i] * i;
        rfact[i - 1] %= mod;
    }
    return;
}
ll c3(ll n, ll r) {
    return (((fact[n] * rfact[r]) % mod ) * rfact[n - r]) % mod;
}
*/
struct Timer {
  int64_t start;
  const int64_t CYCLES_PER_SEC = 2800000000;
  Timer() { reset(); }
  void reset() { start = getCycle(); }
  inline double get() { return (double)(getCycle() - start) / CYCLES_PER_SEC; }
  inline int64_t getCycle() {
    uint32_t low, high;
    __asm__ volatile("rdtsc" : "=a"(low), "=d"(high));
    return ((int64_t)low) | ((int64_t)high << 32);
  }
};

Timer timer;

bool multicase = false;

ll n, m, k, w, num, a, b, c, d, e, h, q, ans, sum, l;

typedef vector<ll> vec;
typedef vector<vector<ll>> mat;

string s, t;
P table[40];
ll cost[50];
ll ans_table[50];

void solve() {
  cin >> n;
  for (int i = 0; i < n; i++) {
    cin >> table[i].first >> table[i].second >> cost[i];
  }
  for (int i = 0; i <= n; i++) {
    ans_table[i] = INF;
  }
  ll bit = 1;
  for (int i = 0; i < n; i++)
    bit *= 3;
  vec v;
  for (int i = 0; i < n; i++)
    v.pb(0);
  vec vx, vy;
  for (int i = 0; i < bit; i++) {
    ll cnt = 0;
    vx.clear();
    vy.clear();
    vx.pb(0);
    vy.pb(0);
    for (int j = 0; j < n; j++) {
      if (v[j] == 1) {
        vx.pb(table[j].first);
      } else if (v[j] == 2) {
        vy.pb(table[j].second);
      }
    }
    cnt = vx.size() + vy.size() - 2;
    ll allcost = 0;
    for (int j = 0; j < n; j++) {
      ll mcost = INF;
      if (v[j] != 0) {
        mcost = 0;
      } else {
        for (int k = 0; k < vx.size(); k++) {
          mcost = min(mcost, cost[j] * abs(vx[k] - table[j].first));
        }
        for (int k = 0; k < vy.size(); k++) {
          mcost = min(mcost, cost[j] * abs(vy[k] - table[j].second));
        }
      }
      allcost += mcost;
    }
    ans_table[cnt] = min(ans_table[cnt], allcost);

    // v をなんとかする処理
    ll back_cnt = 0;
    while (v.back() == 2) {
      v.pop_back();
      back_cnt++;
      if (v.size() == 0)
        break;
    }
    v.back() += 1;
    for (int j = 0; j < back_cnt; j++) {
      v.pb(0);
    }
    if (v.size() == 0)
      break;
  }

  for (int i = 0; i <= n; i++) {
    p(ans_table[i]);
  }

  return;
}

int main() {
  // init();
  ios::sync_with_stdio(false);
  cin.tie(nullptr);
  ll q, testcase = 1;
  if (multicase) {
    cin >> q;
  } else {
    q = 1;
  }
  while (q--) {
    // pu("Case ");pu("#");pu(testcase);pu(": ");
    solve();
    testcase++;
  }
  // solve();
  return 0;
}