結果

問題 No.1921 Range LthKth Query
ユーザー りあん
提出日時 2022-05-01 22:39:17
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 3,856 ms / 4,000 ms
コード長 3,173 bytes
コンパイル時間 5,114 ms
コンパイル使用メモリ 282,256 KB
実行使用メモリ 61,172 KB
最終ジャッジ日時 2024-07-01 01:21:36
合計ジャッジ時間 46,671 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20 TLE * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC target("avx2")
#pragma GCC optimize("O3")
// #pragma GCC optimize("unroll-loops")
#include<bits/stdc++.h>
using namespace std;

using P = pair<int, int>;
const int M = 1000000007;


template<typename T> class BIT {
private:
    int n,m; vector<vector<T> > bit;
public:
    // (i, j) に val を加算する
    void add(int i, int j, T val){
        for(int i_ = i+1; i_ < n; i_ += i_ & -i_)
            for(int j_ = j+1; j_ < m; j_ += j_ & -j_)
                bit[i_][j_] += val;
    }
    // [0,i]×[0,j]の範囲の和を求める
    T sum(int i, int j){
        T s = 0;
        for(int i_ = i+1; i_ > 0; i_ -= i_ & -i_)
            for(int j_ = j+1; j_ > 0; j_ -= j_ & -j_)
                s += bit[i_][j_];
        return s;
    }
    // [lx, rx)×[ly, ry)の範囲の和を求める
    T sum(int lx, int rx, int ly, int ry){
        return sum(rx-1, ry-1) - sum(lx-1, ry-1) - sum(rx-1, ly-1) + sum(lx-1, ly-1);
    }
    BIT(int sz1, int sz2){
        n = sz1 + 1, m = sz2 + 1;
        bit.resize(n, vector<T>(m, 0));
    }
    BIT(vector<vector<T> >& v){
        n = (int)v.size() + 1, m = (int)v[0].size() + 1;
        bit.resize(n, vector<T>(m, 0));
        for(int i = 0; i < n - 1; i++)
            for(int j = 0; j < m - 1; j++)
                add(i, j, v[i][j]);
    }
    void print(){
        for(int i = 0; i < n; i++){
            for(int j = 0; j < m; j++){
                cout<< sum(i-1, i, j-1, j) << " ";
            }
            cout << "\n";
        }
    }
    void print_sum(){
        for(int i = 0; i < n; i++){
            for(int j = 0; j < m; j++){
                cout<< sum(i-1, j-1) << " ";
            }
            cout << "\n";
        }
    }
};

int main() {
    cin.tie(0);
    ios::sync_with_stdio(0);

    int n, k, l;
    cin >> n >> k >> l;
    vector<int> a(n);
    for (int i = 0; i < n; ++i) {
        cin >> a[i];
    }
    vector<vector<P>> ps(n + 1);
    for (int i = 0; i < n; ++i) {
        priority_queue<int> lower;
        for (int j = i; j < n; ++j) {
            lower.push(a[j]);
            if ((int)lower.size() > k) lower.pop();
            if ((int)lower.size() == k) {
                ps[lower.top()].emplace_back(i, j);
            }
        }
    }
    int q;
    cin >> q;
    vector<int> ls(q), rs(q);
    for (int i = 0; i < q; ++i) {
        cin >> ls[i] >> rs[i];
        --ls[i];
    }

    const int MAXLOG = 12;
    vector<int> ng(q, 0), ok(q, n + 1);
    for (int _ = 0; _ < MAXLOG; ++_) {
        vector<vector<int>> qs(n + 1);
        for (int i = 0; i < q; ++i) {
            int m = (ng[i] + ok[i]) / 2;
            qs[m].push_back(i);
        }
        BIT<int> bit(n + 1, n + 1);
        for (int i = 0; i < n + 1; ++i) {
            for (auto& p : ps[i]) {
                bit.add(p.first, p.second, 1);
            }
            for (auto& j : qs[i]) {
                int cnt = bit.sum(ls[j], rs[j], ls[j], rs[j]);
                if (cnt >= l) {
                    ok[j] = i;
                }
                else {
                    ng[j] = i;
                }
            }
        }
    }
    for (auto& i : ok) {
        cout << i << '\n';
    }
    return 0;
}
0