結果

問題 No.404 部分門松列
ユーザー r_dream0r_dream0
提出日時 2017-02-18 17:52:32
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 1,509 ms / 2,000 ms
コード長 2,398 bytes
コンパイル時間 1,386 ms
コンパイル使用メモリ 92,500 KB
実行使用メモリ 29,240 KB
最終ジャッジ日時 2023-09-11 02:48:29
合計ジャッジ時間 20,312 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,384 KB
testcase_01 AC 1 ms
4,384 KB
testcase_02 AC 1 ms
4,380 KB
testcase_03 AC 2 ms
4,384 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 2 ms
4,380 KB
testcase_07 AC 2 ms
4,384 KB
testcase_08 AC 2 ms
4,384 KB
testcase_09 AC 3 ms
4,380 KB
testcase_10 AC 2 ms
4,380 KB
testcase_11 AC 3 ms
4,380 KB
testcase_12 AC 2 ms
4,384 KB
testcase_13 AC 591 ms
12,572 KB
testcase_14 AC 464 ms
7,520 KB
testcase_15 AC 419 ms
6,172 KB
testcase_16 AC 788 ms
19,136 KB
testcase_17 AC 976 ms
21,252 KB
testcase_18 AC 373 ms
6,044 KB
testcase_19 AC 330 ms
10,444 KB
testcase_20 AC 494 ms
10,424 KB
testcase_21 AC 872 ms
18,188 KB
testcase_22 AC 460 ms
6,000 KB
testcase_23 AC 550 ms
11,016 KB
testcase_24 AC 773 ms
10,996 KB
testcase_25 AC 1,509 ms
29,240 KB
testcase_26 AC 1,347 ms
21,480 KB
testcase_27 AC 284 ms
4,384 KB
testcase_28 AC 689 ms
11,848 KB
testcase_29 AC 1,060 ms
19,944 KB
testcase_30 AC 531 ms
6,428 KB
testcase_31 AC 165 ms
4,380 KB
testcase_32 AC 929 ms
22,128 KB
testcase_33 AC 849 ms
16,572 KB
testcase_34 AC 641 ms
7,796 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cassert>
#include <algorithm>
#include <vector>
#include <map>
using namespace std;
struct BIT{
  vector<int64_t> bit;int64_t size;
  BIT(int64_t n){
    size=n;bit=vector<int64_t>(n+1);
  }
  //a[i]=x
  void add(int64_t i,int64_t x){
    i+=1;
    while(i<=size)bit[i]+=x,i+=i&-i;
  }

  int64_t get(int64_t i) {
    return sum(i + 1) - sum(i);
  }
  //a[k] (k<i)の和
  int64_t sum(int64_t i){
    int64_t s=0;
    while(i>0)s+=bit[i],i-=i&-i;
    return s;
  }
};

int N;
vector<int64_t> a;
vector<int64_t> ca;
map<int64_t, int64_t> compress_table;
map<int64_t, int64_t> answer_count;
// aの値を圧縮するテーブルを作成する
void init_compress_table() {
  ca = a;
  sort(ca.begin(), ca.end());
  ca.erase(unique(ca.begin(), ca.end()), ca.end());
  for(int i = 0; i < ca.size(); i++) {
    compress_table[ca[i]] = i;
  }
}

// aの値を圧縮する
int compress(int64_t b) {
  assert(compress_table.count(b));
  return compress_table[b];
}

int main() {
  cin >> N;
  a.resize(N);
  for(auto &num: a) cin >> num;
  init_compress_table();
  BIT left(N + 1), right(N + 1);
  for(int i = 0; i < N; i++) {
    right.add(compress(a[i]), 1);
  }
  int64_t same = 0;
  for(int i = 0; i < N; i++) {
    int cpa = compress(a[i]);
    same -= left.get(cpa) * right.get(cpa);
    right.add(compress(a[i]), -1);
    same += left.get(cpa) * right.get(cpa);
    // 小さいもの
    int64_t l_sum = left.sum(cpa);
    int64_t r_sum = right.sum(cpa);
    answer_count[a[i]] += l_sum * r_sum;
    // 大きいもの
    l_sum = left.sum(N + 1) - left.sum(cpa + 1);
    r_sum = right.sum(N + 1) - right.sum(cpa + 1);
    answer_count[a[i]] += l_sum * r_sum;
    // 同じもの(後で引く)
    l_sum = left.get(cpa);
    r_sum = right.get(cpa);
    answer_count[a[i]] += l_sum * r_sum;
    // 左右が同じものを引く
    answer_count[a[i]] -= same;
    same -= left.get(cpa) * right.get(cpa);
    left.add(compress(a[i]), 1);
    same += left.get(cpa) * right.get(cpa);
  }
  BIT answer(N + 1);
  for(auto t: answer_count) {
    answer.add(compress(t.first), t.second);
  }
  int64_t Q;
  cin >> Q;
  for(int i = 0; i < Q; i++) {
    int64_t L, H;
    cin >> L >> H;
    H ++;
    int64_t a = -answer.sum(lower_bound(ca.begin(), ca.end(), L) - ca.begin());
    a += answer.sum(lower_bound(ca.begin(), ca.end(), H) - ca.begin());
    cout << a << endl;
  }
}
0