#include using namespace std; //#define int long long #define REP(i,m,n) for(int i=(m);i<(n);i++) #define rep(i,n) REP(i,0,n) #define pb push_back #define all(a) a.begin(),a.end() #define rall(c) (c).rbegin(),(c).rend() #define mp make_pair #define endl '\n' #define vec vector #define mat vector > #define fi first #define se second typedef long long ll; typedef unsigned long long ull; typedef pair pll; typedef long double ld; typedef complex comp; const ll INF=1e9+7; const ll inf=INF*INF; const ll MOD=998244353; const ll mod=MOD; const ll MAX=200010; const double PI=acos(-1.0); //BIT閉区間 template< class T > struct BIT { vector< T > data; BIT(int sz) { data.assign(++sz, 0); } T sum(int k) { T ret = 0; for(++k; k > 0; k -= k & -k) ret += data[k]; return (ret); } void add(int k, T x) { for(++k; k < data.size(); k += k & -k) data[k] += x; } }; signed main(){ ll n;cin>>n; vectora(n); setst; rep(i,n){ cin>>a[i]; st.insert(a[i]); } ll now=0; mapidx; for(auto e:st){ idx[e]=now; now++; } BITbit(now); BITbit1(now); BITbit2(now); BITbit3(now); ll ans=0; vectorcntl(n),cntr(n),suml(n),sumr(n); rep(i,n){ ll id=idx[a[i]]; cntl[i]=bit.sum(now-1)-bit.sum(id); bit.add(id,1); } for(int i=n-1;i>=0;i--){ ll id=idx[a[i]]; if(id>0)cntr[i]=bit1.sum(id-1); bit1.add(id,1); } rep(i,n){ ll id=idx[a[i]]; suml[i]=bit2.sum(now-1)-bit2.sum(id); suml[i]%=mod; bit2.add(id,a[i]); } for(int i=n-1;i>=0;i--){ ll id=idx[a[i]]; if(id>0)sumr[i]=bit3.sum(id-1); sumr[i]%=mod; bit3.add(id,a[i]); } rep(i,n){ ans=(ans+a[i]%mod*cntl[i]%mod*cntr[i]%mod)%mod; ans=(ans+suml[i]*cntr[i]%mod)%mod; ans=(ans+sumr[i]*cntl[i]%mod)%mod; } cout<