#include using namespace std; inline int toInt(string s) {int v; istringstream sin(s);sin>>v;return v;} template inline string toString(T x) {ostringstream sout;sout< VI; typedef vector VVI; typedef vector VS; typedef pair P; typedef tuple tpl; #define ALL(a) (a).begin(),(a).end() #define RALL(a) (a).rbegin(), (a).rend() #define SORT(c) sort((c).begin(),(c).end()) #define REVERSE(c) reverse((c).begin(),(c).end()) #define LB(a,x) lower_bound((a).begin(), (a).end(), x) - (a).begin() #define UB(a,x) upper_bound((a).begin(), (a).end(), x) - (a).begin() #define FOR(i,a,b) for(int i=(a);i<(b);++i) #define REP(i,n) FOR(i,0,n) #define RFOR(i,a,b) for(int i=(a)-1;i>=(b);--i) #define RREP(i,n) RFOR(i,n,0) #define en "\n" const double EPS = 1e-9; const double PI = acos(-1.0); const int INT_INF = 2147483647; const long long LL_INF = 1LL<<60; const long long MOD = 1000000007; #define CLR(a) memset((a), 0, sizeof(a)) #define dump(x) cerr << #x << " = " << (x) << endl; #define debug(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ")" << " " << __FILE__ << endl; template inline bool chmax(T& a, T b) { if (a < b) { a = b; return true; } return false; } template inline bool chmin(T& a, T b) { if (a > b) { a = b; return true; } return false; } template struct LazySegmentTree{ typedef function F; typedef function G; typedef function H; typedef function P; int n; F f; G g; H h; P p; T d1; E d0; vector val; vector lazy; LazySegmentTree(vector a, F f_, G g_, H h_, T d1_, E d0_, P p_=[](E a,int b){return a;}): f(f_),g(g_),h(h_),d1(d1_),d0(d0_),p(p_){ int sz = a.size(); n = 1; while (n < sz) n <<= 1; val.clear(); val.resize(2 * n - 1, d1); lazy.clear(); lazy.resize(2 * n - 1, d0); for (int i = 0; i < sz; i++) val[i + n - 1] = a[i]; for (int i = n - 2; i >= 0; i--) val[i] = f(val[i * 2 + 1], val[i * 2 + 2]); } LazySegmentTree() {}; void init(vector a, F f_, G g_, H h_, T d1_, E d0_, P p_=[](E a,int b){return a;}){ f=f_; g=g_; h=h_; d1=d1_; d0=d0_; p=p_; int sz = a.size(); n = 1; while (n < sz) n <<= 1; val.clear(); val.resize(2 * n - 1, d1); lazy.clear(); lazy.resize(2 * n - 1, d0); for (int i = 0; i < sz; i++) val[i + n - 1] = a[i]; for (int i = n - 2; i >= 0; i--) val[i] = f(val[i * 2 + 1], val[i * 2 + 2]); } inline void eval(int k, int len){ if(lazy[k] == d0) return; if(2*k+1 < 2*n-1){ lazy[2*k+1] = h(lazy[2*k+1], lazy[k]); lazy[2*k+2] = h(lazy[2*k+2], lazy[k]); } val[k] = g(val[k], p(lazy[k], len)); lazy[k] = d0; } void update(int a, int b, E x, int k=0, int l=0, int r=-1){ if(r < 0) r = n; eval(k, r-l); if(r <=a || b <= l) return; if(a <= l && r <= b){ lazy[k] = h(lazy[k], x); eval(k, r-l); } else{ update(a, b, x, 2*k+1, l, (l+r)/2); update(a, b, x, 2*k+2, (l+r)/2, r); val[k] = f(val[2*k+1], val[2*k+2]); } } T query(int a, int b, int k=0, int l=0, int r=-1){ if(r < 0) r = n; eval(k, r-l); if(r <= a || b <= l) return d1; if(a <= l && r <= b) return val[k]; T val_l = query(a, b , 2*k+1, l, (l+r)/2); T val_r = query(a, b , 2*k+2, (l+r)/2, r); return f(val_l, val_r); } T get(int pos){ int k = pos + n - 1; eval(k, 1); return val[k]; } }; struct mint { long long x; mint(long long x=0):x((x%MOD+MOD)%MOD){} mint operator-() const { return mint(-x);} mint& operator+=(const mint a) { if ((x += a.x) >= MOD) x -= MOD; return *this; } mint& operator-=(const mint a) { if ((x += MOD-a.x) >= MOD) x -= MOD; return *this; } mint& operator*=(const mint a) { (x *= a.x) %= MOD; return *this; } mint operator+(const mint a) const { mint res(*this); return res+=a; } mint operator-(const mint a) const { mint res(*this); return res-=a; } mint operator*(const mint a) const { mint res(*this); return res*=a; } mint pow(long long t) const { if (!t) return 1; mint a = pow(t>>1); a *= a; if (t&1) a *= *this; return a; } // for prime MOD mint inv() const { return pow(MOD-2); } mint& operator/=(const mint a) { return (*this) *= a.inv(); } mint operator/(const mint a) const { mint res(*this); return res/=a; } }; int main(void){ auto f = [](mint a, mint b){return a+b;}; auto g = [](mint a, ll b){return a+b;}; auto h = [](ll a, ll b){return a+b;}; auto p = [](ll a, int b){return a*(ll)b;}; ll N,d,K; cin >> N >> d >> K; vector a(K+1); REP(i,K+1) a[i] = 0; LazySegmentTree st[N+1]; REP(i,N+1) st[i].init(a,f,g,h,mint(0),0,p); st[0].update(0,1,1); REP(i,N)REP(j,K){ mint t = st[i].query(j,j+1); if(t.x > 0){ st[i+1].update(j+1,min(K,j+d)+1,t.x); } } ll ans = st[N].query(K,K+1).x; cout << ans << en; return 0; }