結果

問題 No.404 部分門松列
ユーザー rickythetarickytheta
提出日時 2016-08-05 21:20:25
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 1,204 ms / 2,000 ms
コード長 2,258 bytes
コンパイル時間 1,680 ms
コンパイル使用メモリ 161,788 KB
実行使用メモリ 30,472 KB
最終ジャッジ日時 2023-09-11 02:34:10
合計ジャッジ時間 13,734 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,380 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,380 KB
testcase_08 AC 2 ms
4,384 KB
testcase_09 AC 2 ms
4,376 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 1 ms
4,384 KB
testcase_12 AC 2 ms
4,380 KB
testcase_13 AC 400 ms
13,052 KB
testcase_14 AC 117 ms
6,992 KB
testcase_15 AC 138 ms
5,888 KB
testcase_16 AC 391 ms
19,648 KB
testcase_17 AC 744 ms
22,228 KB
testcase_18 AC 220 ms
6,228 KB
testcase_19 AC 190 ms
10,876 KB
testcase_20 AC 160 ms
9,540 KB
testcase_21 AC 605 ms
18,632 KB
testcase_22 AC 263 ms
6,420 KB
testcase_23 AC 178 ms
10,136 KB
testcase_24 AC 263 ms
10,144 KB
testcase_25 AC 1,204 ms
30,472 KB
testcase_26 AC 908 ms
22,020 KB
testcase_27 AC 54 ms
4,380 KB
testcase_28 AC 540 ms
12,388 KB
testcase_29 AC 668 ms
20,720 KB
testcase_30 AC 174 ms
6,352 KB
testcase_31 AC 28 ms
4,380 KB
testcase_32 AC 642 ms
23,076 KB
testcase_33 AC 683 ms
17,196 KB
testcase_34 AC 212 ms
7,208 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;
typedef vector<int> vi;
typedef vector<ll> vl;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;

typedef int _loop_int;
#define REP(i,n) for(_loop_int i=0;i<(_loop_int)(n);++i)
#define FOR(i,a,b) for(_loop_int i=(_loop_int)(a);i<(_loop_int)(b);++i)
#define FORR(i,a,b) for(_loop_int i=(_loop_int)(b)-1;i>=(_loop_int)(a);--i)

#define DEBUG(x) cout<<#x<<": "<<x<<endl
#define DEBUG_VEC(v) cout<<#v<<":";REP(i,v.size())cout<<" "<<v[i];cout<<endl
#define ALL(a) (a).begin(),(a).end()

#define CHMIN(a,b) a=min((a),(b))
#define CHMAX(a,b) a=max((a),(b))

// mod
const ll MOD = 1000000007ll;
#define FIX(a) ((a)%MOD+MOD)%MOD

// floating
typedef double Real;
const Real EPS = 1e-11;
#define EQ0(x) (abs(x)<EPS)
#define EQ(a,b) (abs(a-b)<EPS)
typedef complex<Real> P;

struct Fenwick {
  int n;
  vl bit;
  Fenwick(int _n):n(_n+1){
    bit.assign(n+1,0);
  }
  void add(int a,ll w){
    for(int x=a;x<=n;x+=x&-x)
      bit[x] += w;
  }
  ll sum(int a){
    ll r=0;
    for(int x=a;x>0;x-=x&-x)
      r+=bit[x];
    return r;
  }
};

int n;
int a[252521];

int main(){
  scanf("%d",&n);
  REP(i,n)scanf("%d",a+i);
  // 座圧
  set<int> zat;
  REP(i,n)zat.insert(a[i]);
  zat.insert(1252583025);
  map<int,int> mp;
  map<int,int> revmp;
  int it = 1;
  for(int p:zat){
    mp[p] = it;
    revmp[it] = p;
    it++;
  }
  REP(i,n){
    a[i] = mp[a[i]];
  }

  Fenwick ans = Fenwick(n+10);
  Fenwick left = Fenwick(n+10);
  Fenwick mul = Fenwick(n+10);
  Fenwick right = Fenwick(n+10);
  REP(i,n)right.add(a[i],1);
  REP(i,n){
    right.add(a[i],-1);
    mul.add(a[i],-(left.sum(a[i])-left.sum(a[i]-1)));
    ll sm = left.sum(a[i]-1) * right.sum(a[i]-1);
    sm += (left.sum(n+5)-left.sum(a[i])) * (right.sum(n+5)-right.sum(a[i]));
    sm -= mul.sum(a[i]-1);
    sm -= (mul.sum(n+5)-mul.sum(a[i]));
    // DEBUG(sm);
    ans.add(a[i],sm);
    left.add(a[i],1);
    mul.add(a[i],right.sum(a[i])-right.sum(a[i]-1));
  }

  // REP(i,n+1)DEBUG(ans.sum(i));

  int q;
  scanf("%d",&q);
  while(q--){
    int x,y;
    scanf("%d%d",&x,&y);
    x = *zat.lower_bound(x); x = mp[x];
    y = *zat.upper_bound(y); y = mp[y]-1;
    ll p = ans.sum(y) - ans.sum(x-1);
    printf("%lld\n",p);
  }
  return 0;
}
0