結果

問題 No.404 部分門松列
ユーザー msm1993
提出日時 2021-03-23 15:34:47
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 200 ms / 2,000 ms
コード長 1,888 bytes
コンパイル時間 818 ms
コンパイル使用メモリ 82,544 KB
実行使用メモリ 17,668 KB
最終ジャッジ日時 2024-11-25 07:22:30
合計ジャッジ時間 8,586 ms
ジャッジサーバーID
(参考情報)
judge5 / judge1
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 4
other AC * 31
権限があれば一括ダウンロードができます

ソースコード

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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0