class Binary_Indexed_Tree_Exception(Exception): pass class Binary_Indexed_Tree(): def __init__(self,L): """calcを演算とするN項のBinary Indexed Treeを作成 calc:演算(2変数関数,群) unit:群calcの単位元(xe=ex=xを満たすe) inv:群calcの逆元(1変数関数) """ N=len(L) d=max(1,(N-1).bit_length()) k=2**d X=[None]+[0]*k self.num=k self.depth=d if L: for i in range(len(L)): p=i+1 while p<=k: X[p]=(X[p]+L[i]) p+=p&(-p) self.data=X def index(self,k,index=1): """第k要素の値を出力する. k:数列の要素 index:先頭の要素の番号 """ p=k+(1-index) return self.sum(p,p,index) def add(self,k,x,index=1,right=False): """第k要素にxを左から加え,更新を行う. k:数列の要素 x:更新後の値 index:先頭の要素の番号 right:「左から」が「右から」になる """ p=k+(1-index) while p<=self.num: if right==False: #左から self.data[p]=x+self.data[p] else: #右から self.data[p]=self.data[p]+x p+=p&(-p) def update(self,k,x,index=1,right=False): """第k要素をxに変え,更新を行う. k:数列の要素 x:更新後の値 """ a=self.index(k,index) if right==False: #左から y=x-a else: #右から y=-a+x self.add(k,y,index,right) def sum(self,From,To,index=1): """第From要素から第To要素までの総和を求める. ※From!=1を使うならば,群でなくてはならない. From:始まり To:終わり index:先頭の要素の番号 """ alpha=max(1,From+(1-index)) beta=min(self.num,To+(1-index)) if alpha==1: return self.__section(beta) else: return -self.__section(alpha-1)+self.__section(beta) def __section(self,To): S=0 x=To while x>0: S+=self.data[x] x-=x&(-x) return S def all_sum(self): return self.data[-1] #================================================ from collections import deque,defaultdict import sys input=sys.stdin.readline N=int(input()) A=list(map(int,input().split())) Mod=998244353 Count=defaultdict(int) for a in A: Count[a]+=1 S=list(set(A)) S.sort() T=len(S) D={x:i for i,x in enumerate(S)} L=[0]*T for x in D: L[D[x]]=Count[x] X=Binary_Indexed_Tree([0]*T) Y=Binary_Indexed_Tree(L) P=Binary_Indexed_Tree([0]*T) Q=Binary_Indexed_Tree([S[i]*Count[S[i]] for i in range(T)]) K=0 S_min,S_max=S[-1],S[0] for a in A: if a==S_min or a==S_max: X.add(D[a],1,0) Y.add(D[a],-1,0) P.add(D[a],a,0) Q.add(D[a],-a,0) continue alpha=X.sum(D[a]+1,T-1,0) beta =Y.sum(0,D[a]-1,0) mu=P.sum(D[a]+1,T-1,0) nu=Q.sum(0,D[a]-1,0) K+=alpha*nu+beta*mu+a*(alpha*beta) K%=Mod X.add(D[a],1,0) Y.add(D[a],-1,0) P.add(D[a],a,0) Q.add(D[a],-a,0) print(K)