Submission #882807

# Submission time Handle Problem Language Result Execution time Memory
882807 2023-12-03T18:17:54 Z MilosMilutinovic Chorus (JOI23_chorus) C++14
61 / 100
1386 ms 12368 KB
#include <bits/stdc++.h>
 
using namespace std;

struct line {
  long long k, n;
  int f;
  void init() {
    k = 0;
    n = (long long) 1e18;
    f = 0;
  }
  long long val(long long x) {
    return k * x + n;
  }
};

const int MAX = 5e6;

line li[MAX];
int ls[MAX], rs[MAX], tsz, root;

void Update(int& x, int l, int r, line ln) {
  if (!x) {
    x = ++tsz;
    li[x].init();
  }
  int mid = l + r >> 1;
  if (li[x].val(mid) > ln.val(mid) || (ln.val(mid) == li[x].val(mid) && ln.f < li[x].f)) {
    swap(li[x], ln);
  }
  if (l == r) {
    return;
  }
  if (ln.val(l) < li[x].val(l)) {
    Update(ls[x], l, mid, ln);
  } else {
    Update(rs[x], mid + 1, r, ln);
  }
}

pair<long long, int> Query(int x, int l, int r, int i) {
  if (!x) {
    return make_pair((long long) 1e18, 0);
  }
  if (l == r) {
    return make_pair(li[x].val(i), li[x].f);
  }
  int mid = l + r >> 1;
  if (i <= mid) {
    return min(Query(ls[x], l, mid, i), make_pair(li[x].val(i), li[x].f));
  } else {
    return min(Query(rs[x], mid + 1, r, i), make_pair(li[x].val(i), li[x].f));
  }
}

void clear_seg() {
  while (tsz > 0) {
    ls[tsz] = 0;
    rs[tsz] = 0;
    li[tsz].init();
    tsz -= 1;
  }
  root = 0;
}
 
int main() {
  ios::sync_with_stdio(false);
  cin.tie(0);
  int n, k;
  cin >> n >> k;
  string s;
  cin >> s;
  vector<int> a, b;
  for (int i = 0; i < 2 * n; i++) {
    if (s[i] == 'A') {
      a.push_back(i);
    } else {
      b.push_back(i);
    }
  }
  vector<int> ia(2 * n);
  for (int i = 0; i < n; i++) {
    ia[a[i]] = i + 1;
  }
  vector<int> pref(2 * n);
  for (int i = 0; i < 2 * n; i++) {
    if (i > 0) {
      pref[i] = pref[i - 1];
    }
    pref[i] += (s[i] == 'A' ? 1 : 0);
  }
  vector<int> ib(2 * n);
  for (int i = 0; i < n; i++) {
    ib[b[i]] = pref[b[i]];
  }
  for (int i = 0; i < n; i++) {
    if (i > 0) {
      pref[i] = pref[i - 1];
    } else {
      pref[i] = 0;
    }
    pref[i] += ib[b[i]];
  }
  auto Cost = [&](int L, int R) {
    long long ret = 0;
    for (int i = L; i <= R; i++) {
      ret += max(0, ia[a[R]] - ib[b[i]]);
    }
    return ret;
  }; 
  auto Solve = [&](long long mid) {
    const long long inf = (long long) 1e18;
    vector<pair<long long, int>> dp(n + 1, {inf, 0});
    dp[0] = {0, 0};
    int ptr = -1;
    multiset<pair<long long, int>> st;
    clear_seg();
    for (int i = 0; i < n; i++) {
      st.insert(dp[i]);
      while (ptr + 1 <= i && ia[a[i]] >= ib[b[ptr + 1]]) {
        ptr += 1;
        line cur;
        cur.k = -ptr;
        cur.n = dp[ptr].first + (ptr == 0 ? 0 : pref[ptr - 1]) + mid;
        cur.f = dp[ptr].second;
        Update(root, 0, n, cur);
        st.erase(st.find(dp[ptr]));
      }
      if (!st.empty()) {
        pair<long long, int> val = *st.begin();
        dp[i + 1] = min(dp[i + 1], {val.first + mid, val.second + 1});
      }
      if (ptr != -1) {
        pair<long long, int> val = Query(root, 0, n, ia[a[i]]);
        dp[i + 1] = min(dp[i + 1], {val.first + ia[a[i]] * 1LL * (ptr + 1) - pref[ptr], val.second + 1});
      }
    }
    return dp[n];
  };
  long long low = 0, high = 1e12, p = 0;
  while (low <= high) {
    long long mid = low + high >> 1;
    pair<long long, int> res = Solve(mid);
    if (res.second <= k) {
      p = mid;
      high = mid - 1;
    } else {
      low = mid + 1;
    }
  }
  pair<long long, int> res = Solve(p);
  cout << res.first - p * k << '\n';
  return 0;
}

