結果
| 問題 |
No.2327 Inversion Sum
|
| コンテスト | |
| ユーザー |
shobonvip
|
| 提出日時 | 2023-05-28 14:41:13 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 61 ms / 2,000 ms |
| コード長 | 1,368 bytes |
| コンパイル時間 | 4,353 ms |
| コンパイル使用メモリ | 253,760 KB |
| 最終ジャッジ日時 | 2025-02-13 11:53:27 |
|
ジャッジサーバーID (参考情報) |
judge3 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 30 |
ソースコード
#include<bits/stdc++.h>
#include<atcoder/all>
using namespace std;
using namespace atcoder;
typedef modint998244353 mint;
typedef long long ll;
//defmodfact
const int COMinitMAX = 998244;
mint fact[COMinitMAX+1], factinv[COMinitMAX+1];
void modfact(){
fact[0] = 1;
for (int i=1; i<=COMinitMAX; i++){
fact[i] = fact[i-1] * i;
}
factinv[COMinitMAX] = fact[COMinitMAX].inv();
for (int i=COMinitMAX-1; i>=0; i--){
factinv[i] = factinv[i+1] * (i+1);
}
}
mint cmb(int a, int b){
if (a<b || b<0) return mint(0);
return fact[a]*factinv[b]*factinv[a-b];
}
//--------
int main(){
modfact();
int n, m; cin >> n >> m;
vector<int> kettei(n, -1);
for (int i=0; i<m; i++){
int p, k; cin >> p >> k;
p--;
k--;
kettei[k] = p;
}
fenwick_tree<mint> f(n);
for (int i=0; i<n; i++){
if (kettei[i] >= 0){
f.add(kettei[i], -1);
}
f.add(i, 1);
}
fenwick_tree<mint> g(n);
mint fw = 0;
mint ans = 0;
mint kunum = 0;
for (int i=0; i<n; i++){
if (kettei[i] >= 0){
ans += fact[n-m] * g.sum(kettei[i], n);
ans += fact[n-m-1] * kunum * f.sum(kettei[i], n);
fw += f.sum(0, kettei[i]);
g.add(kettei[i], 1);
}else{
ans += fact[n-m-1] * fw;
ans += fact[n-m-2] * mint(n-m) * mint(n-m-1) / 2 * kunum;
kunum += 1;
}
//cout << i << " " << fw.val() << " " << kunum.val() << " " << ans.val() << endl;
}
cout << ans.val() << endl;
}
shobonvip