#include <bits/stdc++.h>
using namespace std;
const int p[5] = {8, 6, 4, 2, 4}, h[5] = {4, 2, 1, 0, 0};
int i, j, k, n, m;
int f[210][210][5][3], v[210][210][5][3];
long long C[210][210];
inline void A(long long &x, long long y) { x = (x + y) % 1000000007; }
inline long long M(long long x, long long y) { return x * y % 1000000007; }
long long cal(int n, int k, int cur, int type) {
  if (!k) return 1;
  if (cur == 5) return 0;
  if (v[n][k][cur][type]) return f[n][k][cur][type];
  v[n][k][cur][type] = 1;
  int i, j;
  long long an = 0;
  if (type == 0) {
    if (!n) return 0;
    int x = cur == 4 ? 2 : 4 - cur;
    for (i = 0; i <= min(x, k); i++)
      A(an, M(cal(n - 1, k - i, cur, 2), C[x][i]));
  } else if (type == 1) {
    if (!n) return 0;
    int x = p[cur], y = h[cur];
    for (i = 0; i <= x; i++)
      for (j = 0; j <= y; j++)
        if (i + 2 * j <= x && i + j <= k)
          A(an, M(cal(n, k - i - j, cur, 0),
                  C[y][j] * C[x - 2 * j][i] % 1000000007));
  } else {
    an = cal(n, k, cur, 1);
    if (cur == 0) {
      if (k > 0) A(an, cal(n, k - 1, 1, 1) * 4 + cal(n, k - 1, 2, 1) * 4);
      if (k > 1)
        A(an, cal(n, k - 2, 2, 1) * 4 + cal(n, k - 2, 4, 1) * 2 +
                  cal(n, k - 2, 3, 1) * 8 + cal(n, k - 2, 5, 1) * 2);
      if (k > 2) A(an, cal(n, k - 3, 3, 1) * 4 + cal(n, k - 3, 5, 1) * 4);
      if (k > 3) A(an, cal(n, k - 4, 5, 1));
    } else if (cur == 1) {
      if (k > 0)
        A(an, cal(n, k - 1, 2, 1) * 2 + cal(n, k - 1, 3, 1) * 2 +
                  cal(n, k - 1, 4, 1));
      if (k > 1) A(an, cal(n, k - 2, 3, 1) * 3 + cal(n, k - 2, 5, 1) * 2);
      if (k > 2) A(an, cal(n, k - 3, 5, 1));
    } else if (cur == 2) {
      if (k > 0) A(an, cal(n, k - 1, 3, 1) * 2 + cal(n, k - 1, 5, 1));
      if (k > 1) A(an, cal(n, k - 2, 5, 1));
    } else if (cur == 4) {
      if (k > 0) A(an, cal(n, k - 1, 3, 1) * 2);
      if (k > 1) A(an, cal(n, k - 2, 5, 1));
    } else if (k > 0)
      A(an, cal(n, k - 1, 5, 1));
  }
  f[n][k][cur][type] = an;
  return an;
}
int main() {
  scanf("%d%d", &n, &k);
  for (i = 0; i <= 200; i++) {
    C[i][0] = 1;
    for (j = 1; j <= i; j++)
      C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % 1000000007;
  }
  long long an = cal(n, k, 0, 2);
  for (i = 1; i <= k; i++) an = an * i % 1000000007;
  printf("%I64d\n", an);
  return 0;
}

### Prompt

Please provide a CPP coded solution to the problem described below:

Vasily the bear has two favorite integers n and k and a pencil. Besides, he's got k jars with different water color paints. All jars are numbered in some manner from 1 to k, inclusive. The jar number i contains the paint of the i-th color. 

Initially the bear took a pencil and drew four segments on the coordinate plane. All of them end at point (0, 0). They begin at: (0, 2n), (0, - 2n), (2n, 0), ( - 2n, 0). Then for each i = 1, 2, ..., n, the bear drew two squares. The first square has the following vertex coordinates: (2i, 0), ( - 2i, 0), (0, - 2i), (0, 2i). The second square has the following vertex coordinates: ( - 2i - 1, - 2i - 1), ( - 2i - 1, 2i - 1), (2i - 1, - 2i - 1), (2i - 1, 2i - 1). After that, the bear drew another square: (1, 0), ( - 1, 0), (0, - 1), (0, 1). All points mentioned above form the set of points A.

<image>

The sample of the final picture at n = 0

<image>

The sample of the final picture at n = 2

