結果

問題 No.1248 Kth Sum
ユーザー hitonanodehitonanode
提出日時 2020-10-02 23:20:11
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 592 ms / 2,000 ms
コード長 6,864 bytes
コンパイル時間 3,003 ms
コンパイル使用メモリ 256,344 KB
実行使用メモリ 32,256 KB
最終ジャッジ日時 2024-07-17 23:50:25
合計ジャッジ時間 13,101 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,812 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 2 ms
6,944 KB
testcase_04 AC 2 ms
6,944 KB
testcase_05 AC 2 ms
6,940 KB
testcase_06 AC 2 ms
6,940 KB
testcase_07 AC 2 ms
6,944 KB
testcase_08 AC 1 ms
6,940 KB
testcase_09 AC 2 ms
6,940 KB
testcase_10 AC 1 ms
6,944 KB
testcase_11 AC 2 ms
6,944 KB
testcase_12 AC 2 ms
6,944 KB
testcase_13 AC 130 ms
21,248 KB
testcase_14 AC 377 ms
28,032 KB
testcase_15 AC 103 ms
18,048 KB
testcase_16 AC 467 ms
32,000 KB
testcase_17 AC 461 ms
32,028 KB
testcase_18 AC 442 ms
31,988 KB
testcase_19 AC 442 ms
32,128 KB
testcase_20 AC 453 ms
31,924 KB
testcase_21 AC 451 ms
31,616 KB
testcase_22 AC 362 ms
27,724 KB
testcase_23 AC 2 ms
6,940 KB
testcase_24 AC 421 ms
32,256 KB
testcase_25 AC 592 ms
32,144 KB
testcase_26 AC 433 ms
32,128 KB
testcase_27 AC 46 ms
10,296 KB
testcase_28 AC 48 ms
10,240 KB
testcase_29 AC 542 ms
32,128 KB
testcase_30 AC 319 ms
32,136 KB
testcase_31 AC 264 ms
32,128 KB
testcase_32 AC 258 ms
32,256 KB
testcase_33 AC 2 ms
6,944 KB
testcase_34 AC 531 ms
32,128 KB
testcase_35 AC 516 ms
32,128 KB
testcase_36 AC 495 ms
32,044 KB
testcase_37 AC 457 ms
31,104 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using lint = long long;
using pint = pair<int, int>;
using plint = pair<lint, lint>;
struct fast_ios {
    fast_ios() { cin.tie(nullptr), ios::sync_with_stdio(false), cout << fixed << setprecision(20); };
} fast_ios_;
#define FOR(i, begin, end) for (int i = (begin), i##_end_ = (end); i < i##_end_; i++)
#define IFOR(i, begin, end) for (int i = (end)-1, i##_begin_ = (begin); i >= i##_begin_; i--)
#define REP(i, n) FOR(i, 0, n)
#define IREP(i, n) IFOR(i, 0, n)
#define ALL(x) (x).begin(), (x).end()
//
template <typename T, typename V>
void ndarray(vector<T>& vec, const V& val, int len) { vec.assign(len, val); }
template <typename T, typename V, typename... Args>
void ndarray(vector<T>& vec, const V& val, int len, Args... args)
{
    vec.resize(len), for_each(begin(vec), end(vec), [&](T& v) { ndarray(v, val, args...); });
}
template <typename T>
bool chmax(T& m, const T q) { return m < q ? (m = q, true) : false; }
template <typename T>
bool chmin(T& m, const T q) { return m > q ? (m = q, true) : false; }
template <typename T1, typename T2>
pair<T1, T2> operator+(const pair<T1, T2>& l, const pair<T1, T2>& r) { return make_pair(l.first + r.first, l.second + r.second); }
template <typename T1, typename T2>
pair<T1, T2> operator-(const pair<T1, T2>& l, const pair<T1, T2>& r) { return make_pair(l.first - r.first, l.second - r.second); }
template <typename T>
vector<T> srtunq(vector<T> vec) { return sort(vec.begin(), vec.end()), vec.erase(unique(vec.begin(), vec.end()), vec.end()), vec; }
template <typename T>
istream& operator>>(istream& is, vector<T>& vec)
{
    return for_each(begin(vec), end(vec), [&](T& v) { is >> v; }), is;
}

