結果

問題 No.2327 Inversion Sum
ユーザー keisuke6
提出日時 2023-05-28 15:36:40
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,768 bytes
コンパイル時間 1,892 ms
コンパイル使用メモリ 204,916 KB
最終ジャッジ日時 2025-02-13 13:37:08
ジャッジサーバーID
（参考情報） 		judge3 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 29 TLE * 1
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
long long modinv(long long a, long long m) {
    long long b = m, u = 1, v = 0;
    while (b) {
        long long t = a / b;
        a -= t * b; swap(a, b);
        u -= t * v; swap(u, v);
    }
    u %= m;
    if (u < 0) u += m;
    return u;
}
ll merge_count(vector<int> &a) {
    int n = a.size();
    if (n <= 1) { return 0; }
 
    ll cnt = 0;
    vector<int> b(a.begin(), a.begin()+n/2);
    vector<int> c(a.begin()+n/2, a.end());
 
    cnt += merge_count(b);
    cnt += merge_count(c);
 
    int ai = 0, bi = 0, ci = 0;
 
    while (ai < n) {
        if ( bi < b.size() && (ci == c.size() || b[bi] <= c[ci]) ) {
            a[ai++] = b[bi++];
        } else {
            cnt += n / 2 - bi;
            a[ai++] = c[ci++];
        }
    }
    return cnt;
}
signed main(){
  int N,M;
  cin>>N>>M;
  set<int> s;
  vector<int> A(N,-1);
  for(int i=0;i<M;i++){
  	int p,k;
  	cin>>p>>k;
  	p--;
  	k--;
  	A[k] = p;
  	s.insert(p);
  }
  vector<int> S = {}, T = {};
  for(int i=0;i<N;i++){
  	if(!s.count(i)) T.push_back(i);
  	if(A[i] != -1) S.push_back(A[i]);
  }
  int mod = 998244353;
  int a = 1, b = (T.size())*(T.size()-1)/2%mod, c = 0;
  for(int i=0;i<T.size();i++){
  	a *= (i+1);
  	a %= mod;
  }
  a = a*merge_count(S)%mod;
  for(int i=3;i<=T.size();i++){
  	b *= i;
  	b %= mod;
  }
  int x = 1, count = 0;
  for(int i=0;i<T.size()-1;i++){
  	x *= (i+1);
  	x %= mod;
  }
  int n = T.size();
  for(int i=0;i<N;i++){
  	if(A[i] == -1){
  		count++;
  		continue;
  	}
  	int ind = lower_bound(T.begin(),T.end(),A[i])-T.begin();
  	c += (count*(n-ind)+(n-count)*ind)%mod*x%mod;
  	c %= mod;
  }
  //cout<<a<<' '<<b<<' '<<c<<endl;
  cout<<(a+b+c+10*mod)%mod<<endl;
}
0