#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; using ll = long long; using P = pair; constexpr int INF = 1001001001; // constexpr int mod = 1000000007; constexpr int mod = 998244353; template inline bool chmax(T& x, T y){ if(x < y){ x = y; return true; } return false; } template inline bool chmin(T& x, T y){ if(x > y){ x = y; return true; } return false; } struct mint { int x; mint() : x(0) {} mint(int64_t y) : x(y >= 0 ? y % mod : (mod - (-y) % mod) % mod) {} mint& operator+=(const mint& p){ if((x += p.x) >= mod) x -= mod; return *this; } mint& operator-=(const mint& p){ if((x -= p.x) < 0) x += mod; return *this; } mint& operator*=(const mint& p){ x = (int)(1LL * x * p.x % mod); return *this; } mint& operator/=(const mint& p){ *this *= p.inverse(); return *this; } mint operator-() const { return mint(-x); } mint operator+(const mint& p) const { return mint(*this) += p; } mint operator-(const mint& p) const { return mint(*this) -= p; } mint operator*(const mint& p) const { return mint(*this) *= p; } mint operator/(const mint& p) const { return mint(*this) /= p; } bool operator==(const mint& p) const { return x == p.x; } bool operator!=(const mint& p) const { return x != p.x; } mint pow(int64_t n) const { mint res = 1, mul = x; while(n > 0){ if(n & 1) res *= mul; mul *= mul; n >>= 1; } return res; } mint inverse() const { return pow(mod - 2); } friend ostream& operator<<(ostream& os, const mint& p){ return os << p.x; } friend istream& operator>>(istream& is, mint& p){ int64_t val; is >> val; p = mint(val); return is; } }; template struct BinaryIndexedTree{ vector data; int sz; BinaryIndexedTree(int n){ sz = n + 2; data.assign(sz, 0); } void add(int i, T x){ ++i; while(i < sz){ data[i] += x; i += i & -i; } } // [0, i] の区間和 T sum(int i){ T res = 0; ++i; while(i > 0){ res += data[i]; i -= i & -i; } return res; } // [l, r) の区間和 T sum(int l, int r){ return sum(r - 1) - sum(l - 1); } // val 以上を満たす 0-indexed の位置を返す int lower_bound(T val){ // if(val <= 0) return -1; int res = 0; int n = 1; while(n < sz) n <<= 1; for(int k = n >> 1; k > 0; k >>= 1){ if(res + k < sz && data[res + k] < val){ val -= data[res + k]; res += k; } } return res; // 1-indexed であれば、res + 1 を返す } }; int main(){ ios::sync_with_stdio(false); cin.tie(nullptr); int N; cin >> N; vector

dat; for(int i = 0; i < N; ++i){ int A; cin >> A; dat.emplace_back(A, i); } sort(dat.begin(), dat.end()); BinaryIndexedTree bit(N), bit2(N); BinaryIndexedTree bit3(N); BinaryIndexedTree bit4(N); vector sum(N); vector cnt(N); for(auto [A, i] : dat){ cnt[i] = bit3.sum(i, N); sum[i] = bit.sum(i, N) + mint(A) * cnt[i]; bit.add(i, A); bit3.add(i, 1); } mint ans = 0; for(auto [A, i] : dat){ ans += bit2.sum(i, N) + mint(A) * bit4.sum(i, N); bit2.add(i, sum[i]); bit4.add(i, cnt[i]); } cout << ans << endl; }