結果

問題 No.404 部分門松列
ユーザー r_dream0r_dream0
提出日時 2017-02-18 17:52:32
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,454 ms / 2,000 ms
コード長 2,398 bytes
コンパイル時間 1,490 ms
コンパイル使用メモリ 94,428 KB
実行使用メモリ 29,312 KB
最終ジャッジ日時 2024-06-28 17:18:44
合計ジャッジ時間 8,087 ms
ジャッジサーバーID
（参考情報） 		judge5 / 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 2 ms
5,376 KB
testcase_04 AC 2 ms
5,376 KB
testcase_05 AC 2 ms
5,376 KB
testcase_06 AC 3 ms
5,376 KB
testcase_07 AC 3 ms
5,376 KB
testcase_08 AC 2 ms
5,376 KB
testcase_09 AC 2 ms
5,376 KB
testcase_10 AC 2 ms
5,376 KB
testcase_11 AC 2 ms
5,376 KB
testcase_12 AC 2 ms
5,376 KB
testcase_13 AC 588 ms
12,800 KB
testcase_14 AC 476 ms
7,552 KB
testcase_15 AC 421 ms
6,528 KB
testcase_16 AC 721 ms
19,328 KB
testcase_17 AC 945 ms
21,376 KB
testcase_18 AC 376 ms
6,144 KB
testcase_19 AC 315 ms
10,752 KB
testcase_20 AC 500 ms
10,496 KB
testcase_21 AC 855 ms
18,304 KB
testcase_22 AC 461 ms
6,144 KB
testcase_23 AC 596 ms
11,008 KB
testcase_24 AC 805 ms
11,264 KB
testcase_25 AC 1,454 ms
29,312 KB
testcase_26 AC 1,279 ms
21,760 KB
testcase_27 AC 287 ms
5,376 KB
testcase_28 AC 660 ms
11,904 KB
testcase_29 AC 1,001 ms
20,224 KB
testcase_30 AC 531 ms
6,784 KB
testcase_31 AC 167 ms
5,376 KB
testcase_32 AC 905 ms
22,272 KB
testcase_33 AC 844 ms
16,640 KB
testcase_34 AC 639 ms
8,064 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