#include <bits/stdc++.h>
using namespace std;
const int N = 5e5, mo = 998244353;
inline void add(int& a, const int& b) {
    a += b - mo;
    a += a >> 31 & mo;
}
inline int poww(int x, int y) {
    int ans = 1;
    for (; y; y >>= 1, x = 1ll * x * x % mo)
        if (y & 1)
            ans = 1ll * ans * x % mo;
    return ans;
}
struct mat {
    int a[3][3];
};
mat I;
inline mat operator*(const mat& a, const mat& b) {
    static mat d;
    static long long c[3][3];
    memset(c, 0, sizeof c);
    c[0][0] += 1ll * a.a[0][0] * b.a[0][0];
    c[0][1] += 1ll * a.a[0][0] * b.a[0][1];
    c[0][2] += 1ll * a.a[0][0] * b.a[0][2];
    c[0][0] += 1ll * a.a[0][1] * b.a[1][0];
    c[0][1] += 1ll * a.a[0][1] * b.a[1][1];
    c[0][2] += 1ll * a.a[0][1] * b.a[1][2];
    c[0][0] += 1ll * a.a[0][2] * b.a[2][0];
    c[0][1] += 1ll * a.a[0][2] * b.a[2][1];
    c[0][2] += 1ll * a.a[0][2] * b.a[2][2];
    c[1][0] += 1ll * a.a[1][0] * b.a[0][0];
    c[1][1] += 1ll * a.a[1][0] * b.a[0][1];
    c[1][2] += 1ll * a.a[1][0] * b.a[0][2];
    c[1][0] += 1ll * a.a[1][1] * b.a[1][0];
    c[1][1] += 1ll * a.a[1][1] * b.a[1][1];
    c[1][2] += 1ll * a.a[1][1] * b.a[1][2];
    c[1][0] += 1ll * a.a[1][2] * b.a[2][0];
    c[1][1] += 1ll * a.a[1][2] * b.a[2][1];
    c[1][2] += 1ll * a.a[1][2] * b.a[2][2];
    c[2][0] += 1ll * a.a[2][0] * b.a[0][0];
    c[2][1] += 1ll * a.a[2][0] * b.a[0][1];
    c[2][2] += 1ll * a.a[2][0] * b.a[0][2];
    c[2][0] += 1ll * a.a[2][1] * b.a[1][0];
    c[2][1] += 1ll * a.a[2][1] * b.a[1][1];
    c[2][2] += 1ll * a.a[2][1] * b.a[1][2];
    c[2][0] += 1ll * a.a[2][2] * b.a[2][0];
    c[2][1] += 1ll * a.a[2][2] * b.a[2][1];
    c[2][2] += 1ll * a.a[2][2] * b.a[2][2];
    d.a[0][0] = c[0][0] % mo;
    d.a[0][1] = c[0][1] % mo;
    d.a[0][2] = c[0][2] % mo;
    d.a[1][0] = c[1][0] % mo;
    d.a[1][1] = c[1][1] % mo;
    d.a[1][2] = c[1][2] % mo;
    d.a[2][0] = c[2][0] % mo;
    d.a[2][1] = c[2][1] % mo;
    d.a[2][2] = c[2][2] % mo;
    return d;
}
inline mat poww(mat a, int k) {
    mat ans = I;
    for (; k; k >>= 1, a = a * a)
        if (k & 1)
            ans = ans * a;
    return ans;
}
int n, K, ans, tcnt = 1;
struct node {
    int l, r, m, fa, o;
} t[N];
void build(int i) {
    static int l, r, fa, o;
    l = t[i].l;
    r = t[i].r;
    fa = t[i].fa;
    o = t[i].o;
    mat A;
    memset(A.a, 0, sizeof A.a);
    int v1, v2, v3, v4, v5;
    v1 = 1ll * l * (n - r + 1) % mo;
    v5 = i > 1 ? 1ll * t[fa].l * (n - t[fa].r + 1) % mo : 0;
    v3 = (1ll * n * (n + 1) / 2 - 1ll * l * (l - 1) / 2 - 1ll * (n - r) * (n - r + 1) / 2) % mo;
    v3 = (v3 + mo - v1) % mo;
    v1 = (v1 + mo - v5) % mo;
    int z = fa ? (t[fa].r - t[fa].l + 1) - (t[i].r - t[i].l + 1) : 0;
    v2 = 1ll * (z + 1) * z / 2 % mo;
    if (o && fa) {
        v2 = (v2 + 1ll * (t[fa].l - 1) * z) % mo;
    } else {
        v2 = (v2 + 1ll * (n - t[fa].r) * z) % mo;
    }
    v4 = (4ll * mo + 1ll * n * (n + 1) / 2 - v1 - v2 - v3 - v5) % mo;
    add(A.a[0][2], v1);
    add(A.a[1][2], v1);
    add(A.a[2][2], v1);
    add(A.a[0][0], v2);
    add(A.a[1][2], v2);
    add(A.a[2][2], v2);
    add(A.a[0][0], v3);
    add(A.a[1][0], v3);
    add(A.a[2][0], v3);
    add(A.a[0][0], v4);
    add(A.a[1][1], v4);
    add(A.a[2][2], v4);
    add(A.a[0][1], v5);
    add(A.a[1][1], v5);
    add(A.a[2][2], v5);
    ans = (ans + poww(A, K).a[0][2]) % mo;
    if (l == r)
        return;
    cin >> t[i].m;
    t[++tcnt] = (node){ t[i].l, t[i].m, 0, i, 0 };
    build(tcnt);
    t[++tcnt] = (node){ t[i].m + 1, t[i].r, 0, i, 1 };
    build(tcnt);
}
int main() {
    for (int i = 0; i < 3; ++i) I.a[i][i] = 1;
    ios::sync_with_stdio(0);
    cin.tie(0);
    cin >> n >> K;
    t[1].l = 1;
    t[1].r = n;
    build(1);
    ans = 1ll * ans * poww(poww(1ll * n * (n + 1) / 2 % mo, K), mo - 2) % mo;
    cout << ans << endl;
}