結果

問題 No.924 紲星
ユーザー betrue12betrue12
提出日時 2019-11-08 23:15:46
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 2,478 bytes
コンパイル時間 2,240 ms
コンパイル使用メモリ 189,008 KB
実行使用メモリ 26,624 KB
最終ジャッジ日時 2024-09-15 02:22:50
合計ジャッジ時間 8,332 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
5,248 KB
testcase_01 AC 2 ms
5,376 KB
testcase_02 AC 2 ms
5,376 KB
testcase_03 AC 3 ms
5,376 KB
testcase_04 AC 3 ms
5,376 KB
testcase_05 AC 3 ms
5,376 KB
testcase_06 AC 3 ms
5,376 KB
testcase_07 AC 2 ms
5,376 KB
testcase_08 TLE -
testcase_09 -- -
testcase_10 -- -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
testcase_17 -- -
testcase_18 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

template<typename T>
struct BIT {
    int n;
    vector<T> dat;

    BIT(int n=0){
        initialize(n);
    }

    void initialize(int nin){
        n = nin;
        dat.resize(n, 0);
    }

    T sum(int i){
        T s = 0;
        while(i >= 0){
            s += dat[i];
            i = (i & (i+1)) - 1;
        }
        return s;
    }

    T sum_between(int i, int j){
        return sum(j) - sum(i-1);
    }

    void plus(int i, T x){
        while(i < n){
            dat[i] += x;
            i |= i+1;
        }
    }

    // a[0]+...+a[ret] >= x
    int lower_bound(T x){
        int ret = -1;
        int k = 1;
        while(2*k <= n) k <<= 1;
        for( ;k>0; k>>=1){
            if(ret+k < n && dat[ret+k] < x){
                x -= dat[ret+k];
                ret += k;
            }
        }
        return ret + 1;
    }
};

int main(){
    int N, Q;
    cin >> N >> Q;
    vector<int> A(N);
    for(int i=0; i<N; i++){
        scanf("%d", &A[i]);
    }
    vector<int> cmp = A;
    sort(cmp.begin(), cmp.end());
    cmp.erase(unique(cmp.begin(), cmp.end()), cmp.end());
    int sz = cmp.size();
    vector<int> Ai(N);
    for(int i=0; i<N; i++) Ai[i] = lower_bound(cmp.begin(), cmp.end(), A[i]) - cmp.begin();

    const int B = 450;
    vector<vector<int>> q_rli[B];
    for(int i=0; i<Q; i++){
        int l, r;
        scanf("%d %d", &l, &r);
        l--; r--;
        q_rli[l/B].push_back({r, l, i});
    }

    BIT<int> bit_num(sz);
    BIT<int64_t> bit_sum(sz);
    int L = 0, R = -1;
    auto add = [&](int i, int d){
        bit_num.plus(Ai[i], d);
        bit_sum.plus(Ai[i], d*A[i]);
    };
    auto adjust = [&](int L2, int R2){
        while(L2 < L) L--, add(L, 1);
        while(R < R2) R++, add(R, 1);
        while(L < L2) add(L, -1), L++;
        while(R2 < R) add(R, -1), R--;
    };

    vector<int64_t> ans(Q);
    for(int b=0; b<B; b++){
        auto& qs = q_rli[b];
        sort(qs.begin(), qs.end());
        for(auto& q : qs){
            int r = q[0], l = q[1], i = q[2];
            adjust(l, r);
            int mi = bit_num.lower_bound((r-l+2)/2);
            int num1 = bit_num.sum_between(0, mi);
            int num2 = r-l+1 - num1;
            int64_t med = cmp[mi];
            ans[i] += num1 * med - bit_sum.sum_between(0, mi);
            ans[i] += bit_sum.sum_between(mi+1, sz-1) - num2 * med;
        }
    }
    for(int64_t a : ans) printf("%lld\n", a);
    return 0;
}
0