結果

問題 No.404 部分門松列
ユーザー msm1993msm1993
提出日時 2021-03-23 15:34:47
言語 C++14
(gcc 13.2.0 + boost 1.83.0)
結果
AC  
実行時間 206 ms / 2,000 ms
コード長 1,888 bytes
コンパイル時間 923 ms
コンパイル使用メモリ 82,348 KB
実行使用メモリ 17,560 KB
最終ジャッジ日時 2023-08-16 17:44:19
合計ジャッジ時間 9,718 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 1 ms
4,376 KB
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 1 ms
4,380 KB
testcase_05 AC 2 ms
4,376 KB
testcase_06 AC 1 ms
4,376 KB
testcase_07 AC 1 ms
4,376 KB
testcase_08 AC 2 ms
4,376 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 2 ms
4,380 KB
testcase_11 AC 2 ms
4,376 KB
testcase_12 AC 1 ms
4,380 KB
testcase_13 AC 91 ms
8,784 KB
testcase_14 AC 65 ms
6,764 KB
testcase_15 AC 72 ms
5,884 KB
testcase_16 AC 94 ms
12,244 KB
testcase_17 AC 134 ms
13,080 KB
testcase_18 AC 66 ms
4,696 KB
testcase_19 AC 47 ms
7,212 KB
testcase_20 AC 77 ms
9,076 KB
testcase_21 AC 127 ms
11,728 KB
testcase_22 AC 81 ms
4,736 KB
testcase_23 AC 90 ms
9,424 KB
testcase_24 AC 137 ms
9,648 KB
testcase_25 AC 206 ms
17,560 KB
testcase_26 AC 194 ms
13,944 KB
testcase_27 AC 37 ms
4,380 KB
testcase_28 AC 111 ms
7,756 KB
testcase_29 AC 144 ms
13,044 KB
testcase_30 AC 91 ms
6,212 KB
testcase_31 AC 21 ms
4,380 KB
testcase_32 AC 126 ms
13,484 KB
testcase_33 AC 132 ms
10,416 KB
testcase_34 AC 110 ms
7,072 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
int main() {
	cin.tie(0);
	ios_base::sync_with_stdio(false);
	int N;
	cin >> N;
	vector<int> A(N);
	for (int i = 0; i < N; ++i) {
		cin >> A[i];
	}
	vector<int> SA = A;
	sort(SA.begin(), SA.end());
	SA.erase(unique(SA.begin(), SA.end()), SA.end());
	for (int i = 0; i < N; ++i) {
		A[i] = lower_bound(SA.begin(), SA.end(), A[i]) - SA.begin();
	}
	int S = SA.size();
	vector<vector<int> > G(S);
	for (int i = 0; i < N; ++i) {
		G[A[i]].push_back(i);
	}
	vector<int> lb(N), rb(N);
	for (int i = 0; i < S; ++i) {
		for (int j = 0; j < G[i].size(); ++j) {
			lb[G[i][j]] = j;
			rb[G[i][j]] = G[i].size() - j - 1;
		}
	}
	vector<int> bit1(N + 1);
	vector<long long> ans(S + 1);
	vector<long long> sum(N + 1);
	for (int i = 0; i < S; ++i) {
		for (int j = 0; j < G[i].size(); ++j) {
			long long diff = 1LL * (j + 1) * (G[i].size() - j - 1) - 1LL * j * (G[i].size() - j);
			sum[G[i][j] + 1] += diff;
		}
	}
	for (int i = 0; i < N; ++i) {
		sum[i + 1] += sum[i];
		ans[A[i] + 1] -= sum[i];
	}
	int cnt = 0;
	for (int i = 0; i < S; ++i) {
		for (int j : G[i]) {
			int lx = 0;
			for (int k = j; k >= 1; k -= k & (-k)) {
				lx += bit1[k];
			}
			int ly = j - lb[j] - lx;
			int rx = cnt - lx;
			int ry = (N - j - 1) - rb[j] - rx;
			ans[i + 1] += 1LL * lx * rx + 1LL * ly * ry;
		}
		for (int j = 0; j < G[i].size(); ++j) {
			ans[i + 1] += 1LL * j * (G[i].size() - j);
		}
		cnt += G[i].size();
		for (int j : G[i]) {
			for (int k = j + 1; k <= N; k += k & (-k)) {
				++bit1[k];
			}
		}
	}
	for (int i = 0; i < S; ++i) {
		ans[i + 1] += ans[i];
	}
	int Q;
	cin >> Q;
	for (int i = 0; i < Q; ++i) {
		int L, H;
		cin >> L >> H;
		int pl = lower_bound(SA.begin(), SA.end(), L) - SA.begin();
		int pr = lower_bound(SA.begin(), SA.end(), H + 1) - SA.begin();
		cout << ans[pr] - ans[pl] << '\n';
	}
	return 0;
}
0