#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int MAXN=1e6+7;
int n,m;
ll a[MAXN],b[MAXN];
ll c[MAXN];
int vis[MAXN];
bool check()
{
	int tot=0;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			c[++tot]=(a[i]*b[j])%(n*m);
		}
	}
	sort(c+1,c+1+n*m);
	for(int i=1;i<=n*m;i++)
	{
		if(c[i]!=i-1) return 0; 
	}
	return 1;
}
int num=0;
void dfs2(int x,int pre)
{
	if(x==m+1)
	{
		// if(!check()) return;
		num++;
		for(int i=1;i<=n;i++) cout<<a[i]<<' ';cout<<'\n';
		for(int i=1;i<=m;i++) cout<<b[i]<<" ";cout<<'\n';
		cout<<endl;
		if(num>=100||a[1]!=1) exit(0);
		
		return;
	}
	for(int i=pre+1;i<=n*m;i++)
	{

		bool flag=1;
		for(int j=1;j<=n;j++)
		{
			ll sss=(i*a[j])%(n*m);
			if(vis[sss]) flag=0;
			vis[sss]++;
		}
		if(flag)
		{
			b[x]=i;
			dfs2(x+1,i);
		}
		for(int j=1;j<=n;j++)
		{
			ll sss=(i*a[j])%(n*m);
			vis[sss]--;
		}
	}
	return;
}
void dfs1(int x,int pre)
{
	if(x==n+1)
	{
		dfs2(1,0);
		return;
	}
	for(int i=pre+1;i<=n*m;i++)
	{
		a[x]=i;
		dfs1(x+1,i);
	}
	return;
}
bool work(int x)
{
	for(int i=2;i*i<=x;i++)
	{
		if(x%i==0) return 0;
	}
	return 1;
}
int main()
{

	// ios::sync_with_stdio(false);
	// cin.tie(0);
	// cout.tie(0);
	
	// cin>>n>>m;
	// dfs1(1,0);
	// cout<<"-1";
	int t;
	cin>>t;
	while(t--)
	{
		cin>>n>>m;
		int nn=n,mm=m;
		if(n>m) swap(n,m);
		if(gcd(n,m)!=1)
		{
			cout<<"No\n";
			continue;
		}
		for(int i=1;i<=n;i++) a[i]=1LL+(i-1)*m;
		for(int i=1;i<=m;i++) b[i]=1LL+(i-1)*n;
		// if(check()) cout<<"Yes\n";
		// else cout<<"No\n";
		cout<<"Yes\n";
		if(nn>mm)
		{
			for(int i=1;i<=m;i++) cout<<b[i]<<" ";cout<<'\n';
			for(int i=1;i<=n;i++) cout<<a[i]<<" ";cout<<'\n';
		}
		else
		{
			
			for(int i=1;i<=n;i++) cout<<a[i]<<" ";cout<<'\n';
			for(int i=1;i<=m;i++) cout<<b[i]<<" ";cout<<'\n';
		}
		
	}
    return 0;
}
/*

*/