#include<bits/stdc++.h>
#define ll long long
#define fi first
#define se second
using namespace std;
const ll mod = 1e9 + 7;
ll dp[2][301][301];
ll hiy[3001][3001];
int main() {
dp[0][0][0] = 1;
for(int i=0; i<=3000; i++) {
for(int j=hiy[i][0]=1; j<=i; j++) {
hiy[i][j] = hiy[i-1][j-1] + hiy[i-1][j];
hiy[i][j]%=mod;
}
}
int H,W;
cin >> H >> W;
for(int i=0; i<H; i++) {
for(int j=0; j<=W; j++) {
for(int k=0; k<=W; k++) {
if(!dp[0][j][k]) continue;
int xp = W-k-j;
if(k) {
dp[1][j+1][k-1]+=dp[0][j][k] * k % mod;
dp[1][j+1][k-1]%=mod;
}
if(xp) {
dp[1][j+1][k]+=dp[0][j][k] * 3 * xp % mod;
dp[1][j+1][k]%=mod;
dp[1][j][k+1]+=dp[0][j][k] * xp%mod;
dp[1][j][k+1]%=mod;
}
if(xp>=2) {
dp[1][j+2][k]+=dp[0][j][k] * hiy[xp][2] % mod;
dp[1][j+2][k]%=mod;
}
dp[1][j][k]+=dp[0][j][k];
dp[1][j][k]%=mod;
}
}
for(int j=0; j<=W; j++) {
for(int k=0; k<=W; k++) {
dp[0][j][k] = dp[1][j][k];
dp[1][j][k] = 0;
}
}
}
ll jwb = 0;
for(int i=0; i<=W; i++) {
for(int j=0; j<=W; j++) {
jwb+=dp[0][i][j];
jwb%=mod;
}
}
jwb+=(mod-1);
if(jwb>=mod) jwb-=mod;
cout << jwb << endl;
return 0;
}
