#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const int mod = 1e9 + 7;
int n, m;
ll dp[3005][3005];
ll coef(ll x) {
x = x * (x - 1) / 2;
return x % mod;
}
ll f(int r, int c) {
if(r == n + 1) return 1;
if(dp[r][c] != -1) return dp[r][c];
ll ans = f(r+1, c);
//random nebitno
if(c >= 1) ans = (ans + f(r+1, c-1) * 4ll * c) % mod;
//EW
if(c >= 2) ans = (ans + f(r+1, c-2) * coef(c)) % mod;
//NS
if(c >= 1) ans = (ans + f(r+2, c-1) * 1ll * c * (n - r)) % mod;
return dp[r][c] = ans;
}
signed main() {
memset(dp, -1, sizeof(dp));
cin >> n >> m;
cout << (f(1, m) - 1 + mod) % mod << '\n';
}
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |
# | Verdict | Execution time | Memory | Grader output |
---|
Fetching results... |