#include <bits/stdc++.h>
using namespace std;
#define Re register
#define Mod 998244353
#define mkp(a,b) make_pair(a,b)
typedef pair<int, int> pii;
inline int read() {
    int x = 0;
    char ch = getchar();

    while (!isdigit(ch))
        ch = getchar();

    while (isdigit(ch))
        x = x * 10 + ch - '0', ch = getchar();

    return x;
}
#define N 500010
int f[N], g[N], fac[N], inv[N], Inv[N];
inline int init(int n) {
    fac[0] = fac[1] = inv[0] = inv[1] = Inv[0] = Inv[1] = 1;

    for (int i = 2; i <= n; i++)
        fac[i] = 1LL * fac[i - 1] * i % Mod;

    for (int i = 2; i <= n; i++)
        inv[i] = 1LL * (Mod - Mod / i) * inv[Mod % i] % Mod;

    for (int i = 2; i <= n; i++)
        Inv[i] = 1LL * Inv[i - 1] * inv[i] % Mod;
}
inline int C(int n, int m) {
    if (n < m)
        return 0;

    return 1LL * fac[n] * Inv[m] % Mod * Inv[n - m] % Mod;
}
inline void Add(int &x, int y) {
    x + y < Mod ? x += y : x += y - Mod;
}
inline vector<int> GetP(int n, int k) {
    vector<int> vec(n + 1);

    for (int i = 0; i <= n / (k + 1); i++) {
        for (int j = 0; j <= i; j++) {
            int tmp = 1LL * C(n - i * k, i) * C(i, j) % Mod;
            Add(vec[i * k + j], (i - j) & 1 ? Mod - tmp : tmp);
        }
    }

    // while (!vec.empty() && !vec.back()) vec.pop_back();
    // for (int i=0;i<vec.size();i++) if (vec[i]) vec[i]=Mod-vec[i];
    return vec;
}
inline vector<int> GetPoly(int n, int k) {
    vector<int> v1 = GetP(n, k);
    int siz1 = v1.size();
    vector<int> v2 = GetP(n - k, k);
    int siz2 = v2.size();

    for (int i = 0; i < siz2; i++)
        Add(v1[i + k], Mod - v2[i]);

    while (!v1.empty() && !v1.back())
        v1.pop_back();

    // for (int i=0;i<v1.size();i++) if (v1[i]) v1[i]=Mod-v1[i];
    // for (auto i:v1) printf("%d ",i); puts("");
    return v1;
    // vector<int> vec(n+2); ++n;
    // for (Re int i=0;i*k<=n;i++)
    // for (Re int j=0;j*(k+1)<=n;j++) {
    // int tmp=1LL*C(i+j,i)*C(n-i*k-j*(k+1),i+j)%Mod;
    // Add(vec[i*k+j*(k+1)],i&1?Mod-tmp:tmp);
    // }
    // while (!vec.empty() && !vec.back()) vec.pop_back();
    // for (int i=0;i<vec.size();i++) if (vec[i]) vec[i]=Mod-vec[i];
    // for (int i=0;i<vec.size();i++)
    // printf("%d ",vec[i]);
    // puts("");
    // return vec;
}
inline int Pow(int a, int b, int p = Mod) {
    int res = 1;

    for (; b; b >>= 1, a = 1LL * a * a % p)
        if (b & 1)
            res = 1LL * res * a % p;

    return res;
}
const int G = 3;
int A[N], B[N], rev[N];
inline void NTT(int *a, int len, int f) {
    int k = 0;

    while ((1 << k) < len)
        k++;

    for (int i = 0; i < len; i++) {
        rev[i] = (rev[i >> 1] >> 1) | ((i & 1) << (k - 1));

        if (i < rev[i])
            swap(a[i], a[rev[i]]);
    }

    for (int l = 2; l <= len; l <<= 1) {
        int m = l >> 1, w = Pow(G, (Mod - 1) / l);

        if (f)
            w = Pow(w, Mod - 2);

        for (int i = 0; i < len; i += l) {
            int wk = 1;

            for (int j = 0; j < m; j++, wk = 1LL * wk * w % Mod) {
                int p = a[i + j], q = 1LL * wk * a[i + j + m] % Mod;
                a[i + j] = (p + q) % Mod, a[i + j + m] = (Mod + p - q) % Mod;
            }
        }
    }
}
inline void DFT(int *a, int len) {
    NTT(a, len, 0);
}
inline void IDFT(int *a, int len) {
    NTT(a, len, 1);
    int inv = Pow(len, Mod - 2);

    for (int i = 0; i < len; i++)
        a[i] = 1LL * a[i] * inv % Mod;
}
inline vector<int> Mul(vector<int> a, vector<int> b) {
    int n = a.size(), m = b.size();
    vector<int> ans(n + m - 1);

    for (int i = 0; i < n; i++)
        for (int j = 0; j < m; j++)
            Add(ans[i + j], 1LL * a[i]*b[j] % Mod)/* ,printf("%d %d %d %d %d\n",i,j,1LL*a[i]*b[j]%Mod) */;

    // for (auto i:ans) printf("%d ",i); puts("");
    return ans;
    // int n=a.size()+1,m=b.size()+1,len=1; while (len<n+m) len<<=1;
    // for (int i=0;i<n-1;i++) A[i]=a[i]; for (int i=n-1;i<len;i++) A[i]=0;
    // for (int i=0;i<m-1;i++) B[i]=b[i]; for (int i=m-1;i<len;i++) B[i]=0;
    // DFT(A,len),DFT(B,len); for (int i=0;i<len;i++) A[i]=1LL*A[i]*B[i]%Mod; IDFT(A,len);
    // vector<int> ans(len); for (int i=0;i<len;i++) ans[i]=A[i];
    // while (!ans.empty() && !ans.back()) ans.pop_back(); return ans;
}
vector<vector<int>> vec;
inline vector<int> Merge() {
    priority_queue<pii, vector<pii>, greater<pii>> heap;

    for (int i = 0; i < vec.size(); i++)
        heap.push(mkp(vec[i].size(), i));

    while (!heap.empty()) {
        pii vec1 = heap.top();
        heap.pop();

        if (heap.empty())
            return vec[vec1.second];

        pii vec2 = heap.top();
        heap.pop();
        // puts("Mul");
        // for (int i=0;i<vec[vec1.second].size();i++) printf("%d ",vec[vec2.second][i]); puts("");
        // for (int i=0;i<vec[vec2.second].size();i++) printf("%d ",vec[vec2.second][i]); puts("");
        vec[vec1.second] = Mul(vec[vec1.second], vec[vec2.second]);
        // puts("=");
        // for (int i=0;i<vec[vec1.second].size();i++) printf("%d ",vec[vec1.second][i]); puts("");
        heap.push(mkp(vec[vec1.second].size(), vec1.second));
    }
}
int vis[N];
int main() {
    int n = read(), k = read(), m = n << 1, cnt = 0;
    init(m);

    for (int i = 1; i <= n; i++)
        g[i] = (g[i - 1] + 1LL * n * inv[i] % Mod) % Mod;

    for (int i = 1; i <= n; i++)
        vis[read()] = true;

    for (int i = 1; i <= m; i++) {
        if (vis[i]) {
            ++cnt;
            continue;
        }

        if (cnt >= k)
            vec.push_back(GetPoly(cnt, k));

        cnt = 0;
    }

    vector<int> ans = Merge();
    int res = 0; // puts("");

    // for (int i=0;i<ans.size();i++) printf("%d ",ans[i]); puts("");
    for (int i = 0; i < ans.size(); i++)
        (res += 1LL * (Mod - ans[i]) * g[i] % Mod) %= Mod;

    return printf("%d\n", res), 0;
}