Problem: Yaroslav calls an array of r integers a1, a2, ..., ar good, if it meets the following conditions: |a1 - a2| = 1, |a2 - a3| = 1, ..., |ar - 1 - ar| = 1, |ar - a1| = 1, at that <image>. 

An array of integers b1, b2, ..., br is called great, if it meets the following conditions:

  1. The elements in it do not decrease (bi ≤ bi + 1). 
  2. If the inequalities 1 ≤ r ≤ n and 1 ≤ bi ≤ m hold. 
  3. If we can rearrange its elements and get at least one and at most k distinct good arrays. 



Yaroslav has three integers n, m, k. He needs to count the number of distinct great arrays. Help Yaroslav! As the answer may be rather large, print the remainder after dividing it by 1000000007 (109 + 7).

Two arrays are considered distinct if there is a position in which they have distinct numbers.

Input

The single line contains three integers n, m, k (1 ≤ n, m, k ≤ 100).

Output

In a single line print the remainder after dividing the answer to the problem by number 1000000007 (109 + 7).

Examples

Input

1 1 1


Output

0


Input

3 3 3


Output

2

Solution (in CPP):

Solution:
#include <bits/stdc++.h>
using namespace std;
const int intmax = 0x3f3f3f3f;
const long long lldmax = 0x3f3f3f3f3f3f3f3fll;
double eps = 1e-8;
template <class T>
inline void checkmin(T &a, T b) {
  if (b < a) a = b;
}
template <class T>
inline void checkmax(T &a, T b) {
  if (b > a) a = b;
}
template <class T>
inline T sqr(T x) {
  return x * x;
}
template <class T>
inline T lowbit(T n) {
  return (n ^ (n - 1)) & n;
}
template <class T>
inline int countbit(T n) {
  return (n == 0) ? 0 : (1 + countbit(n & (n - 1)));
}
template <class T>
inline T checkmod(T n, T m) {
  return (n % m + m) % m;
}
template <class T>
inline T gcd(T a, T b) {
  if (a < 0) return gcd(-a, b);
  if (b < 0) return gcd(a, -b);
  return (b == 0) ? a : gcd(b, a % b);
}
template <class T>
inline T euclid(T a, T b, T &x, T &y) {
  if (a < 0) {
    T d = euclid(-a, b, x, y);
    x = -x;
    return d;
  }
  if (b < 0) {
    T d = euclid(a, -b, x, y);
    y = -y;
    return d;
  }
  if (b == 0) {
    x = 1;
    y = 0;
    return a;
  } else {
    T d = euclid(b, a % b, x, y);
    T t = x;
    x = y;
    y = t - (a / b) * y;
    return d;
  }
}
template <class T>
inline int dblcmp(T a, const T b) {
  a -= b;
  return a > eps ? 1 : (a < -eps ? -1 : 0);
}
template <class T>
inline int sgn(T a) {
  return a > eps ? 1 : (a < -eps ? -1 : 0);
}
const int N = 102;
const long long P = 1000000007;
long long dp[53][N][52][N];
int comb[N][N];
int main() {
  for (int i = 0; i < N; ++i) {
    comb[i][0] = comb[i][i] = 1;
    for (int j = 1; j < i; ++j) {
      comb[i][j] = min(N, comb[i - 1][j] + comb[i - 1][j - 1]);
    }
  }
  int n, m, kind;
  cin >> n >> m >> kind;
  long long ans = -m;
  dp[1][0][1][1] = 1;
  for (int i = 1; i <= n / 2 + 1 && i <= m; ++i)
    for (int j = 0; j <= n; ++j)
      for (int l = 0; l <= n / 2 + 1; ++l)
        for (int k = 1; k <= kind; ++k) {
          dp[i][j][l][k] %= P;
          if (dp[i][j][l][k] == 0) continue;
          ans += dp[i][j][l][k] * (m - i + 1);
          ans %= P;
          for (int nl = 1; 2 * nl + j <= n && comb[l + nl - 1][nl] * k <= kind;
               ++nl) {
            dp[i + 1][j + 2 * nl][nl][comb[l + nl - 1][nl] * k] +=
                dp[i][j][l][k];
            dp[i + 1][j + 2 * nl][nl][comb[l + nl - 1][nl] * k] %= P;
          }
        }
  cout << ans << endl;
  return 0;
}