#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using P = pair<int, ll>;
using T = tuple<int, int, int>;
#define al(a) a.begin(), a.end()
#define ral(a) a.rbegin(), a.rend()
#define sz(a) (int)a.size()
#define rep(i, n) for (int i = 0; i < (n); ++i)
#define rrep(i, n) for (int i = 1; i <= (n); ++i)
#define drep(i, n) for (int i = (n)-1; i >= 0; --i)
#define db(a, b) cout << #a << ": " << a << " " << #b << ": " << b << endl;

#include <atcoder/modint>
using namespace atcoder;
using mint = modint1000000007;

int main() {
  ll n, s, K;
  cin >> n >> s >> K;
  s -= n * (n - 1) / 2 * K;

  map<P, mint> memo;
  auto dfs = [&](auto &f, int n, ll s) -> mint {
    if (n == 0 && s == 0) return 1;
    if (n == 0 && s != 0) return 0;
    if (s < 0) return 0;
    if (memo.find({n, s}) != memo.end()) return memo[{n, s}];

    mint res;
    for (int i = 0; i * n <= s; ++i) {
      ll tmp = n * i;
      res += f(f, n - 1, s - tmp);
    }

    return memo[{n, s}] = res;
  };
  mint ans = dfs(dfs, n, s);
  cout << ans.val() << endl;
}
// int main() {
//   int n, s, K;
//   cin >> n >> s >> K;

//   vector<mint> dp(s + 5);
//   dp[0] = 1;
//   rep(i, n) {
//     vector<mint> sum(s + 5);
//     rep(j, s + 1) {
//       if (j) sum[j] = sum[j - 1];
//       sum[j] += dp[j];
//     }
//     dp = vector<mint>(s + 5);
//     rep(j, s + 1) {
//       //   for (int k = K; k < s; ++k) {
//       //     if (j + k > s) break;
//       //   }
//       //   if (j + K > s) break;
//       dp[j] = sum[j];
//       if (j - K >= 0) dp[j] -= sum[j - K];
//     }
//   }
//   cout << dp[s].val() << endl;
//   return 0;
// }