#include <cstdio>
using namespace std;
template<typename T>
T read(){
	char c=getchar();T x=0;
	while(c<48||c>57) c=getchar();
	do x=(x<<1)+(x<<3)+(c^48),c=getchar();
	while(c>=48&&c<=57);
	return x;
}
typedef long long ll;
const int P=1000000007;
int qp(int a,int b=P-2){
	int res=1;
	while(b){
		if(b&1) res=(ll)res*a%P;
		a=(ll)a*a%P;b>>=1;
	}
	return res;
}
int fac[1000003],fiv[1000003];
int comb(int a,int b){
	if(b<0||b>a) return 0;
	return (ll)fac[a]*fiv[b]%P*fiv[a-b]%P;
}
void inc(int &x,int v){if((x+=v)>=P) x-=P;}
void dec(int &x,int v){if((x-=v)<0) x+=P;}
int main(){
	int tc=read<int>();
	fac[0]=1;
	for(int i=1;i<=1000000;++i) fac[i]=(ll)fac[i-1]*i%P;
	fiv[1000000]=qp(fac[1000000]);
	for(int i=1000000;i;--i) fiv[i-1]=(ll)fiv[i]*i%P;
	while(tc--){
		int n=read<int>();
		ll m=read<ll>();
		if(m>n){puts("0");continue;}
		if(m==n){printf("%d\n",fac[n-1]);continue;}
		m=n-m;
		int pw=1,res=0;
		for(int i=0;i<=m;++i){
			if((m^i)&1) dec(res,(ll)comb(m,i)*pw%P*comb(n-m+i-1,m-1)%P);
			else inc(res,(ll)comb(m,i)*pw%P*comb(n-m+i-1,m-1)%P);
			inc(pw,pw);
		}
		res=(ll)res*fiv[m]%P*fac[n]%P*qp((P+1)>>1,m)%P;
		printf("%d\n",res);
	}
	return 0;
}