結果
問題 | No.2327 Inversion Sum |
ユーザー | pockyny |
提出日時 | 2023-05-28 14:30:13 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 123 ms / 2,000 ms |
コード長 | 2,125 bytes |
コンパイル時間 | 777 ms |
コンパイル使用メモリ | 77,240 KB |
実行使用メモリ | 8,328 KB |
最終ジャッジ日時 | 2023-08-27 09:59:07 |
合計ジャッジ時間 | 3,507 ms |
ジャッジサーバーID (参考情報) |
judge13 / judge12 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 83 ms
8,328 KB |
testcase_01 | AC | 123 ms
8,048 KB |
testcase_02 | AC | 110 ms
8,140 KB |
testcase_03 | AC | 69 ms
7,300 KB |
testcase_04 | AC | 114 ms
7,200 KB |
testcase_05 | AC | 37 ms
5,120 KB |
testcase_06 | AC | 106 ms
8,328 KB |
testcase_07 | AC | 42 ms
4,868 KB |
testcase_08 | AC | 15 ms
4,380 KB |
testcase_09 | AC | 113 ms
7,480 KB |
testcase_10 | AC | 28 ms
4,380 KB |
testcase_11 | AC | 40 ms
5,760 KB |
testcase_12 | AC | 28 ms
4,912 KB |
testcase_13 | AC | 10 ms
4,380 KB |
testcase_14 | AC | 65 ms
5,440 KB |
testcase_15 | AC | 107 ms
6,928 KB |
testcase_16 | AC | 81 ms
7,356 KB |
testcase_17 | AC | 54 ms
6,516 KB |
testcase_18 | AC | 16 ms
4,380 KB |
testcase_19 | AC | 76 ms
7,220 KB |
testcase_20 | AC | 1 ms
4,376 KB |
testcase_21 | AC | 1 ms
4,380 KB |
testcase_22 | AC | 2 ms
4,376 KB |
testcase_23 | AC | 2 ms
4,376 KB |
testcase_24 | AC | 1 ms
4,380 KB |
testcase_25 | AC | 2 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 | 2 ms
4,380 KB |
testcase_30 | AC | 2 ms
4,376 KB |
testcase_31 | AC | 1 ms
4,376 KB |
testcase_32 | AC | 1 ms
4,376 KB |
ソースコード
#include <iostream> #include <vector> using namespace std; typedef long long ll; ll INF = 1000000000000000; struct segment_tree{ int n; vector<ll> seg; ll e(){return 0;} ll op(ll a,ll b){return (a + b);} segment_tree(){} segment_tree(int sz){n = sz; seg.resize(2*n);} segment_tree(vector<ll> v){ n = v.size(); seg.resize(2*n); for(int i=0;i<n;i++) seg[i + n] = v[i]; for(int i=n - 1;i>0;i--) seg[i] = op(seg[i<<1],seg[i<<1|1]); } void init(int sz){n = sz; seg.resize(2*n);} void update(int p,ll val){ for(seg[p += n] = val;p>1;p>>=1) seg[p>>1] = op(seg[p],seg[p^1]); } ll get(int p){return seg[p + n];} ll query(int l,int r){ ll res = e(); for(l += n,r += n; l<r;l>>=1,r>>=1){ if(l&1) res = op(res,seg[l++]); if(r&1) res = op(res,seg[--r]); } return res; } void clear(){for(int i=0;i<2*n;i++) seg[i] = 0;} }; ll mod = 998244353; ll pw(ll a,ll x){ ll ret = 1; while(x){ if(x&1) (ret *= a) %= mod; (a *= a) %= mod; x /= 2; } return ret; } int p[200010]; int main(){ int i,n,m; cin >> n >> m; for(i=0;i<n;i++) p[i] = -1; for(i=0;i<m;i++){ int a,b; cin >> a >> b; a--; b--; p[b] = a; } ll c = 0,ans = 0; segment_tree seg(n),seg_clear(n); for(i=0;i<n;i++) seg_clear.update(i,1); for(i=0;i<n;i++){ if(p[i]!=-1) seg_clear.update(p[i],0); } segment_tree seg2(n); for(i=0;i<n;i++){ if(p[i]!=-1){ seg.update(p[i],1); (ans += seg.query(p[i] + 1,n)) %= mod; ll num = seg_clear.query(0,n); (ans += seg_clear.query(p[i],n)*c%mod*pw(num,mod - 2)%mod) %= mod; seg2.update(p[i],num - c); }else{ c++; } } ll num = seg_clear.query(0,n); (ans += (num*(num - 1)/2)%mod*(mod + 1)/2%mod) %= mod; for(i=0;i<n;i++){ if(seg_clear.get(i)){ (ans += seg2.query(i,n)*pw(num,mod - 2)) %= mod; } } for(i=1;i<=num;i++) (ans *= i) %= mod; cout << ans << "\n"; }