結果
問題 | No.1300 Sum of Inversions |
ユーザー | Example0911 |
提出日時 | 2020-11-29 11:13:28 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 229 ms / 2,000 ms |
コード長 | 1,735 bytes |
コンパイル時間 | 2,349 ms |
コンパイル使用メモリ | 202,556 KB |
最終ジャッジ日時 | 2025-01-16 10:08:03 |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 34 |
ソースコード
#include "bits/stdc++.h" using namespace std; using ll = long long; ll mod = 998244353; template <typename T> struct BIT { int n; vector<T>bit; BIT(int n_) : n(n_ + 1), bit(n, 0) {} void add(int i, T x) { for (int idx = i; idx < n; idx += (idx & -idx)) { bit[idx] += x; } } T sum(int i) { T s(0); for (int idx = i; idx > 0; idx -= (idx & -idx)) { s += bit[idx]; } return s; } }; signed main() { int N; cin >> N; // 座圧パート vector<ll>A(N), cA(N), zA(N); for (int i = 0; i < N; i++)cin >> A[i]; reverse(A.begin(), A.end()); for (int i = 0; i < N; i++)cA[i] = A[i]; sort(A.begin(), A.end()); A.erase(unique(A.begin(), A.end()), A.end()); for (int i = 0; i < N; i++) { zA[i] = lower_bound(A.begin(), A.end(), cA[i]) - A.begin(); zA[i]++; } A.resize(N); for (int i = 0; i < N; i++) { A[i] = cA[i]; } BIT<ll> leftsum(200010), rightsum(200010); for (int i = 0; i < N; i++) { rightsum.add(zA[i], A[i]); } BIT<ll>leftcnt(200010), rightcnt(200010); for (int i = 0; i < N; i++) { rightcnt.add(zA[i], 1); } // setとBITを用いて答えを求める ll ans = 0; for (int i = 0; i < N; i++) { rightsum.add(zA[i], -A[i]); rightcnt.add(zA[i], -1); ll nowans = 0; ll sum1 = leftsum.sum(zA[i] - 1); ll sum2 = rightsum.sum(200009) - rightsum.sum(zA[i]); ll cnt1 = rightcnt.sum(200009) - rightcnt.sum(zA[i]); ll cnt2 = leftcnt.sum(zA[i] - 1); sum1 %= mod; sum2 %= mod; nowans += sum1 * cnt1; nowans %= mod; nowans += sum2 * cnt2; nowans %= mod; ll tmp = A[i]; tmp *= cnt1; tmp %= mod; tmp *= cnt2; tmp %= mod; nowans += tmp; nowans %= mod; ans += nowans; ans %= mod; leftsum.add(zA[i], A[i]); leftcnt.add(zA[i], 1); } cout << ans << endl; }