結果

問題 No.1300 Sum of Inversions
ユーザー ysystem7ysystem7
提出日時 2020-11-27 22:07:38
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 937 ms / 2,000 ms
コード長 4,888 bytes
コンパイル時間 3,279 ms
コンパイル使用メモリ 237,420 KB
実行使用メモリ 43,020 KB
最終ジャッジ日時 2024-07-26 12:57:24
合計ジャッジ時間 25,586 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,816 KB
testcase_01 AC 2 ms
6,940 KB
testcase_02 AC 2 ms
6,944 KB
testcase_03 AC 718 ms
37,764 KB
testcase_04 AC 694 ms
37,148 KB
testcase_05 AC 536 ms
25,984 KB
testcase_06 AC 808 ms
39,824 KB
testcase_07 AC 769 ms
38,948 KB
testcase_08 AC 878 ms
40,944 KB
testcase_09 AC 871 ms
40,848 KB
testcase_10 AC 433 ms
23,552 KB
testcase_11 AC 429 ms
23,680 KB
testcase_12 AC 714 ms
37,312 KB
testcase_13 AC 674 ms
36,756 KB
testcase_14 AC 937 ms
42,712 KB
testcase_15 AC 854 ms
40,556 KB
testcase_16 AC 733 ms
38,100 KB
testcase_17 AC 413 ms
23,040 KB
testcase_18 AC 479 ms
24,832 KB
testcase_19 AC 620 ms
35,340 KB
testcase_20 AC 727 ms
35,516 KB
testcase_21 AC 597 ms
35,592 KB
testcase_22 AC 520 ms
25,984 KB
testcase_23 AC 824 ms
39,660 KB
testcase_24 AC 548 ms
26,496 KB
testcase_25 AC 454 ms
24,320 KB
testcase_26 AC 447 ms
24,192 KB
testcase_27 AC 511 ms
25,472 KB
testcase_28 AC 875 ms
41,456 KB
testcase_29 AC 622 ms
35,512 KB
testcase_30 AC 864 ms
40,812 KB
testcase_31 AC 517 ms
26,112 KB
testcase_32 AC 571 ms
26,624 KB
testcase_33 AC 18 ms
6,940 KB
testcase_34 AC 29 ms
6,940 KB
testcase_35 AC 502 ms
43,020 KB
testcase_36 AC 507 ms
42,936 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#pragma GCC optimize("Ofast")
#include <bits/stdc++.h>
#define rep(i,n) for(int i=0;i<n;i++)
#define cinf(n,x) for(int i=0;i<(n);i++)cin>>x[i];
#define ft first
#define sc second
#define pb push_back
#define lb lower_bound
#define ub upper_bound
#define all(v) (v).begin(),(v).end()
#define LB(a,x) lb(all(a),x)-a.begin()
#define UB(a,x) ub(all(a),x)-a.begin()
//#define mod 1000000007
#define mod 998244353
#define FS fixed<<setprecision(15)
using namespace std;
typedef long long ll;
const double pi=3.141592653589793;
template<class T> using V=vector<T>;
using Graph = vector<vector<int>>;
using P=pair<ll,ll>;
typedef unsigned long long ull;
typedef long double ldouble;
template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; }
template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; }
template<class T> inline void out(T a){ cout << a << '\n'; }
void YN(bool ok){if(ok) cout << "Yes" << endl; else cout << "No" << endl;}
//void YN(bool ok){if(ok) cout << "YES" << endl; else cout << "NO" << endl;}


const ll INF=1e18;
const int mx=200005;

class Segtree{
public:
    ll N=1;
    V<ll> dat;
    explicit Segtree(ll sz){
        while(N<sz) N*=2;
        dat.reserve(2*N);
        rep(i,2*N)dat[i]=0;
    }

    // Range Minimum Query -------------------------

    void min_init(){
        for(int i=0;i<N;i++)min_update(i,INF);
    }

    void min_update(ll k,ll a){
        k+=N-1;
        dat[k]=a;
        while(k>0){
            k=(k-1)/2;
            dat[k]=min(dat[2*k+1],dat[2*k+2]);
        }
    }

    ll min_query(ll a,ll b,ll k,ll l,ll r){
        if(r<=a||b<=l)return INF;
        if(a<=l&&r<=b)return dat[k];
        ll m=(l+r)/2;
        ll u=min_query(a,b,2*k+1,l,m);
        ll v=min_query(a,b,2*k+2,m,r);
        return min(u,v);
    }

