#include using namespace std; typedef long long int ll; typedef pair P; typedef vector VI; typedef vector VVI; #define REP(i,n) for(int i=0;i<(n);i++) #define ALL(v) v.begin(),v.end() template bool chmax(T& x, const T& y){return (x bool chmin(T& x, const T& y){return (x>y)?(x=y,true):false;}; constexpr ll MOD=998244353; constexpr ll INF=2e18; template struct Segtree{ using F=function; int n; F f; X ex; vector dat; Segtree(int n_, F f_, X ex_):f(f_), ex(ex_){ n=1; while(n_>n){n*=2;} dat.assign(2*n,ex); } void set(int i, X x){ dat[i+n-1]=x; } void build(){ for(int k=n-2;k>=0;k--) dat[k]=f(dat[2*k+1],dat[2*k+2]); } void update(int i, X x){ i+=n-1; dat[i]=x; while(i>0){ i=(i-1)/2; dat[i]=f(dat[i*2+1],dat[i*2+2]); } } X query(int a, int b){ return query_sub(a,b,0,0,n); } X query_sub(int a, int b, int k, int l, int r){ if(r<=a||b<=l) return ex; else if(a<=l&&r<=b) return dat[k]; else{ X vl=query_sub(a,b,k*2+1,l,(l+r)/2); X vr=query_sub(a,b,k*2+2,(l+r)/2,r); return f(vl,vr); } } }; long long power(long long x, long long y){ x%=MOD; long long ret=1; while(y){ if(y&1) ret=ret*x%MOD; x=x*x%MOD; y>>=1; } return ret; } long long divid(long long x, long long y){ x%=MOD; return x*power(y,MOD-2)%MOD; } long long inv(int n){ static vector a={1,1}; if((int)a.size()<=n){ for(int i=a.size();i<=n;i++) a.push_back(MOD-a[MOD%i]*(MOD/i)%MOD); } return a[n]; } long long fact(int n){ static vector a={1,1}; if((int)a.size()<=n){ for(int i=a.size();i<=n;i++) a.push_back(a.back()*i%MOD); } return a[n]; } long long finv(int n){ static vector a={1,1}; if((int)a.size()<=n){ for(int i=a.size();i<=n;i++) a.push_back(a.back()*inv(i)%MOD); } return a[n]; } long long com(int n, int k){ if(n> n >> m; VI p(m), k(m); REP(i,m) cin >> p[i] >> k[i], p[i]--, k[i]--; ll ans=0; ans=(ans+divid(fact(n-m)*(n-m)%MOD*(n-m-1)%MOD,4))%MOD; VI x(n,-1); REP(i,m) x[k[i]]=p[i]; Segtree seg(n,[](int l, int r){return l+r;},0); REP(i,m) seg.set(p[i],1); seg.build(); ll t=0; REP(i,n){ if(x[i]!=-1){ t+=seg.query(0,x[i]); seg.update(x[i],0); } } ans=(ans+t*fact(n-m)%MOD)%MOD; REP(i,m) seg.update(p[i],1); ll cnt=0; REP(i,n){ if(x[i]!=-1){ ll l=i-cnt; ll r=n-m-l; ll f=x[i]-seg.query(0,x[i]); ll g=n-x[i]-seg.query(x[i],n); ans=(ans+f*fact(n-m-1)%MOD*r%MOD)%MOD; ans=(ans+g*fact(n-m-1)%MOD*l%MOD)%MOD; cnt++; } } cout << ans << endl; return 0; }