#include <bits/stdc++.h> using namespace std; const int mod = 1000000007; int mem[3003][3003]; long long dp(int H, int W) { if(H < 0) return 0; if(W < 0) return 0; if(H == 0 || W == 0) return 1; if(mem[H][W] != -1) return mem[H][W]; long long ans = 0; ans += dp(H - 1, W); ans += dp(H - 1, W - 2) * ((W * (W - 1)) / 2); ans += 4LL * dp(H - 1, W - 1) * W; ans += dp(H - 2, W - 1) * W * (H - 1); ans %= mod; return mem[H][W] = ans; } int main() { memset(mem, -1, sizeof mem); int H, W; cin >> H >> W; long long ans = dp(H, W) - 1; if(ans < 0) ans += mod; cout << ans << endl; return 0; }