#include <bits/stdc++.h>
using namespace std;

// vector<int> ar,br;
// int n,m;
// void dfs2(int num,int p)
// {
//     if(p==m){
//         set<int> kk;
//         for(int i=0;i<n;++i){
//             for(int j=0;j<m;++j){
//                 kk.insert(ar[i]*br[j]%(n*m));
//             }
//         }
//         if(kk.size()==n*m){
//             for(int i=0;i<n;++i)printf("%d ",ar[i]);
//             puts("");
//             for(int i=0;i<m;++i)printf("%d ",br[i]);
//             puts("");
            
//         }
//         return;
//     }
//     for(int i=num;i<n*m;++i){
//         br.push_back(i);
//         dfs2(i,p+1);
//         br.pop_back();
//     }
// }
// void dfs1(int num,int p){
//     if(p==n){
//         dfs2(0,0);
//         return;
//     }
//     for(int i=num;i<n*m;++i){
//         ar.push_back(i);
//         dfs1(i,p+1);
//         ar.pop_back();
//     }
// }
void solve()
{
    int n,m;
    cin>>n>>m;
    vector<int> ar,br;
    int p=n*m;
    for(int i=1,j=1;j<=n;i+=n-1,++j){
        ar.push_back(i%p);
    }
    for(int i=1,j=1;j<=m;i+=m-1,++j){
        br.push_back(i%p);
    }
    vector<int> table(n*m);
    int ok=1;
    for(int i=0;i<n;++i){
        if(!ok)
        {
            break;
        }
        for(int j=0;j<m;++j){
            if(++table[ar[i]*br[j]%(n*m)]>1){
                ok=0;
                break;
            }
            
        }
    }
    if(!ok)puts("No");
    else {
        puts("Yes");
        for(auto v:ar){
            printf("%d ",v);
        }
        puts("");
        for(auto v:br){
            printf("%d ",v);
        }
        puts("");
    }
}

int main() {
    int t=1;
    cin>>t;
    while(t--)solve();
}