#include<bits/stdc++.h>
using namespace std;
using ll=long long;
const int N=1e3+9,mod=1e9+7;
int T,n,c,ans,a[N],fac[N],C[(N<<1)][N],A[N];
inline int inv(int x){
	if(x<2)return 1;
	return (ll)(mod-mod/x)*inv(mod%x)%mod;
}
void dfs(int u,int s){
	if(u==c){
		int cur=A[n]*fac[n];
		if(s==0){
			for(int i=0;i<c;i++)cur=(ll)cur*inv((ll)A[a[i]]*fac[a[i]]%mod)%mod;
			ans=(ans+cur)%mod;
			// cerr<<">> "<<cur<<endl;
			// for(int i=0;i<c;i++)cerr<<a[i]<<" \n"[i+1==c];
		}
		return;
	}
	for(int x=0;x<=s;x++){
		a[u]=x;
		// cerr<<"cho: "<<u<<" "<<x<<endl;
		dfs(u+1,s-x);
	}
}
int solve(){
	ans=0;
	for(int i=0;i<=n;i++)A[i]=(ll)C[i<<1][i]*inv(i+1)%mod;
	// for(int i=0;i<=n;i++)cerr<<A[i]<<" \n"[i==n];
	dfs(0,n);
	return ans;
}
int main(){
#ifdef memset0
	freopen("G.in","r",stdin);
#endif
	ios::sync_with_stdio(0),cin.tie(0);
	fac[0]=1;
	for(int i=1;i<N;i++)fac[i]=(ll)fac[i-1]*i%mod;
	for(int i=0;i<(N<<1);i++){
		C[i][0]=1;
		for(int j=1;j<=i&&j<N;j++)C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
	}
	cin>>T;
	while(T--){
		cin>>n>>c;
		cout<<solve()<<endl;
	}
}