The bear decided to paint the resulting picture in k moves. The i-th move consists of the following stages: 

  1. The bear chooses 3 distinct points in set А so that any pair of the chosen points has a segment on the picture between them. The chosen points and segments mark the area that mustn't contain any previously painted points. 
  2. The bear paints the area bounded by the chosen points and segments the i-th color. 



Note that after the k-th move some parts of the picture can stay unpainted.

The bear asked you to calculate, how many distinct ways there are to paint his picture. A way to paint the picture is a sequence of three-element sets of points he chose on each step. Two sequences are considered distinct if there is such number i (1 ≤ i ≤ k), that the i-th members of these sequences do not coincide as sets. As the sought number can be rather large, you only need to calculate the remainder after dividing it by number 1000000007 (109 + 7).

Input

The first line contains two integers n and k, separated by a space (0 ≤ n, k ≤ 200).

Output

Print exactly one integer — the answer to the problem modulo 1000000007 (109 + 7).

Examples

Input

0 0


Output

1


Input

0 1


Output

8


Input

0 2


Output

32


Input

1 1


Output

32

### Response

```cpp
#include <bits/stdc++.h>
using namespace std;
const int p[5] = {8, 6, 4, 2, 4}, h[5] = {4, 2, 1, 0, 0};
int i, j, k, n, m;
int f[210][210][5][3], v[210][210][5][3];
long long C[210][210];
inline void A(long long &x, long long y) { x = (x + y) % 1000000007; }
inline long long M(long long x, long long y) { return x * y % 1000000007; }
long long cal(int n, int k, int cur, int type) {
  if (!k) return 1;
  if (cur == 5) return 0;
  if (v[n][k][cur][type]) return f[n][k][cur][type];
  v[n][k][cur][type] = 1;
  int i, j;
  long long an = 0;
  if (type == 0) {
    if (!n) return 0;
    int x = cur == 4 ? 2 : 4 - cur;
    for (i = 0; i <= min(x, k); i++)
      A(an, M(cal(n - 1, k - i, cur, 2), C[x][i]));
  } else if (type == 1) {
    if (!n) return 0;
    int x = p[cur], y = h[cur];
    for (i = 0; i <= x; i++)
      for (j = 0; j <= y; j++)
        if (i + 2 * j <= x && i + j <= k)
          A(an, M(cal(n, k - i - j, cur, 0),
                  C[y][j] * C[x - 2 * j][i] % 1000000007));
  } else {
    an = cal(n, k, cur, 1);
    if (cur == 0) {
      if (k > 0) A(an, cal(n, k - 1, 1, 1) * 4 + cal(n, k - 1, 2, 1) * 4);
      if (k > 1)
        A(an, cal(n, k - 2, 2, 1) * 4 + cal(n, k - 2, 4, 1) * 2 +
                  cal(n, k - 2, 3, 1) * 8 + cal(n, k - 2, 5, 1) * 2);
      if (k > 2) A(an, cal(n, k - 3, 3, 1) * 4 + cal(n, k - 3, 5, 1) * 4);
      if (k > 3) A(an, cal(n, k - 4, 5, 1));
    } else if (cur == 1) {
      if (k > 0)
        A(an, cal(n, k - 1, 2, 1) * 2 + cal(n, k - 1, 3, 1) * 2 +
                  cal(n, k - 1, 4, 1));
      if (k > 1) A(an, cal(n, k - 2, 3, 1) * 3 + cal(n, k - 2, 5, 1) * 2);
      if (k > 2) A(an, cal(n, k - 3, 5, 1));
    } else if (cur == 2) {
      if (k > 0) A(an, cal(n, k - 1, 3, 1) * 2 + cal(n, k - 1, 5, 1));
      if (k > 1) A(an, cal(n, k - 2, 5, 1));
    } else if (cur == 4) {
      if (k > 0) A(an, cal(n, k - 1, 3, 1) * 2);
      if (k > 1) A(an, cal(n, k - 2, 5, 1));
    } else if (k > 0)
      A(an, cal(n, k - 1, 5, 1));
  }
  f[n][k][cur][type] = an;
  return an;
}
int main() {
  scanf("%d%d", &n, &k);
  for (i = 0; i <= 200; i++) {
    C[i][0] = 1;
    for (j = 1; j <= i; j++)
      C[i][j] = (C[i - 1][j - 1] + C[i - 1][j]) % 1000000007;
  }
  long long an = cal(n, k, 0, 2);
  for (i = 1; i <= k; i++) an = an * i % 1000000007;
  printf("%I64d\n", an);
  return 0;
}

```