#include <bits/stdc++.h>
using namespace std;
const int INF = 0x3f3f3f3f;
const long long LINF = 0x3f3f3f3f3f3f3f3f;
const double eps = 1e-12;
const double PI = acos(-1.0);
const double angcst = PI / 180.0;
const long long mod = 998244353;
long long gcd(long long a, long long b) { return b == 0 ? a : gcd(b, a % b); }
long long qmul(long long a, long long b) {
  long long r = 0;
  while (b) {
    if (b & 1) r = (r + a) % mod;
    b >>= 1;
    a = (a + a) % mod;
  }
  return r;
}
long long qpow(long long a, long long n) {
  long long r = 1;
  while (n) {
    if (n & 1) r = (r * a) % mod;
    n >>= 1;
    a = (a * a) % mod;
  }
  return r;
}
long long qpow(long long a, long long n, long long p) {
  long long r = 1;
  while (n) {
    if (n & 1) r = (r * a) % p;
    n >>= 1;
    a = (a * a) % p;
  }
  return r;
}
long long fac[5050], inv[5050], n2[5050];
long long tinv[5050];
void init() {
  fac[0] = inv[0] = n2[0] = 1;
  for (int i = (1); i <= (5000); i++)
    fac[i] = fac[i - 1] * i % mod, n2[i] = (n2[i - 1] << 1) % mod,
    tinv[i] = qpow(i, mod - 2);
  inv[5000] = qpow(fac[5000], mod - 2);
  for (int i = (4999); i >= (1); i--) inv[i] = inv[i + 1] * (i + 1) % mod;
}
long long getC(int n, int m) {
  if (m < 0 || n < 0 || n - m < 0) return 0;
  return fac[n] * inv[m] % mod * inv[n - m] % mod;
}
long long tmp[5050];
void solve() {
  init();
  long long l, n, k;
  cin >> n >> k >> l;
  for (int i = (k); i <= (n); i++)
    for (int j = (0); j <= (i); j++)
      tmp[i] = (tmp[i] + getC(i, j) * tinv[i + j + 1] % mod * (j % 2 ? -1 : 1) +
                mod) %
               mod;
  long long ans = 0;
  for (int i = (k); i <= (n); i++) {
    long long d = 0;
    for (int j = (0); j <= (n - i); j++)
      d = (d +
           n2[i + j] * getC(n - i, j) % mod * tmp[i + j] % mod *
               (j % 2 ? -1 : 1) +
           mod) %
          mod;
    ans = (ans + d * getC(n, i)) % mod;
  }
  cout << ans * l % mod << '\n';
}
int main() {
  ios::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  { solve(); }
  return 0;
}
