#include #include using namespace std; #define rep(i,n) for(int i=0;i=0;i--) #define rng(i,c,n) for(int i=c;i #define _3PawRSC ios::sync_with_stdio(0),cin.tie(0) typedef long long ll; using vi=vector; using pii=pair; void print(){cout<<'\n';} template void print(const h&v,const t&...u){cout< struct modint{ ll x; modint(ll x=0):x((x%mod+mod)%mod){} modint operator-()const{return modint(-x);} modint& operator+=(const modint a){if((x+=a.x)>=mod) x-=mod; return *this;} modint& operator-=(const modint a){if((x+=mod-a.x)>=mod) x-=mod; return *this;} modint& operator*=(const modint a){(x*=a.x)%=mod; return *this;} modint operator+(const modint a)const{modint res(*this); return res+=a;} modint operator-(const modint a)const{modint res(*this); return res-=a;} modint operator*(const modint a)const{modint res(*this); return res*=a;} modint pow(ll n)const{ modint res=1,x(*this); while(n){ if(n&1)res*=x; x*=x; n>>=1; } return res; } modint inv()const{return pow(mod-2);} }; using mint=modint<998244353>; const int _n=3e5; mint fact[_n],invfact[_n]; void initfact(){ fact[0]=1; rng(i,1,_n){ fact[i]=fact[i-1]*mint(i); } invfact[_n-1]=fact[_n-1].inv(); per(i,_n-1){ invfact[i]=invfact[i+1]*mint(i+1); } } mint cnk(ll k,ll n){ return k>n?0:fact[n]*invfact[k]*invfact[n-k]; } int op(int l,int r){ return l+r; } int e(){ return 0; } void slv(){ int n; cin>>n; vi a(n); rep(i,n){ cin>>a[i]; } rep(i,n){ a[i]-=1; } vec(vi) cl(2,vi(n)); { atcoder::segtree seg(n); rep(i,n){ int x=a[i]; cl[0][i]=seg.prod(0,x); cl[1][i]=seg.prod(x,n); seg.set(x,1); } } vec(vi) cr(2,vi(n)); { atcoder::segtree seg(n); per(i,n){ int x=a[i]; cr[0][i]=seg.prod(0,x); cr[1][i]=seg.prod(x,n); seg.set(x,1); } } mint ans=0; rep(i,n){ int x=min(cl[1][i],cr[0][i]); mint wys=cnk(x,cl[1][i]+cr[0][i]); int y=min(cl[0][i],cr[1][i]); wys*=cnk(y,cl[0][i]+cr[1][i]); ans+=wys; } print(ans.x); } signed main(){ _3PawRSC; initfact(); slv(); }