제출 #1128621

#제출 시각아이디문제언어결과실행 시간메모리
1128621vladiliusCake 3 (JOI19_cake3)C++20
100 / 100
1138 ms17028 KiB
#include <bits/stdc++.h> using namespace std; using ll = long long; using pii = pair<int, int>; using pli = pair<ll, int>; #define pb push_back #define ff first #define ss second const ll infm = -1e18; struct ST{ vector<pli> t; int n, m; ST(int ns, int ms){ n = ns; m = ms; t.resize(4 * n); } void upd(int v, int tl, int tr, int& p, int& x){ if (tl == tr){ t[v] = {x, (x > 0)}; return; } int tm = (tl + tr) / 2, vv = 2 * v; if (p <= tm){ upd(vv, tl, tm, p, x); } else { upd(vv + 1, tm + 1, tr, p, x); } t[v].ff = t[vv].ff + t[vv + 1].ff; t[v].ss = t[vv].ss + t[vv + 1].ss; } void upd(int p, int x){ upd(1, 1, n, p, x); } ll get(int v, int tl, int tr, int k){ if (!k) return 0; if (t[v].ss <= k) return t[v].ff; int tm = (tl + tr) / 2, vv = 2 * v; if (t[vv + 1].ss >= k){ return get(vv + 1, tm + 1, tr, k); } return get(vv, tl, tm, k - t[vv + 1].ss) + t[vv + 1].ff; } ll get(){ return get(1, 1, n, m - 2); } }; int main(){ ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0); int n, M; cin>>n>>M; vector<pii> all; for (int i = 1; i <= n; i++){ int x, y; cin>>x>>y; all.pb({y, x}); } sort(all.begin(), all.end()); vector<int> a(n + 1), b(n + 1); for (int i = 1; i <= n; i++){ a[i] = all[i - 1].ss; b[i] = all[i - 1].ff; } all.clear(); for (int i = 1; i <= n; i++){ all.pb({a[i], i}); } sort(all.begin(), all.end()); vector<int> p(n + 1); for (int i = 0; i < n; i++){ p[all[i].ss] = i + 1; } ST T(n, M); auto f = [&](int l, int r){ if ((r - l + 1) < M) return infm; return a[l] + a[r] - 2 * (b[r] - b[l]) + T.get(); }; auto f1 = [&](int l, int r){ if ((r - l + 1) < M) return infm; ll out = a[l] + a[r] - 2 * (b[r] - b[l]); multiset<int> st; for (int i = l + 1; i < r; i++){ st.insert(a[i]); } int k = M - 2; while (k--){ auto it = prev(st.end()); out += *it; st.erase(it); } return out; }; ll out = infm; int tl = 1, tr = 0; function<void(int, int, int, int)> solve = [&](int l, int r, int l1, int r1){ if (l > r) return; int m = (l + r) / 2; pli opt = {infm, 0}; for (int i = max(l1, m + M - 1); i <= r1; i++){ while (tr < (i - 1)){ tr++; T.upd(p[tr], a[tr]); } while (tl > (m + 1)){ tl--; T.upd(p[tl], a[tl]); } while (tr > (i - 1)){ T.upd(p[tr], 0); tr--; } while (tl < (m + 1)){ T.upd(p[tl], 0); tl++; } opt = max(opt, {f(m, i), i}); } out = max(out, opt.ff); solve(l, m - 1, l1, opt.ss); solve(m + 1, r, opt.ss, r1); }; solve(1, n - M + 1, 1, n); cout<<out<<"\n"; }
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...
#Verdict Execution timeMemoryGrader output
Fetching results...