#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){ cntl[i]=bit.sum(now-1)-bit.sum(idx[a[i]]); bit.add(idx[a[i]],1); } for(int i=n-1;i>=0;i--){ if(idx[a[i]]>0)cntr[i]=bit1.sum(idx[a[i]]-1); bit1.add(idx[a[i]],1); } rep(i,n){ suml[i]=bit2.sum(now-1)-bit2.sum(idx[a[i]]); suml[i]%=mod; bit2.add(idx[a[i]],a[i]); } for(int i=n-1;i>=0;i--){ if(idx[a[i]]>0)sumr[i]=bit3.sum(idx[a[i]]-1); sumr[i]%=mod; bit3.add(idx[a[i]],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<