結果

問題 No.1300 Sum of Inversions
ユーザー rodearodea
提出日時 2020-12-06 16:53:29
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 773 ms / 2,000 ms
コード長 3,035 bytes
コンパイル時間 1,722 ms
コンパイル使用メモリ 132,672 KB
実行使用メモリ 45,372 KB
最終ジャッジ日時 2024-09-17 13:08:25
合計ジャッジ時間 20,470 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 2 ms
6,944 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 594 ms
39,372 KB
testcase_04 AC 568 ms
38,784 KB
testcase_05 AC 437 ms
27,372 KB
testcase_06 AC 645 ms
41,756 KB
testcase_07 AC 609 ms
40,780 KB
testcase_08 AC 709 ms
43,076 KB
testcase_09 AC 712 ms
42,808 KB
testcase_10 AC 365 ms
24,676 KB
testcase_11 AC 409 ms
24,916 KB
testcase_12 AC 600 ms
39,036 KB
testcase_13 AC 559 ms
38,436 KB
testcase_14 AC 773 ms
45,040 KB
testcase_15 AC 687 ms
42,632 KB
testcase_16 AC 587 ms
39,956 KB
testcase_17 AC 340 ms
24,356 KB
testcase_18 AC 402 ms
26,272 KB
testcase_19 AC 498 ms
36,828 KB
testcase_20 AC 530 ms
37,248 KB
testcase_21 AC 499 ms
37,252 KB
testcase_22 AC 431 ms
27,484 KB
testcase_23 AC 650 ms
41,780 KB
testcase_24 AC 451 ms
27,868 KB
testcase_25 AC 361 ms
25,768 KB
testcase_26 AC 371 ms
25,464 KB
testcase_27 AC 420 ms
27,012 KB
testcase_28 AC 716 ms
43,612 KB
testcase_29 AC 500 ms
37,136 KB
testcase_30 AC 726 ms
42,940 KB
testcase_31 AC 448 ms
27,596 KB
testcase_32 AC 480 ms
28,124 KB
testcase_33 AC 143 ms
26,620 KB
testcase_34 AC 157 ms
26,680 KB
testcase_35 AC 457 ms
45,324 KB
testcase_36 AC 460 ms
45,372 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize("O3")
#include <iostream>
#include <iomanip>
#include <cstdio>
#include <string>
#include <cstring>
#include <deque>
#include <list>
#include <queue>
#include <stack>
#include <vector>
#include <utility>
#include <algorithm>
#include <map>
#include <set>
#include <complex>
#include <cmath>
#include <limits>
#include <cfloat>
#include <climits>
#include <ctime>
#include <cassert>
#include <numeric>
#include <fstream>
#include <functional>
#include <bitset>
using namespace std;

using ll = long long;
using P = pair<int, int>;
using T = tuple<int, int, int>;

template <class T> inline T chmax(T &a, const T b) {return a = (a < b) ? b : a;}
template <class T> inline T chmin(T &a, const T b) {return a = (a > b) ? b : a;}

constexpr int mod = 998244353;
constexpr int MOD = 1e9 + 7;
constexpr int inf = 1e9;
constexpr long long INF = 1e18;

#define all(a) (a).begin(), (a).end()

int dx[] = {1, 0, -1, 0};
int dy[] = {0, 1, 0, -1};

// Range Sum Query
class SegmentTree{
public:
    int n;
    vector<ll> node;

    SegmentTree(vector<ll> v) {
        int sz = v.size();
        n = 1; while(n < sz) n *= 2;
        node.resize(2*n-1, 0);
        for(int i=0; i<sz; i++) node[i+n-1] = v[i];
        for(int i=n-2; i>=0; i--) node[i] = node[i*2+1] + node[i*2+2];
    }

    void add(int x, ll val){
        x += (n - 1);
        node[x] += val;
        while(0 < x){
            x = (x - 1) / 2;
            node[x] = node[2*x+1] + node[2*x+2];
        }
    }

    ll getsum(int a, int b, int k=0, int l=0, int r=-1){
        if(r < 0) r = n;
        if(b <= l || r <= a) return 0;
        if(a <= l && r <= b) return node[k];

        ll vl = getsum(a, b, 2*k+1, l, (l+r)/2);
        ll vr = getsum(a, b, 2*k+2, (l+r)/2, r);

        return vl + vr;
    }
};

int main(){
    cin.tie(0);
    ios::sync_with_stdio(false);

    int n; cin>>n;
    vector<int> a(n);
    for(int i=0; i<n; i++) cin>>a[i];

    set<int> st;
    for(int i=0; i<n; i++) st.emplace(a[i]);

    map<int, int> num2idx;
    int idx = 0;
    for(auto i:st){
        num2idx[i] = idx++;
    }

    SegmentTree clseg(vector<ll>(n, 0)), slseg(vector<ll>(n, 0));
    vector<ll> cntl(n, 0), suml(n, 0);
    for(int i=0; i<n; i++){
        cntl[i] = clseg.getsum(num2idx[a[i]] + 1, n);
        suml[i] = slseg.getsum(num2idx[a[i]] + 1, n);
        clseg.add(num2idx[a[i]], 1);
        slseg.add(num2idx[a[i]], a[i]);
    }

    SegmentTree crseg(vector<ll>(n, 0)), srseg(vector<ll>(n, 0));
    vector<ll> cntr(n, 0), sumr(n, 0);
    for(int i=n-1; i>=0; i--){
        cntr[i] = crseg.getsum(0, num2idx[a[i]]);
        sumr[i] = srseg.getsum(0, num2idx[a[i]]);
        crseg.add(num2idx[a[i]], 1);
        srseg.add(num2idx[a[i]], a[i]);
    }

    ll ans = 0;
    for(int i=1; i+1<n; i++){
        ll L = (suml[i] % mod) * cntr[i] % mod;
        ll R = (sumr[i] % mod) * cntl[i] % mod;
        ll mid = a[i] * cntl[i] % mod * cntr[i] % mod;
        ans += (L + mid + R) % mod;
        ans %= mod;
    }

    cout << ans << endl;

    return 0;
}
0