#include <bits/stdc++.h>
using namespace std;
#define x first
#define y second
#define mp(Tx, Ty) make_pair(Tx, Ty)
#define For(Ti, Ta, Tb) for(auto Ti = (Ta); Ti <= (Tb); Ti++)
#define Dec(Ti, Ta, Tb) for(auto Ti = (Ta); Ti >= (Tb); Ti--)
#define debug(...) fprintf(stderr, __VA_ARGS__)
#define range(Tx) begin(Tx),end(Tx) 
const int N = 2005, mod = 1e9 + 7;
long long f[N][N][2];
int main() {
	//freopen("count.in", "r", stdin);
	cin.tie(nullptr)->sync_with_stdio(false);
	int t;
	cin >> t;
	int n, k;
	f[1][0][0] = 1;
	f[1][0][1] = 0;
	For(i, 2, 2001 - 1) {
		For(j, 0, i - 1) {
		 	f[i][j][0] = (f[i - 1][j][0] * max((i - 2 - j), 0) % mod + f[i - 1][j][1] * max(0, (i - 2 - (j - 1))) % mod + f[i - 1][j + 1][0] * (j + 1) % mod + f[i - 1][j + 1][1] * j % mod) % mod;
			f[i][j][1] = (f[i - 1][j - 1][0] * 2 % mod + f[i - 1][j - 1][1] + f[i - 1][j][1]) % mod;
		}
	}
//	For(i, 1, 5) {
//		For(j, 0, 5) {
//			cout << f[i][j][0] + f[i][j][1] << ' ';
//		}
//		cout << endl;
//	}
	while (t--) {
		cin >> n >> k;
		cout << (f[n][k][0] + f[n][k][1]) % mod << '\n';
	}
	return 0;
}