#include #include #include #include using namespace std; const long long MOD = 998244353; //B template struct BIT{//1_Indexed int n; vector bit; BIT(int n_):n(n_+1),bit(n,0){} T sum(int a){ T ret = 0; for(int i = a; i > 0; i -= i & -i) ret += bit[i]; return ret; } void add(int a,int w){ for(int i = a; i <= n; i += i & -i)bit[i] += w; } T query(int r,int l){ return sum(l-1)-sum(r-1); } }; //E int main(){ int n;cin>>n; vector A(n),B(n); for(int i = 0; n > i; i++)cin>>A[i],B[i] = A[i]; sort(B.begin(),B.end()); B.erase(unique(B.begin(),B.end()),B.end()); map C; for(int i = 0; B.size() > i; i++){ C[B[i]] = i+1; } BIT X(B.size()+1),Y(B.size()+1),Z(B.size()+1),W(B.size()+1); for(int i = 0; n > i; i++){ Y.add(C[A[i]],1); W.add(C[A[i]],A[i]); } long long ans = 0; for(int i = 0; n > i; i++){ auto x = X.query(C[A[i]]+1,B.size()+1); auto y = Y.query(1,C[A[i]]); auto z = Z.query(C[A[i]]+1,B.size()+1); auto w = W.query(1,C[A[i]]); ans = (ans + z*y + ((A[i]*x)%MOD)*y + w*x)%MOD; X.add(C[A[i]],1); Y.add(C[A[i]],-1); Z.add(C[A[i]],A[i]); W.add(C[A[i]],-A[i]); } cout << ans << endl; }