結果

問題 No.1300 Sum of Inversions
コンテスト
ユーザー Example0911
提出日時 2020-11-29 11:04:31
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,676 bytes
コンパイル時間 2,872 ms
コンパイル使用メモリ 202,624 KB
最終ジャッジ日時 2025-01-16 10:07:06
ジャッジサーバーID
（参考情報） 		judge4 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 33 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
		nowans += A[i] * cnt1*cnt2;
		nowans %= mod;
		ans += nowans;
		ans %= mod;
		leftsum.add(zA[i], A[i]);
		leftcnt.add(zA[i], 1);
	}
	cout << ans << endl;
}
0