#include <bits/stdc++.h>
using namespace std;
const int mod = 1000000007;
inline int add(int a, int b) {
  if ((a += b) >= mod) a -= mod;
  return a;
}
inline int mult(int a, int b) {
  long long t = 1ll * a * b;
  if (t >= mod) t %= mod;
  return t;
}
int c[55][55], n, m, f[55][55][2];
inline void add(int &a, int b, int c, int d) {
  a = add(a, mult(b, mult(c, d)));
}
int main() {
  c[0][0] = 1;
  for (int i = 1; i <= 50; i++) {
    c[i][0] = 1;
    for (int j = 1; j <= 50; j++) c[i][j] = add(c[i - 1][j - 1], c[i - 1][j]);
  }
  scanf("%d%d", &n, &m);
  if (n <= 1) {
    puts("0");
    return 0;
  }
  f[1][0][0] = 1;
  f[0][0][1] = 1;
  for (int i = 2; i <= n; i++) {
    for (int j = 1; j <= (i >> 1); j++) {
      for (int k1 = 0, k2 = i - 1; k1 <= k2; k1++, k2--) {
        for (int p = 0; p <= j; p++) {
          int C;
          if (k1 != k2)
            C = c[i - 1][k1];
          else
            C = c[i - 2][k1 - 1];
          if (k1)
            C = mult(C, mult(k1, k2));
          else
            C = mult(C, k2);
          add(f[i][j][0], f[k1][p][1], f[k2][j - p][1], C);
          add(f[i][j][1], f[k1][p][1], f[k2][j - p - 1][0], C);
          add(f[i][j][1], f[k1][p][0], f[k2][j - p - 1][1], C);
          add(f[i][j][1], f[k1][p][0], f[k2][j - p - 1][0], C);
        }
      }
    }
  }
  printf("%d\n", add(f[n][m][0], f[n][m][1]));
  return 0;
}
