#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++) #define chmin(x, y) (x) = min((x), (y)) #define chmax(x, y) (x) = max((x), (y)) #define all(x) (x).begin(),(x).end() #define inf 1e18 #define mod 998244353 using namespace std; typedef long long llint; typedef long long ll; typedef pair P; struct BIT{ int size; vector bit; BIT(){size = 0;} BIT(int s){ size = s; bit.resize(size+1); init(); } void init(){ for(int i = 1; i <= size; i++) bit[i] = 0; } llint query(int i){ llint ret = 0; while(i > 0){ ret += bit[i], ret %= mod; i -= i&(-i); } return ret; } void add(int i, llint x){ while(i <= size){ bit[i] += x, bit[i] %= mod; i += i&(-i); } } }; ll n; ll a[200005]; vector comp; ll lnum[200005], lsum[200005]; ll rnum[200005], rsum[200005]; BIT bit(200005), bit2(200005); int main(void) { ios::sync_with_stdio(0); cin.tie(0); cin >> n; rep(i, 1, n){ cin >> a[i]; comp.push_back(a[i]); } sort(all(comp)); comp.erase(unique(all(comp)), comp.end()); rep(i, 1, n) a[i] = lower_bound(all(comp), a[i]) - comp.begin() + 1; rep(i, 1, n){ lnum[i] = bit.query(n) - bit.query(a[i]); lsum[i] = bit2.query(n) - bit2.query(a[i]) + mod, lsum[i] %= mod; bit.add(a[i], 1), bit2.add(a[i], comp[a[i]-1]); } bit.init(), bit2.init(); for(int i = n; i >= 1; i--){ rnum[i] = bit.query(a[i]-1); rsum[i] = bit2.query(a[i]-1) + mod, rsum[i] %= mod; bit.add(a[i], 1), bit2.add(a[i], comp[a[i]-1]); } ll ans = 0; rep(j, 1, n){ ans += comp[a[j]-1] * lnum[j] % mod * rnum[j] % mod, ans %= mod; ans += lsum[j] * rnum[j] % mod, ans %= mod; ans += rsum[j] * lnum[j] % mod, ans %= mod; } cout << ans << endl; return 0; }