#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e5+100,mod=1e9+7,maxn=1e5;

int fl[N],nf[N],Mi[N];

int mi(int x,int t){
    int d=1;
    while(t){
        if(t%2) d=(ll)d*x%mod;
        x=(ll)x*x%mod;t/=2;
    }
    return d;
}
int ni(int x) {return mi(x,mod-2);}

int C(int n,int m){
    if(m<0||m>n) return 0;
    return (ll)fl[n]*nf[m]%mod*nf[n-m]%mod;
}

int main()
{

    fl[0]=1;for(int i=1;i<=maxn;i++) fl[i]=(ll)fl[i-1]*i%mod;
    nf[maxn]=ni(fl[maxn]);for(int i=maxn-1;i>=0;i--) nf[i]=(ll)nf[i+1]*(i+1)%mod;

    int Ntwo=ni(2);
    Mi[0]=1;for(int i=1;i<=maxn;i++) Mi[i]=(ll)Mi[i-1]*Ntwo%mod;

    int T;scanf("%d",&T);

    while(T--){
        int n;ll m;
        scanf("%d%lld",&n,&m);

        if(n==m){
            printf("%d\n",fl[n-1]);
            continue;
        }

        if(m>n){
            printf("0\n");
            continue;
        }

        if(m==0){
            printf("1\n");
            continue;
        }

        int Ans=0;
        for(int i=1;i<=n/2&&i<=n-m&&i<=m;i++){
            int t=(ll)C(n,m+i)*C(m+i,2*i)%mod*fl[2*i]%mod*Mi[i]%mod*fl[m-i]%mod*C(m-1,i-1)%mod*nf[i]%mod;
            Ans=(Ans+t)%mod;
        }
        printf("%d\n",Ans);

    }

    return 0;
}