# include <bits/stdc++.h>
using namespace std;
using ll=long long;
using vi=vector<int>;
constexpr int mod=1e9+7,iv[]={0,1,500000004,333333336,250000002};
int n,N=0,X[210],Y[210],val[410],f[210][210];
vi F[210][210],M[410];
bool flag[210][210];
template<typename T>
inline void upd(int &x,const T &y){x=(x+y)%mod;}
int power(int a,int b)
{
	int ans=1;
	for(;b;b>>=1,a=(ll)a*a%mod)
		if(b&1) ans=(ll)ans*a%mod;
	return ans;
}
vi operator+(vi A,const vi &B)
{
	for(int i=0;i<=N;i++) upd(A[i],B[i]);
	return A;
}
vi operator-(vi A,const vi &B)
{
	for(int i=0;i<=N;i++) upd(A[i],mod-B[i]);
	return A;
}
vi operator*(vi A,int x)
{
	for(int &i:A) i=(ll)i*x%mod;
	return A;
}
vi &operator+=(vi &A,const vi &B)
{
	for(int i=0;i<=N;i++) upd(A[i],B[i]);
	return A;
}
inline vi getvi(int x,int y){vi v(N+1);v[x]=y;return v;}
inline int getD(int x,int y){return (x>1)+(x<n)+(y>1)+(y<n);}
inline int getP(int x,int y){return iv[getD(x,y)];}
vi calcP(int x,int y)
{
	vi ans=getvi(N,x==X[0] && y==Y[0]);
	if(x>1 && !flag[x-1][y]) ans+=F[x-1][y]*getP(x-1,y);
	if(x<n && !flag[x+1][y]) ans+=F[x+1][y]*getP(x+1,y);
	if(y>1 && !flag[x][y-1]) ans+=F[x][y-1]*getP(x,y-1);
	if(y<n && !flag[x][y+1]) ans+=F[x][y+1]*getP(x,y+1);
	return ans;
}
vi calcS(int x,int y)
{
	vi ans(N+1);
	if(x>1 && !flag[x-1][y]) ans+=(F[x-1][y]+getvi(N,f[x-1][y]))*getP(x-1,y);
	if(x<n && !flag[x+1][y]) ans+=(F[x+1][y]+getvi(N,f[x+1][y]))*getP(x+1,y);
	if(y>1 && !flag[x][y-1]) ans+=(F[x][y-1]+getvi(N,f[x][y-1]))*getP(x,y-1);
	if(y<n && !flag[x][y+1]) ans+=(F[x][y+1]+getvi(N,f[x][y+1]))*getP(x,y+1);
	return ans;
}
void Gauss()
{
	for(int i=0;i<N;i++)
	{
		for(int j=i;j<N;j++)
			if(M[j][i])
			{
				if(i!=j) swap(M[i],M[j]);
				break;
			}
		int invn=power(M[i][i],mod-2);
		for(int j=i;j<=N;j++) M[i][j]=(ll)M[i][j]*invn%mod;
		for(int j=i+1;j<N;j++)
		{
			int t=mod-M[j][i];
			for(int k=i;k<=N;k++) upd(M[j][k],(ll)t*M[i][k]);
		}
	}
	for(int i=0;i<N;i++) val[i]=(mod-M[i][N])%mod;
	for(int i=N-2;i>=0;i--)
		for(int j=i+1;j<N;j++)
			upd(val[i],(ll)(mod-M[i][j])*val[j]);
}
int main()
{
	int T,K;
	cin>>T;
	while(T--)
	{
		int cnt=0,cur=0;
		scanf("%d%d",&n,&K);
		for(int i=1;i<=K;i++) scanf("%d%d",&X[i],&Y[i]),flag[X[i]][Y[i]]=1;
		scanf("%d%d",&X[0],&Y[0]);
		N=n;
		for(int i=2;i<=n;i++)
			for(int j=1;j<=n;j++)
				N+=flag[i][j];
		for(int i=1;i<=n;i++) fill(F[i]+1,F[i]+n+1,vi(N+1));
		for(int i=1;i<=n;i++) F[1][i]=getvi(cnt++,1);
		for(int i=1;i<n;i++)
			for(int j=1;j<=n;j++)
				if(flag[i+1][j]) M[cur++]=F[i][j]-calcP(i,j),F[i+1][j]=getvi(cnt++,1);
				else F[i+1][j]=(F[i][j]-calcP(i,j))*getD(i+1,j);
		for(int i=1;i<=n;i++) M[cur++]=F[n][i]-calcP(n,i);
		Gauss();
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
			{
				f[i][j]=F[i][j][N];
				for(int k=0;k<N;k++) upd(f[i][j],(ll)F[i][j][k]*val[k]);
			}
		cnt=cur=0;
		for(int i=1;i<=n;i++) fill(F[i]+1,F[i]+n+1,vi(N+1));
		for(int i=1;i<=n;i++) F[1][i]=getvi(cnt++,1);
		for(int i=1;i<n;i++)
			for(int j=1;j<=n;j++)
				if(flag[i+1][j]) M[cur++]=F[i][j]-calcS(i,j),F[i+1][j]=getvi(cnt++,1);
				else F[i+1][j]=(F[i][j]-calcS(i,j))*getD(i+1,j);
		for(int i=1;i<=n;i++) M[cur++]=F[n][i]-calcS(n,i);
		Gauss();
		for(int i=1;i<=K;i++)
		{
			int ans;
			if(f[X[i]][Y[i]])
			{
				ans=F[X[i]][Y[i]][N];
				for(int j=0;j<N;j++) upd(ans,(ll)F[X[i]][Y[i]][j]*val[j]);
				ans=(ll)ans*power(f[X[i]][Y[i]],mod-2)%mod;
			}
			else ans=-1;
			printf("%d%c",ans," \n"[i==K]);
		}
		for(int i=1;i<=n;i++) memset(flag[i]+1,0,sizeof(bool)*n);
		for(int i=1;i<=n;i++) memset(f[i]+1,0,sizeof(int)*n);
	}
	return 0;
}