const MOD: usize = 998244353;

#[derive(Debug, Clone)]
struct SegmentTree<F, T> {
    n: usize,
    orgn: usize,
    unit: T,
    mem: Vec<T>,
    operator: F
}

impl<F, T> SegmentTree<F, T> where
T: std::fmt::Debug + Copy,
F: Fn(T, T) -> T
{
    fn new(n: usize, unit: T, operator: F) -> Self {
        let mut size = 1usize;
        while size < n { size *= 2; }
        Self {
            n: size,
            orgn: n,
            unit: unit,
            mem: vec![unit; size*2-1],
            operator: operator,
        }
    }

    fn find_closed(&self, l: usize, r: usize) -> T {
        self.find(l, r+1)
    }

    fn find(&self, l: usize, r: usize) -> T {
        self._find(l, r, 0, 0, self.n)
    }

    fn set(&mut self, pos: usize, val: T) {
        let mut pos = pos + self.n - 1;
        self.mem[pos] = val;
        while pos > 0 {
            pos = self._parent(pos);
            self.mem[pos] = (self.operator)(self.mem[self._left(pos)], self.mem[self._right(pos)]);
        }
    }

    fn get(&mut self, pos: usize) -> T {
        self.mem[pos + self.n - 1]
    }

    fn _find(&self, a: usize, b: usize, nodenum: usize, l: usize, r: usize) -> T {
        if r <= a || b <= l {
            return self.unit;
        } else if a <= l && r <= b {
            return self.mem[nodenum];
        }
        let left = self._find(a, b, self._left(nodenum), l, (l+r)/2);
        let right = self._find(a, b, self._right(nodenum), (l+r)/2, r);
        (self.operator)(left, right)
    }
    fn _left(&self, k: usize) -> usize {
        2*k+1
    }
    fn _right(&self, k: usize) -> usize {
        2*k+2
    }
    fn _parent(&self, k: usize) -> usize {
        (k-1)/2
    }
}

fn main() {
    let mut n = String::new();
    std::io::stdin().read_line(&mut n).ok();
    let n: usize = n.trim().parse().unwrap();
    let mut a = String::new();
    std::io::stdin().read_line(&mut a).ok();
    let a: Vec<usize> = a.trim().split_whitespace().map(|s| s.parse().unwrap()).collect();

    let mut a = (0..n).map(|i| (a[i], i)).collect::<Vec<_>>();
    a.sort();
    let mut minvalseg = SegmentTree::new(n, 0usize, |x, y| x + y);
    let mut maxvalseg = SegmentTree::new(n, 0usize, |x, y| x + y);
    let mut mincntseg = SegmentTree::new(n, 0usize, |x, y| x + y);
    let mut maxcntseg = SegmentTree::new(n, 0usize, |x, y| x + y);
    let mut min_cnts = vec![0usize; n];
    let mut max_cnts = vec![0usize; n];
    let mut min_vals = vec![0usize; n];
    let mut max_vals = vec![0usize; n];
    for &(v, i) in a.iter() {
        min_cnts[i] = mincntseg.find(i, n);
        mincntseg.set(i, 1);
        min_vals[i] = minvalseg.find(i, n) % MOD;
        minvalseg.set(i, v);
    }
    a.reverse();
    for &(v, i) in a.iter() {
        max_cnts[i] = maxcntseg.find(0, i);
        maxcntseg.set(i, 1);
        max_vals[i] = maxvalseg.find(0, i) % MOD;
        maxvalseg.set(i, v);
    }
    let mut result = 0usize;
    for &(v, i) in a.iter() {
        result += min_cnts[i] * max_vals[i] + max_cnts[i] * min_vals[i] + v * min_cnts[i] % MOD * max_cnts[i] % MOD;
        result %= MOD;
    }
    println!("{}", result);
}