#include <bits/stdc++.h>
using namespace std;

#define all(x) begin(x), end(x)
#ifdef CKISEKI
#define safe cerr << __PRETTY_FUNCTION__ << " line " << __LINE__ << " safe\n"
#define debug(a...) debug_(#a, a)
#define orange(a...) orange_(#a, a)
void debug_(auto s, auto ...a) {
  cerr << "\e[1;32m(" << s << ") = (";
  int f = 0;
  (..., (cerr << (f++ ? ", " : "") << a));
  cerr << ")\e[0m\n";
}
#include <experimental/iterator>
void orange_(auto s, auto L, auto R) {
  cerr << "\e[1;33m[ " << s << " ] = [ ";
  using namespace experimental;
  copy(L, R, make_ostream_joiner(cerr, ", "));
  cerr << " ]\e[0m\n";
}
#else
#define safe ((void)0)
#define debug(...) safe
#define orange(...) safe
#endif

using llf = long double;

const int maxn = 505;

llf g[maxn][maxn][maxn];
llf fac[maxn * 2];
llf pw[maxn * 2];

int main() {
  cin.tie(nullptr)->sync_with_stdio(false);

  int n, d, r;
  cin >> n >> d >> r;

  fac[0] = 1;
  for (int i = 1; i < maxn * 2; i++)
    fac[i] = fac[i - 1] * i;
  debug((long double)fac[maxn - 1]);
  pw[0] = 1;
  pw[1] = 1;
  {
    llf x = fac[n];
    x = powl(x, 1.0 / llf(n));
    pw[1] *= x;
  }
  {
    llf x = fac[n - 1];
    x = powl(x, 1.0 / llf(n));
    pw[1] *= x;
  }
  {
    llf x = fac[d];
    x = powl(x, 1.0 / llf(n));
    pw[1] *= x;
  }
  {
    llf x = fac[n + d - 1];
    x = powl(x, 1.0 / llf(n));
    pw[1] /= x;
  }
  debug(pw[1]);
  for (int i = 1; i <= n; i++)
  {
    //pw[i] = //powl(fac[n-1] * fac[d] / fac[n+d-1] * fac[n], i / llf(n));
    pw[i] = pw[i - 1] * pw[1];
  }
  orange(pw, pw + 500);

  for (int i = 0; i <= n; i++) {
    g[i][0][0] = pw[i] / fac[i];
  }
  for (int i = 0; i <= n; i++) {
    for (int j = 1; j <= d; j++) {
      for (int k = 0; k <= d; k++) {
        for (int t = 0; t <= i && j * t <= k; t++) {
          g[i][j][k] += g[i-t][j-1][k - j*t] / fac[t] * pw[t];
        }
      }
    }
  }
  
  llf sum = 0;

  for (int i = 0; i <= n; i++) {
    for (int j = 1; j <= d; j++) {
      for (int k = 0; k <= d - i * j; k++) {
        sum += min(i, r) * g[i][d - j][k] * g[n-i][j-1][d - i * j - k];
      }
      // orange(g[i][j], g[i][j] + 10);
    }
  }

  debug((long double)sum);

  // sum *= fac[n-1];

  // sum *= fac[d];
  // sum /= fac[n+d-1];

  sum += r;

  cout << fixed << setprecision(20);
  cout << (long double)sum << '\n';

  return 0;
}

/*
3 3 2
(fac[maxn - 1]) = (7.77944e+1144)
[ g[i][j], g[i][j] + 10 ] = [ 1, 0, 0, 0, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 0, 0, 0, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 0, 0, 0, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 0, 0, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 1, 0, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 1, 1, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 0.5, 0, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 1.5, 1, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 1.5, 2, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 0.5, 0.166667, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 1.5, 1.16667, 0, 0, 0, 0, 0, 0 ]
[ g[i][j], g[i][j] + 10 ] = [ 1, 1, 1.5, 2.16667, 0, 0, 0, 0, 0, 0 ]
(sum) = (8)
6.80000000000000000017
*/