//@formatter:off #include #define rep(i,n) for (int i = 0; i < int(n); ++i) #define rrep(i,n) for (int i = int(n)-1; i >= 0; i--) #define rep2(i,s,n) for (int i = int(s); i < int(n); ++i) #define all(a) a.begin(),a.end() #define rall(a) a.rbegin(),a.rend() #define pb push_back #define eb emplace_back #define vi vector #define vvi vector> #define vl vector #define vvl vector> #define vd vector #define vvd vector> #define vs vector #define vc vector #define vvc vector> #define vb vector #define vvb vector> #define vp vector

#define vvp vector> using namespace std; using ll = long long; using P = pair; using LP = pair; template istream& operator>>(istream &is,pair &p) { return is >> p.first >> p.second; } template ostream& operator<<(ostream &os,const pair &p) { return os<<'{'< istream& operator>>(istream &is,vector &v) { for(T &t:v){is>>t;} return is; } template ostream& operator<<(ostream &os,const vector &v) { os<<'[';rep(i,v.size())os< bool chmin(T& a,T b) {if(a > b){a = b; return true;} return false;} template bool chmax(T& a,T b) {if(a < b){a = b; return true;} return false;} const int inf = 1001001001; const ll linf = 1001001001001001001; //@formatter:on //constexpr int mod = 1000000007; constexpr int mod = 998244353; struct mint { ll x; constexpr mint(ll x = 0) : x((x % mod + mod) % mod) {} constexpr mint operator-() const { return mint(-x); } constexpr mint &operator+=(const mint &a) { if ((x += a.x) >= mod) x -= mod; return *this; } constexpr mint &operator++() { return *this += mint(1); } constexpr mint &operator-=(const mint &a) { if ((x += mod - a.x) >= mod) x -= mod; return *this; } constexpr mint &operator--() { return *this -= mint(1); } constexpr mint &operator*=(const mint &a) { (x *= a.x) %= mod; return *this; } constexpr mint operator+(const mint &a) const { mint res(*this); return res += a; } constexpr mint operator-(const mint &a) const { mint res(*this); return res -= a; } constexpr mint operator*(const mint &a) const { mint res(*this); return res *= a; } constexpr mint pow(ll t) const { mint res = mint(1), a(*this); while (t > 0) { if (t & 1) res *= a; t >>= 1; a *= a; } return res; } // for prime mod constexpr mint inv() const { return pow(mod - 2); } constexpr mint &operator/=(const mint &a) { return *this *= a.inv(); } constexpr mint operator/(const mint &a) const { mint res(*this); return res /= a; } }; ostream &operator<<(ostream &os, const mint &a) { return os << a.x; } bool operator==(const mint &a, const mint &b) { return a.x == b.x; } bool operator!=(const mint &a, const mint &b) { return a.x != b.x; } template class segtree { using T = typename Monoid::value_type; int n; vector val; public: constexpr segtree(int _n) { n = 1; while (n < _n) n *= 2; val = vector(2 * n, Monoid::identity); } constexpr segtree(vector init) { int _n = init.size(); n = 1; while (n < _n) n *= 2; val = vector(2 * n, Monoid::identity); rep(i, _n) val[i + n] = init[i]; rrep(i, n) val[i] = Monoid::operate(val[i * 2], val[i * 2 + 1]); } // segment [l,r) T fold(int l, int r) { l += n, r += n; T fold_l = Monoid::identity; T fold_r = Monoid::identity; while (l < r) { if (l & 1) fold_l = Monoid::operate(fold_l, val[l++]); if (r & 1) fold_r = Monoid::operate(val[--r], fold_r); l >>= 1, r >>= 1; } return Monoid::operate(fold_l, fold_r); } template void update(int i, const F &f) { i += n; val[i] = f(val[i]); while (i > 1) { i >>= 1; val[i] = Monoid::operate(val[i * 2], val[i * 2 + 1]); } } }; class Node { public: using value_type = mint; value_type value; Node(value_type value) : value(value) {} static constexpr value_type identity = 0; static constexpr value_type operate(const value_type &l, const value_type &r) { return l + r; } }; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr); int n; cin >> n; vi a(n); cin >> a; int mx = 0; rep(i, n - 1) chmax(mx, a[i + 1] - a[i]); mint ans = 0; mint pre = 0; rep2(i, mx, 2500) { segtree st(n); st.update(n - 1, [](auto) { return 1; }); rrep(j, n - 2) { auto it = upper_bound(all(a), a[j] + i); mint sum = st.fold(j + 2, it - a.begin()); st.update(j + 1, [&](auto) { return sum; }); } ans += (st.fold(1, n) - pre) * 2 * i; pre = st.fold(1, n); } cout << ans << endl; }