結果

問題 No.2327 Inversion Sum
ユーザー 0214sh70214sh7
提出日時 2023-05-28 15:57:50
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 54 ms / 2,000 ms
コード長 3,959 bytes
コンパイル時間 2,257 ms
コンパイル使用メモリ 211,992 KB
実行使用メモリ 8,692 KB
最終ジャッジ日時 2023-08-27 13:28:12
合計ジャッジ時間 3,985 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge11
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 12 ms
7,144 KB
testcase_01 AC 48 ms
8,428 KB
testcase_02 AC 37 ms
8,692 KB
testcase_03 AC 10 ms
6,696 KB
testcase_04 AC 54 ms
7,988 KB
testcase_05 AC 9 ms
5,036 KB
testcase_06 AC 37 ms
7,744 KB
testcase_07 AC 21 ms
4,756 KB
testcase_08 AC 8 ms
4,380 KB
testcase_09 AC 47 ms
8,016 KB
testcase_10 AC 14 ms
4,452 KB
testcase_11 AC 5 ms
4,904 KB
testcase_12 AC 4 ms
4,616 KB
testcase_13 AC 2 ms
4,376 KB
testcase_14 AC 31 ms
5,940 KB
testcase_15 AC 53 ms
8,096 KB
testcase_16 AC 20 ms
6,956 KB
testcase_17 AC 6 ms
5,932 KB
testcase_18 AC 7 ms
4,376 KB
testcase_19 AC 17 ms
6,876 KB
testcase_20 AC 1 ms
4,376 KB
testcase_21 AC 2 ms
4,380 KB
testcase_22 AC 1 ms
4,376 KB
testcase_23 AC 2 ms
4,376 KB
testcase_24 AC 1 ms
4,376 KB
testcase_25 AC 1 ms
4,376 KB
testcase_26 AC 1 ms
4,380 KB
testcase_27 AC 2 ms
4,380 KB
testcase_28 AC 2 ms
4,376 KB
testcase_29 AC 1 ms
4,376 KB
testcase_30 AC 1 ms
4,380 KB
testcase_31 AC 2 ms
4,376 KB
testcase_32 AC 2 ms
4,380 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<ll,ll> PP;
//#define MOD 1000000007
#define MOD 998244353
#define INF 2305843009213693951
//#define INF 810114514
#define PI 3.141592653589
#define setdouble setprecision
#define REP(i,n) for(ll i=0;i<(n);++i)
#define OREP(i,n) for(ll i=1;i<=(n);++i)
#define RREP(i,n) for(ll i=(n)-1;i>=0;--i)
#define ALL(v) (v).begin(), (v).end()
#define GOODBYE do { cout << "-1" << endl; return 0; } while (false)
#define MM <<" "<<
#define Endl endl
#define debug true
#define debug2 false

class compress{
    /*
    Copyright (c) 2021 0214sh7
    https://github.com/0214sh7/library/
    */
    private:
    std::vector<int> E;
    
    public:
    void init(std::vector<long long> A){
        E.clear();
        sort(A.begin(),A.end());
        for(int i=0;i<A.size();i++){
            if(i==0 || A[i]!=A[i-1]){
                E.push_back(A[i]);
            }
        }
    }
    
    compress(std::vector<long long> A){
        init(A);
    }
    
    int size(){
        return (int)E.size();
    }
    
    int value(int x){
        if(0<=x && x<(int)E.size()){
            return E[x];
        }else{
            return 0;
        }
    }
    
    int index(int X){
        return (upper_bound(E.begin(),E.end(),X))-E.begin()-1;
    }
    
};

class Fenwick_tree{
    /*
    Copyright (c) 2021 0214sh7
    https://github.com/0214sh7/library/
    */
    private:
    std::vector<long long> BIT;
    
    public:
    
    void add(int a,long long w){
        for(int x=a;x<BIT.size();x|=(x+1)){
            BIT[x]+=w;
        }
    }
    
    void init(int n){
        BIT.clear();
        for(int i=0;i<n;i++){
            BIT.push_back(0);
        }
    }
    
    Fenwick_tree(int n){
        init(n);
    }
    
    long long sum(int a){
        long long r=0;
        for(int x=a-1;x>=0;x=(x&(x+1))-1){
            r+=BIT[x];
        }
        return r;
    }
    
    long long inversion(std::vector<long long> V){
        long long r=0;
        init(V.size());
        for(int i=0;i<V.size();i++){
            add(V[i],1);
            r+=i-sum(V[i]);
        }
        return r;
    }
};

long long inverse(long long b,long long mod){
    /*
    Copyright (c) 2021 0214sh7
    https://github.com/0214sh7/library/
    */
    long long r=1,e=mod-2;
    while(e){
        if(e&1){
            r=(r*b)%mod;
        }
        b=(b*b)%mod;
        e >>=1;
    }
    return r;
}

int main(void){
    //cin.tie(nullptr);
    //ios::sync_with_stdio(false);

    ll N,M;
    cin >> N >> M;
    vector<ll> P(M),K(M);
    REP(i,M){
        cin >> P[i] >> K[i];
        P[i]--;K[i]--;
    }

    vector<ll> Y(N,1),Z(N,-1);
    REP(i,M){
        Y[P[i]]=0;
        Z[K[i]]=P[i];
    }

    vector<ll> E = {-INF};
    REP(i,N){
        if(Y[i]==1)E.push_back(i);
    }


    vector<ll> X;
    REP(i,N){if(Z[i]!=-1)X.push_back(Z[i]);}
    compress zaatsu(X);
    REP(i,X.size()){X[i] = zaatsu.index(X[i]);}
    
    Fenwick_tree BIT(N);
    ll Ans = BIT.inversion(X)%MOD;
    REP(i,N-M){Ans = (Ans*(i+1))%MOD;}

    //cout << Ans << endl;
    vector<ll> W(N,0);
    ll sum = 0;
    REP(i,N){
        if(Z[i]!=-1){
            auto it = upper_bound(ALL(E),Z[i]);it--;
            ll u = distance(E.begin(),it);
            W[Z[i]] = u%MOD;
            sum = (sum+N-M-u)%MOD;
        }
    }


    ll num = 0;
    REP(i,N){
        if(Z[i]!=-1){
            sum = (sum+MOD-(N-M-W[Z[i]]))%MOD;
            sum = (sum+W[Z[i]])%MOD;
        }else{
            num = (num+sum)%MOD;
        }
        
    }

    //num = (num*(MOD+1)/2)%MOD;
    num = (num*inverse(N-M,MOD))%MOD;
    REP(i,N-M){num = (num*(i+1))%MOD;}
    
    Ans = (Ans+num)%MOD;

    //cout << num << endl;

    //cout << Ans << endl;

    ll n = N-M;
    n = (n*(n-1)/2)%MOD;
    n = (n*(MOD+1)/2)%MOD;
    REP(i,N-M){n = (n*(i+1))%MOD;}
    //cout << n << endl;
    Ans = (Ans+n)%MOD;

    cout << Ans << endl;


    return 0;
}
0