#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int N = 300, M = 1e9 + 7;
int dp[N][4][N];
int f[4][5] = {1, 1, 0, 0, 0, 1, 3, 1, 0, 0, 1, 5, 5, 1, 0, 1, 8, 16, 8, 1};
int x[4] = {3, 4, 5, 4}, y[4] = {0, 1, 2, 4};
inline void update(int &x, ll y) { x = (x + y) % M; }
int g[N * 2];
int main() {
  int n, k;
  scanf("%d %d", &n, &k);
  ll prod = 1;
  for (int i = 1; i <= k; i++) prod = prod * i % M;
  for (int i = 0; i < 4; i++)
    for (int j = 0; j <= i + 1; j++) dp[0][i][j] = f[i][j];
  for (int i = 1; i <= n; i++) {
    for (int j = 0; j < 4; j++) {
      int x = ::x[j], y = ::y[j];
      for (int t = 0; t <= k; t++) dp[i][j][t] = dp[i - 1][j][t];
      for (int u = 1; u <= x; u++)
        for (int t = k; t; t--) update(dp[i][j][t], dp[i][j][t - 1]);
      for (int u = 1; u <= y; u++) {
        for (int t = k; t >= 2; t--)
          update(dp[i][j][t], dp[i][j][t - 1] * 3ll + dp[i][j][t - 2]);
        update(dp[i][j][1], dp[i][j][0] * 3ll);
      }
    }
    for (int x = 0; x <= k; x++) g[x] = 0;
    for (int x = 0; x <= k; x++)
      for (int y = 0; y <= k; y++)
        update(g[x + y], (ll)dp[i][0][x] * dp[i][0][y]);
    for (int t = k; t; t--) {
      update(dp[i][3][t], dp[i][2][t - 1] * 4ll);
      update(dp[i][3][t], 4ll * dp[i][1][t - 1]);
      if (t >= 2) {
        update(dp[i][3][t], g[t - 2] * 2ll);
        update(dp[i][3][t], 4ll * dp[i][1][t - 2]);
        update(dp[i][3][t], 8ll * dp[i][0][t - 2]);
        if (t >= 3) update(dp[i][3][t], 4ll * dp[i][0][t - 3]);
      }
    }
    update(dp[i][3][2], 2);
    update(dp[i][3][3], 4);
    update(dp[i][3][4], 1);
    for (int t = k; t; t--) {
      update(dp[i][2][t], dp[i][1][t - 1] * 2ll);
      update(dp[i][2][t], g[t - 1]);
      update(dp[i][2][t], dp[i][0][t - 1] * 2ll);
      if (t >= 2) update(dp[i][2][t], dp[i][0][t - 2] * 3ll);
    }
    update(dp[i][2][2], 2);
    update(dp[i][2][3], 1);
    for (int t = k; t; t--) update(dp[i][1][t], dp[i][0][t - 1] * 2ll);
    update(dp[i][1][1], 1);
    update(dp[i][1][2], 1);
    update(dp[i][0][1], 1);
  }
  printf("%lld\n", prod * dp[n][3][k] % M);
}
