import sys
input=lambda: sys.stdin.readline().rstrip()
n,d,k=map(int,input().split())
mod=10**9+7
A=[0]*(n*d+1)
A[0]=1
for _ in range(n):
  B=[0]*(n*d+1)
  for i in range(n*d+1):
    if not A[i]:
      continue
    else:
      B[i+1]+=A[i]
      if i+d+1<=n*d:
        B[i+d+1]-=A[i]
  A[0]=0
  for i in range(1,n*d+1):
    A[i]=(A[i-1]+B[i])%mod
print(A[k])