    ll get_min(ll l,ll r){
        return min_query(l,r,0,0,N);
    }

    //Range Maximum Query ---------------------------

    void max_init(){
        for(int i=0;i<N;i++)max_update(i,-INF);
    }

    void max_update(ll k,ll a){
        k+=N-1;
        dat[k]=a;
        while(k>0){
            k=(k-1)/2;
            dat[k]=max(dat[2*k+1],dat[2*k+2]);
        }
    }

    ll max_query(ll a,ll b,ll k,ll l,ll r){
        if(r<=a||b<=l)return -INF;
        if(a<=l&&r<=b)return dat[k];
        ll m=(l+r)/2;
        ll u=max_query(a,b,2*k+1,l,m);
        ll v=max_query(a,b,2*k+2,m,r);
        return max(u,v);
    }

    ll get_max(ll l,ll r){
        return max_query(l,r,0,0,N);
    }

    //Range Sum Query -----------------------------

    void sum_update(ll k,ll a){
        k+=N-1;
        dat[k]+=a;
        dat[k]%=mod;
        while(k>0){
            k=(k-1)/2;
            dat[k]=(dat[2*k+1]+dat[2*k+2])%mod;
        }
    }
    
    ll sum_query(ll a,ll b,ll k,ll l,ll r){
        if(r<=a||b<=l)return 0LL;
        if(a<=l&&r<=b)return dat[k];
        ll m=(l+r)/2;
        ll u=sum_query(a,b,2*k+1,l,m);
        ll v=sum_query(a,b,2*k+2,m,r);
        return (u+v)%mod;
    }

    ll get_sum(ll l,ll r){
        return sum_query(l,r,0,0,N);
    }

    //その他----------------------------------------
    /*
    ll unit(){
        return ;//単位元
    }
    ll calc(ll a,ll b){
        return ;//演算
    }
    void update(ll k,ll a){
        k+=N-1;
        dat[k]=a;
        while(k>0){
            k=(k-1)/2;
            dat[k]=calc(dat[2*k+1],dat[2*k+2]);
        }
    }
    ll query(ll a,ll b,ll k,ll l,ll r){
        if(r<=a||b<=l)return unit();
        if(a<=l&&r<=b)return dat[k];
        ll m=(l+r)/2;
        ll u=query(a,b,2*k+1,l,m);
        ll v=query(a,b,2*k+2,m,r);
        return calc(u,v);
    }
    ll get_val(ll l,ll r){
        return query(l,r,0,0,N);
    }
    */
};

int main(){
    //オーバーフローは大丈夫ですか？？
    cin.tie(0);ios::sync_with_stdio(false);
    ll n;
    cin>>n;
    V<ll> a(n);
    set<ll> st;
    //Segtree slsum(n),slcnt(n)
    rep(i,n){
        cin>>a[i];
        st.insert(a[i]);
    }
    int sz=(int)st.size();
    Segtree slsum(sz),slcnt(sz),srsum(sz),srcnt(sz);
    //V<ll> v;
    int pos=0;
    map<ll,ll> mp;
    for(ll p:st){
        //v.pb(p);
        mp[p]=pos;
        pos++;
    }
    slsum.sum_update(mp[a[0]],a[0]);
    slcnt.sum_update(mp[a[0]],1);
    for(int i=1;i<n;i++){
        srsum.sum_update(mp[a[i]],a[i]);
        srcnt.sum_update(mp[a[i]],1);
    }

    ll ans=0;
    for(int i=1;i<n-1;i++){
        int p=mp[a[i]];
        ll s=slcnt.get_sum(p+1,sz);
        ll t=srcnt.get_sum(0,p);
        ll ss=slsum.get_sum(p+1,sz);
        ll tt=srsum.get_sum(0,p);
        ans=(ans+a[i]*s%mod*t%mod)%mod;
        ans=(ans+ss*t%mod)%mod;
        ans=(ans+s*tt%mod)%mod;
        slcnt.sum_update(p,1LL);
        slsum.sum_update(p,(ll)a[i]);
        srcnt.sum_update(p,-1LL);
        srsum.sum_update(p,(ll)-a[i]);
    }
    out(ans);
}
0