#include <bits/stdc++.h>
using namespace std;

const int mod = 1000000007;
int dp[2][101][10001];

int main() {
    cin.tie(0)->sync_with_stdio(0);
    int n, m;
    cin >> n >> m;
    dp[0][0][0] = 1;
    for (int i = 1; i <= n; i++) {
        for (int j = 0; j <= i; j++) {
            for (int k = 0; k <= m; k++) {
                if (n-2*i+j+2 >= 2) {
                    dp[i&1][j][k] += (long long)dp[(i&1)^1][j][k] * (n-2*i+j+2) * (n-2*i+j+1) % mod;
                    if (dp[i&1][j][k] >= mod) {
                        dp[i&1][j][k] -= mod;
                    }
                }
                if (j >= 1 && k >= i && (n-2*i+j+1) >= 0) {
                    dp[i&1][j][k] += (long long)dp[(i&1)^1][j-1][k-i] * (n-2*i+j+1) * (2*i-2*j+1) % mod;
                    if (dp[i&1][j][k] >= mod) {
                        dp[i&1][j][k] -= mod;
                    }
                }
                if (j >= 2 && k >= 2*i && i-j+1 >= 0) {
                    dp[i&1][j][k] += (long long)dp[(i&1)^1][j-2][k-2*i] * (i-j+1) * (i-j+1) % mod;
                    if (dp[i&1][j][k] >= mod) {
                        dp[i&1][j][k] -= mod;
                    }
                }
            }
        }
    }
    cout << dp[n&1][n][m] << '\n';
    return 0;
}