#include<bits/stdc++.h>
using namespace std;
const int N = 2e5+10,mod=1e9+7;
typedef long long LL;
int fact[N],infact[N];
int qmi(int a,int b)
{
    int res=1;
    while(b)
    {
        if(b&1)res=(LL)res*a%mod;
        a=(LL)a*a%mod;
        b>>=1;
    }
    return res;
}
void init()
{
    infact[0]=fact[0]=1;
    for(int i=1;i<N;i++)    fact[i]=(LL)fact[i-1]*i%mod;
    infact[N-1]=qmi(fact[N-1],mod-2);
    for(int i= N-2;i>0;i--) infact[i]=(LL)infact[i+1]*(i+1)%mod;
}
int C(int a,int b)
{
    if(b>a) return 0;
    return (LL)fact[a]*infact[b]%mod*infact[a-b]%mod;
}
void solve()
{
    int n,m;
    cin>>n>>m;
    if(m>n)
    {
        puts("0");
        return;
    }
    if(m==n)
    {
        if(n<=2)
        {
            puts("0");
            return;
        }
        cout<<(LL)fact[n-1]*infact[2]%mod<<endl;
        return ;
    }
    if(m==0)
    {
        puts("1");
        return ;
    }
    int t=1,ans=0;
    for(int i=1;i<=n-m;i++)
    {
        t=(LL)t*(2*i-1)%mod;
        ans=(ans+(LL)C(n,n-m-i)*C(m+i,2*i)%mod*t%mod
        *fact[m-i]%mod*C(m-1,i-1)%mod)%mod;
    }
    cout<<ans<<endl;
}
int main()
{
    int T =1;
    cin>>T;
    init();
    while(T--)
    {
        solve();
    }
    return 0;
}