// output
template <typename T, typename V>
ostream& dmpseq(ostream&, const T&, const string&, const string&, const string&);
#if __cplusplus >= 201703L
template <typename... T>
ostream& operator<<(ostream& os, const tuple<T...>& tpl)
{
    return apply([&os](auto&&... args) { ((os << args << ','), ...); }, tpl), os;
}
#endif
//
template <typename T1, typename T2>
ostream& operator<<(ostream& os, const pair<T1, T2>& p) { return os << '(' << p.first << ',' << p.second << ')'; }
template <typename T>
ostream& operator<<(ostream& os, const vector<T>& x) { return dmpseq<vector<T>, T>(os, x, "[", ",", "]"); }
template <typename T>
ostream& operator<<(ostream& os, const deque<T>& x) { return dmpseq<deque<T>, T>(os, x, "deq[", ",", "]"); }
template <typename T>
ostream& operator<<(ostream& os, const set<T>& x) { return dmpseq<set<T>, T>(os, x, "{", ",", "}"); }
template <typename T, typename TH>
ostream& operator<<(ostream& os, const unordered_set<T, TH>& x) { return dmpseq<unordered_set<T, TH>, T>(os, x, "{", ",", "}"); }
template <typename T>
ostream& operator<<(ostream& os, const multiset<T>& x) { return dmpseq<multiset<T>, T>(os, x, "{", ",", "}"); }
template <typename TK, typename T>
ostream& operator<<(ostream& os, const map<TK, T>& x) { return dmpseq<map<TK, T>, pair<TK, T>>(os, x, "{", ",", "}"); }
template <typename TK, typename T, typename TH>
ostream& operator<<(ostream& os, const unordered_map<TK, T, TH>& x) { return dmpseq<unordered_map<TK, T, TH>, pair<TK, T>>(os, x, "{", ",", "}"); }
template <typename T, typename V>
ostream& dmpseq(ostream& os, const T& seq, const string& pre, const string& sp, const string& suf)
{
    return os << pre, for_each(begin(seq), end(seq), [&](V x) { os << x << sp; }), os << suf;
}
template <typename T>
void print(const vector<T>& x) { dmpseq<vector<T>, T>(cout, x, "", " ", "\n"); }
#ifdef HITONANODE_LOCAL
#define dbg(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ") " << __FILE__ << endl
#else
#define dbg(x) {}
#endif

template <typename T>
struct BIT : std::vector<T>
{
    BIT(int len = 0) : std::vector<T>(len + 1) {}
    void reset() { fill(this->begin(), this->end(), 0); }
    void add(int pos, T v) {
        while (pos > 0 and pos < (int)this->size()) (*this)[pos] += v, pos += pos & -pos;
    }
    T sum(int pos) const { // (0, pos]
        T res = 0;
        while (pos > 0) res += (*this)[pos], pos -= pos & -pos;
        return res;
    }
    friend std::ostream &operator<<(std::ostream &os, const BIT &bit) {
        T prv = 0;
        os << '[';
        for (int i = 1; i < (int)bit.size(); i++) {
            T now = bit.sum(i);
            os << now - prv << ",";
            prv = now;
        }
        os << ']';
        return os;
    }
};

#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/pb_ds/tag_and_trait.hpp>
using namespace __gnu_pbds; // find_by_order(), order_of_key()
template <typename TK> using pbds_set = tree<TK, null_type, less<TK>, rb_tree_tag, tree_order_statistics_node_update>;
template <typename TK, typename TV> using pbds_map = tree<TK, TV, less<TK>, rb_tree_tag, tree_order_statistics_node_update>;


int main()
{
    int N, K;
    cin >> N >> K;
    vector<lint> A(N);
    cin >> A;
    vector<plint> A2I(N);
    REP(i, N) A2I[i] = plint(A[i], i);
    sort(ALL(A2I));
    vector<int> ii(N);
    REP(i, N) ii[A2I[i].second] = i + 1;
    BIT<lint> bit(N + 10);

    lint ret = A[K - 1];

    set<pint> cur_set;
    set<pint> rightset;
    pbds_set<pint> all_set;
    FOR(i, K - 1, N) all_set.insert(pint(A[i], i)), bit.add(ii[i], A[i]);
    FOR(i, K - 1, 2 * (K - 1) + 1) if (i < N) cur_set.emplace(A[i], i);
    FOR(i, 2 * K - 1, N) rightset.emplace(A[i], i);

    for (lint m = 2; m * K <= N; m++) {
        if (cur_set.empty() or rightset.empty()) break;
        dbg(m);
        dbg(rightset);
        dbg(cur_set);
        lint cum = rightset.begin()->first + cur_set.begin()->first;
        int ir = rightset.begin()->second, ic = cur_set.begin()->second;
        all_set.erase(pint(A[ir], ir));
        all_set.erase(pint(A[ic], ic));
        bit.add(ii[ir], -A[ir]);
        bit.add(ii[ic], -A[ic]);
        dbg(cum);
        dbg(bit);
        if (int(all_set.size()) < m - 2) break;
        if (int(all_set.size()) > m - 3 and m >= 3)
        {
            auto p = *all_set.find_by_order(m - 3);
            cum += bit.sum(ii[p.second]);
        }
        chmin(ret, cum);

        all_set.insert(pint(A[ir], ir));
        all_set.insert(pint(A[ic], ic));
        bit.add(ii[ir], A[ir]);
        bit.add(ii[ic], A[ic]);

        const int l = (K - 1) * (m - 1);
        FOR(i, l, l + K - 1)
        {
            if (i >= N) break;
            bit.add(ii[i], -A[i]);
            cur_set.erase(pint(A[i], i));
            all_set.erase(pint(A[i], i));
        }
        FOR(i, l + K, l + K * 2 - 1) if (i < N)
        {
            cur_set.insert(pint(A[i], i));
        }
        FOR(i, m * K - 1, (m + 1) * K - 1)
        {
            if (i >= N) break;
            rightset.erase(pint(A[i], i));
        }
    }
    dbg(cur_set);
    dbg(rightset);
    dbg(bit);
    cout << ret << '\n';
}
0