#include using namespace std; #define rep(i,n) REP(i,0,n) #define REP(i,s,e) for(int i=(s); i<(int)(e); i++) #define repr(i, n) REPR(i, n, 0) #define REPR(i, s, e) for(int i=(int)(s-1); i>=(int)(e); i--) #define pb push_back #define all(r) r.begin(),r.end() #define rall(r) r.rbegin(),r.rend() #define fi first #define se second typedef long long ll; typedef vector vi; typedef vector vl; typedef pair pii; typedef pair pll; const ll INF = 1e18; const ll MOD = 1e9 + 7; const double EPS = 1e-8; template T chmax(T& a, const T& b){return a = (a > b ? a : b);} template T chmin(T& a, const T& b){return a = (a < b ? a : b);} template T pow_fast(T x, U n) { T ret = (T)1; while(n) { if(n & (U)1) ret *= x; x *= x; n >>= 1; } return ret; } template struct ModInt { ModType val; ModInt() : val(0) {} ModInt(ModType x) : val(x % MOD) { if(val < 0) val += MOD; }; operator ModType() const { return val; } ModInt &operator+=(const ModInt &m) { val += m.val; if(val >= MOD) val -= MOD; return *this; } ModInt &operator-=(const ModInt &m) { val -= m.val; if(val < 0) val += MOD; return *this; } ModInt &operator*=(const ModInt &m) { (val *= m.val) %= MOD; return *this; } ModInt &operator/=(const ModInt &m) { *this *= m.inverse(); return *this; } ModInt operator+(const ModInt &m) const { return ModInt(*this) += m; } ModInt operator-(const ModInt &m) const { return ModInt(*this) -= m; } ModInt operator*(const ModInt &m) const { return ModInt(*this) *= m; } ModInt operator/(const ModInt &m) const { return ModInt(*this) /= m; } ModInt inverse() const { return pow_fast(*this, MOD - 2); } ModInt pow(ModType n) const {return pow_fast(*this, n); } friend std::ostream &operator<<(std::ostream &os, const ModInt &rhs) { os << rhs.val; return os; } friend std::istream &operator>>(std::istream &is, ModInt &rhs) { ModType value; is >> value; rhs = ModInt(value); return is; } }; int main(){ using mint = ModInt; int n, d, k; cin >> n >> d >> k; vector cur(n*d+1); cur[0] = 1; rep(i, n) { vector nxt(cur.size()); ll sum = cur[0]; REP(j, 1, d) { nxt[j] = sum; sum += cur[j]; } REP(j, d, cur.size()) { nxt[j] = sum; sum += cur[j]; sum -= cur[j-d]; } swap(cur, nxt); } cout << cur[k] << endl; return 0; }