結果

問題 No.1300 Sum of Inversions
ユーザー kcz146kcz146
提出日時 2020-11-27 22:02:06
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 300 ms / 2,000 ms
コード長 2,867 bytes
コンパイル時間 2,399 ms
コンパイル使用メモリ 215,124 KB
実行使用メモリ 22,796 KB
最終ジャッジ日時 2024-07-26 12:45:39
合計ジャッジ時間 10,874 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
6,816 KB
testcase_01 AC 1 ms
6,812 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 224 ms
21,984 KB
testcase_04 AC 207 ms
21,856 KB
testcase_05 AC 172 ms
13,240 KB
testcase_06 AC 259 ms
22,164 KB
testcase_07 AC 252 ms
22,228 KB
testcase_08 AC 270 ms
22,412 KB
testcase_09 AC 267 ms
22,412 KB
testcase_10 AC 136 ms
12,992 KB
testcase_11 AC 138 ms
13,092 KB
testcase_12 AC 211 ms
21,932 KB
testcase_13 AC 210 ms
21,904 KB
testcase_14 AC 300 ms
22,704 KB
testcase_15 AC 277 ms
22,324 KB
testcase_16 AC 230 ms
22,068 KB
testcase_17 AC 133 ms
12,964 KB
testcase_18 AC 163 ms
13,196 KB
testcase_19 AC 186 ms
21,616 KB
testcase_20 AC 193 ms
21,816 KB
testcase_21 AC 183 ms
21,768 KB
testcase_22 AC 172 ms
13,388 KB
testcase_23 AC 253 ms
22,356 KB
testcase_24 AC 179 ms
13,460 KB
testcase_25 AC 156 ms
13,152 KB
testcase_26 AC 151 ms
13,164 KB
testcase_27 AC 157 ms
13,292 KB
testcase_28 AC 296 ms
22,512 KB
testcase_29 AC 184 ms
21,732 KB
testcase_30 AC 269 ms
22,340 KB
testcase_31 AC 168 ms
13,372 KB
testcase_32 AC 177 ms
13,516 KB
testcase_33 AC 229 ms
22,796 KB
testcase_34 AC 244 ms
22,728 KB
testcase_35 AC 205 ms
22,776 KB
testcase_36 AC 210 ms
22,756 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// ＵＴＦ−８

#include<bits/stdc++.h>
/* #include<atcoder/all> */
/* using namespace atcoder; */

using namespace std;
using ll = long long int;
using lc = complex<ll>;

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

template <typename T>
struct SegmentTree {
    using F = function<T(T, T)>;

    const T e;
    const F f;
    size_t sz;
    vector<T> tree;

    SegmentTree(size_t n, const F &f, const T &e = 0) : f(f), e(e) {
        sz = 1;
        while(sz < n) sz <<= 1;
        tree.assign(2*sz, e);
    }

    void set(typename vector<T>::iterator begin, typename vector<T>::iterator end) {
        copy(begin, end, tree.begin() + sz);
        for(size_t k=sz-1; k>0; k--)
            tree[k] = f(tree[2*k+0], tree[2*k+1]);
    }

    void update(size_t k, const T &x) {
        k += sz;
        tree[k] = x;
        while(k >>= 1)
            tree[k] = f(tree[2*k+0], tree[2*k+1]);
    }

    T query(size_t a, size_t b) const {
        b = min(b, sz);
        a = max((size_t)0, a);
        T l = e, r = e;
        for(a += sz, b += sz; a < b; a >>= 1, b >>= 1) {
            if(a & 1) l = f(l, tree[a++]);
            if(b & 1) r = f(tree[--b], r);
        }
        return f(l, r);
    }

    T operator[](const size_t k) const {
        return tree[sz + k];
    }
};


int main(void) {
    constexpr ll MOD = 0 ? 1e9+7 : 998244353;
    constexpr double PI = acos(-1);
    constexpr double eps = 1e-10;
    cout << fixed << setprecision(32);
    cin.tie(0); ios::sync_with_stdio(false);

    if(1) {
        ll n;
        cin >> n;
        vector<ll> a(n);
        for(auto &e: a) cin >> e;
        vector<ll> b = a;
        sort(begin(b), end(b));
        b.erase(unique(begin(b),end(b)), end(b));

        SegmentTree<pair<ll,ll>> st1(n+10, [&](auto a, auto b) -> pair<ll,ll>{
                auto [x, m] = a;
                auto [y, n] = b;
                return {(x+y)%MOD, (m+n)%MOD};
                }, {0ll, 0ll});
        SegmentTree<pair<ll,ll>> st2(n+10, [&](auto a, auto b) -> pair<ll,ll>{
                auto [x, m] = a;
                auto [y, n] = b;
                return {(x+y)%MOD, (m+n)%MOD};
                }, {0ll, 0ll});

        ll r = 0;
        for(ll i=0; i<n; i++) {
            ll j = lower_bound(begin(b), end(b), a[i]) - begin(b);
            auto [x, y] = st1.query(j+1, n+5);
            auto [z, w] = st2.query(j+1, n+5);
            (r += (z + ((w*a[i])%MOD)) % MOD) %= MOD;
            st2.update(j, {(get<0>(st2[j]) + x+((a[i]*y) % MOD))%MOD, (get<1>(st2[j]) + y) %MOD});
            st1.update(j, {(get<0>(st1[j]) + a[i])%MOD, (get<1>(st1[j]) + 1) %MOD});
            // cout << x << ' ' << y << ' ' << z << ' ' << w << ' ' << r << endl;
        }
        cout << r << endl;
    }

    return 0;
}
0