#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define endl '\n'
#define debug(x) cout<<#x<<": "<<x<<endl
const ll mod=1e9+7;
const ll N=1e5+10;

ll qpow(ll a,ll b)
{
	ll ans=1;
	a%=mod;
	while(b)
	{
		if(b&1) ans=ans*a%mod;
		a=a*a%mod;
		b>>=1;
	}
	return ans;
}

ll f[N],g[N],inv[N];
ll t[N];
void init_inverse1()//求1开始的连续逆元
{
	inv[1]=1;
	for(int i=2;i<N;i++)
		inv[i]=(ll)(mod-mod/i)*inv[mod%i]%mod;
}

void init_combinatorial()
{
	init_inverse1();//用init_inverse2()会稍微快一点
	f[0]=g[0]=1;
	for(int i=1;i<N;i++)
		f[i]=(ll)f[i-1]*i%mod;
	for(int i=1;i<N;i++)
		g[i]=(ll)g[i-1]*inv[i]%mod;
}

ll getC(int n,int m)//n大m小
{
	if(n<m) return 0;
	return ((ll)f[n]*g[m]%mod)*g[n-m]%mod;
}

void solve()
{
	ll n,m;
	cin>>n>>m;
	if(m>n)
	{
		cout<<0<<endl;
		return;
	}
	else if(n==m)
	{
		cout<<f[n-1]*qpow(2,mod-2)%mod<<endl;
		return;
	}
	ll k=n-m;
	ll ans=qpow(qpow(2,mod-2),k);
	ll now=0;
	for(int i=0;i<=n-k;i++)
	{
		now+=(t[k-i]*getC(k,i)%mod)*(getC(n-i-1,n-k-i)%mod)*qpow(-1,i%2);
		now%=mod;
	}
	ans=ans*now%mod;
	ans=ans*f[n]%mod;
	ans=ans*g[k]%mod;
	cout<<ans<<endl;
}

int main()
{
	ios::sync_with_stdio(0);
	cin.tie(0),cout.tie(0);
	
	int T=1;
	cin>>T;
	init_combinatorial();
	t[0]=1;
	for(int i=1;i<N;i++) t[i]=t[i-1]*2%mod;
	while(T--)
	{
		solve();
	}

	return 0;
}