Compilation message

chorus.cpp: In function 'void Update(int&, int, int, line)':
chorus.cpp:28:15: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
   28 |   int mid = l + r >> 1;
      |             ~~^~~
chorus.cpp: In function 'std::pair<long long int, int> Query(int, int, int, int)':
chorus.cpp:49:15: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
   49 |   int mid = l + r >> 1;
      |             ~~^~~
chorus.cpp: In function 'int main()':
chorus.cpp:143:25: warning: suggest parentheses around '+' inside '>>' [-Wparentheses]
  143 |     long long mid = low + high >> 1;
      |                     ~~~~^~~~~~
chorus.cpp:105:8: warning: variable 'Cost' set but not used [-Wunused-but-set-variable]
  105 |   auto Cost = [&](int L, int R) {
      |        ^~~~
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2396 KB Output is correct
2 Correct 1 ms 2396 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2392 KB Output is correct
5 Correct 1 ms 4444 KB Output is correct
6 Correct 1 ms 2396 KB Output is correct
7 Correct 1 ms 2396 KB Output is correct
8 Correct 1 ms 2392 KB Output is correct
9 Correct 1 ms 2396 KB Output is correct
10 Correct 1 ms 2396 KB Output is correct
11 Correct 1 ms 2396 KB Output is correct
12 Correct 1 ms 2644 KB Output is correct
13 Correct 1 ms 2396 KB Output is correct
14 Correct 0 ms 2396 KB Output is correct
15 Correct 1 ms 2392 KB Output is correct
16 Correct 1 ms 2392 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2396 KB Output is correct
2 Correct 1 ms 2396 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2392 KB Output is correct
5 Correct 1 ms 4444 KB Output is correct
6 Correct 1 ms 2396 KB Output is correct
7 Correct 1 ms 2396 KB Output is correct
8 Correct 1 ms 2392 KB Output is correct
9 Correct 1 ms 2396 KB Output is correct
10 Correct 1 ms 2396 KB Output is correct
11 Correct 1 ms 2396 KB Output is correct
12 Correct 1 ms 2644 KB Output is correct
13 Correct 1 ms 2396 KB Output is correct
14 Correct 0 ms 2396 KB Output is correct
15 Correct 1 ms 2392 KB Output is correct
16 Correct 1 ms 2392 KB Output is correct
17 Correct 2 ms 2396 KB Output is correct
18 Correct 4 ms 2396 KB Output is correct
19 Correct 5 ms 2556 KB Output is correct
20 Correct 4 ms 2396 KB Output is correct
21 Correct 3 ms 2396 KB Output is correct
22 Correct 3 ms 2396 KB Output is correct
23 Correct 3 ms 2396 KB Output is correct
24 Correct 4 ms 2548 KB Output is correct
25 Correct 4 ms 2544 KB Output is correct
26 Correct 4 ms 2648 KB Output is correct
27 Correct 4 ms 2396 KB Output is correct
28 Correct 4 ms 2392 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2396 KB Output is correct
2 Correct 1 ms 2396 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2392 KB Output is correct
5 Correct 1 ms 4444 KB Output is correct
6 Correct 1 ms 2396 KB Output is correct
7 Correct 1 ms 2396 KB Output is correct
8 Correct 1 ms 2392 KB Output is correct
9 Correct 1 ms 2396 KB Output is correct
10 Correct 1 ms 2396 KB Output is correct
11 Correct 1 ms 2396 KB Output is correct
12 Correct 1 ms 2644 KB Output is correct
13 Correct 1 ms 2396 KB Output is correct
14 Correct 0 ms 2396 KB Output is correct
15 Correct 1 ms 2392 KB Output is correct
16 Correct 1 ms 2392 KB Output is correct
17 Correct 2 ms 2396 KB Output is correct
18 Correct 4 ms 2396 KB Output is correct
19 Correct 5 ms 2556 KB Output is correct
20 Correct 4 ms 2396 KB Output is correct
21 Correct 3 ms 2396 KB Output is correct
22 Correct 3 ms 2396 KB Output is correct
23 Correct 3 ms 2396 KB Output is correct
24 Correct 4 ms 2548 KB Output is correct
25 Correct 4 ms 2544 KB Output is correct
26 Correct 4 ms 2648 KB Output is correct
27 Correct 4 ms 2396 KB Output is correct
28 Correct 4 ms 2392 KB Output is correct
29 Correct 41 ms 4696 KB Output is correct
30 Correct 52 ms 4700 KB Output is correct
31 Correct 53 ms 4844 KB Output is correct
32 Correct 32 ms 4696 KB Output is correct
33 Correct 33 ms 4696 KB Output is correct
34 Correct 40 ms 4700 KB Output is correct
35 Correct 38 ms 4700 KB Output is correct
36 Correct 51 ms 4700 KB Output is correct
37 Correct 47 ms 4700 KB Output is correct
38 Correct 40 ms 4700 KB Output is correct
39 Correct 48 ms 4876 KB Output is correct
40 Correct 45 ms 4696 KB Output is correct
41 Correct 45 ms 4700 KB Output is correct
42 Correct 43 ms 5028 KB Output is correct
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2396 KB Output is correct
2 Correct 1 ms 2396 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2392 KB Output is correct
5 Correct 1 ms 4444 KB Output is correct
6 Correct 1 ms 2396 KB Output is correct
7 Correct 1 ms 2396 KB Output is correct
8 Correct 1 ms 2392 KB Output is correct
9 Correct 1 ms 2396 KB Output is correct
10 Correct 1 ms 2396 KB Output is correct
11 Correct 1 ms 2396 KB Output is correct
12 Correct 1 ms 2644 KB Output is correct
13 Correct 1 ms 2396 KB Output is correct
14 Correct 0 ms 2396 KB Output is correct
15 Correct 1 ms 2392 KB Output is correct
16 Correct 1 ms 2392 KB Output is correct
17 Correct 2 ms 2396 KB Output is correct
18 Correct 4 ms 2396 KB Output is correct
19 Correct 5 ms 2556 KB Output is correct
20 Correct 4 ms 2396 KB Output is correct
21 Correct 3 ms 2396 KB Output is correct
22 Correct 3 ms 2396 KB Output is correct
23 Correct 3 ms 2396 KB Output is correct
24 Correct 4 ms 2548 KB Output is correct
25 Correct 4 ms 2544 KB Output is correct
26 Correct 4 ms 2648 KB Output is correct
27 Correct 4 ms 2396 KB Output is correct
28 Correct 4 ms 2392 KB Output is correct
29 Correct 41 ms 4696 KB Output is correct
30 Correct 52 ms 4700 KB Output is correct
31 Correct 53 ms 4844 KB Output is correct
32 Correct 32 ms 4696 KB Output is correct
33 Correct 33 ms 4696 KB Output is correct
34 Correct 40 ms 4700 KB Output is correct
35 Correct 38 ms 4700 KB Output is correct
36 Correct 51 ms 4700 KB Output is correct
37 Correct 47 ms 4700 KB Output is correct
38 Correct 40 ms 4700 KB Output is correct
39 Correct 48 ms 4876 KB Output is correct
40 Correct 45 ms 4696 KB Output is correct
41 Correct 45 ms 4700 KB Output is correct
42 Correct 43 ms 5028 KB Output is correct
43 Correct 759 ms 7904 KB Output is correct
44 Correct 1386 ms 10320 KB Output is correct
45 Correct 1322 ms 10188 KB Output is correct
46 Correct 799 ms 9956 KB Output is correct
47 Correct 801 ms 9980 KB Output is correct
48 Correct 1053 ms 12216 KB Output is correct
49 Correct 1036 ms 12368 KB Output is correct
50 Incorrect 1306 ms 10048 KB Output isn't correct
51 Halted 0 ms 0 KB -
# Verdict Execution time Memory Grader output
1 Correct 0 ms 2396 KB Output is correct
2 Correct 1 ms 2396 KB Output is correct
3 Correct 1 ms 2396 KB Output is correct
4 Correct 1 ms 2392 KB Output is correct
5 Correct 1 ms 4444 KB Output is correct
6 Correct 1 ms 2396 KB Output is correct
7 Correct 1 ms 2396 KB Output is correct
8 Correct 1 ms 2392 KB Output is correct
9 Correct 1 ms 2396 KB Output is correct
10 Correct 1 ms 2396 KB Output is correct
11 Correct 1 ms 2396 KB Output is correct
12 Correct 1 ms 2644 KB Output is correct
13 Correct 1 ms 2396 KB Output is correct
14 Correct 0 ms 2396 KB Output is correct
15 Correct 1 ms 2392 KB Output is correct
16 Correct 1 ms 2392 KB Output is correct
17 Correct 2 ms 2396 KB Output is correct
18 Correct 4 ms 2396 KB Output is correct
19 Correct 5 ms 2556 KB Output is correct
20 Correct 4 ms 2396 KB Output is correct
21 Correct 3 ms 2396 KB Output is correct
22 Correct 3 ms 2396 KB Output is correct
23 Correct 3 ms 2396 KB Output is correct
24 Correct 4 ms 2548 KB Output is correct
25 Correct 4 ms 2544 KB Output is correct
26 Correct 4 ms 2648 KB Output is correct
27 Correct 4 ms 2396 KB Output is correct
28 Correct 4 ms 2392 KB Output is correct
29 Correct 41 ms 4696 KB Output is correct
30 Correct 52 ms 4700 KB Output is correct
31 Correct 53 ms 4844 KB Output is correct
32 Correct 32 ms 4696 KB Output is correct
33 Correct 33 ms 4696 KB Output is correct
34 Correct 40 ms 4700 KB Output is correct
35 Correct 38 ms 4700 KB Output is correct
36 Correct 51 ms 4700 KB Output is correct
37 Correct 47 ms 4700 KB Output is correct
38 Correct 40 ms 4700 KB Output is correct
39 Correct 48 ms 4876 KB Output is correct
40 Correct 45 ms 4696 KB Output is correct
41 Correct 45 ms 4700 KB Output is correct
42 Correct 43 ms 5028 KB Output is correct
43 Correct 759 ms 7904 KB Output is correct
44 Correct 1386 ms 10320 KB Output is correct
45 Correct 1322 ms 10188 KB Output is correct
46 Correct 799 ms 9956 KB Output is correct
47 Correct 801 ms 9980 KB Output is correct
48 Correct 1053 ms 12216 KB Output is correct
49 Correct 1036 ms 12368 KB Output is correct
50 Incorrect 1306 ms 10048 KB Output isn't correct
51 Halted 0 ms 0 KB -