#include <bits/stdc++.h>
#define int int64_t
using namespace std;
void solve(){
    int N, na, nb, nc;
    cin>>N>>na>>nb>>nc;
    array<int, 3> q={na, nb, nc};
    sort(q.begin(), q.end(), greater<int> ());
    vector<int> v(N);
    for(auto &i:v)cin>>i;
    sort(v.begin(),v.end());
    vector<int> ans(N);
    int tot = accumulate(v.begin(), v.end(), 0LL);
    int s1=0;
    for(int i=0;i<q[0];i++){
        ans[i]=1;
        s1+=v[i];
    }
    for(int i=q[0]-1;i>=0;i--){
        int z=s1+v[N-q[0]+i]-v[i];
        if(z+z<tot){
            ans[i]=0;
            ans[N-q[0]+i]=1;
        }else{
            break;
        }
    }
    vector<int> aa[3];
    vector<int> p;
    for(int i=0;i<N;i++){
        if(ans[i]==0)p.push_back(v[i]);
        else{
            aa[0].push_back(v[i]);
        }
    }

    // cout<<"p: ";
    // for(auto x:p)cout<<x<<" ";
    // cout<<endl;

    vector<int> ans1(N, 0);
    int tot1=tot-s1;
    int s2=0;
    for(int i=0;i<q[1];i++){
        ans1[i]=2;
        s2+=p[i];
    }
    for(int i=q[1]-1;i>=0;i--){
        int z=s2+p[p.size()-q[1]+i]-p[i];
        if(z+z<tot1){
            ans1[i]=0;
            ans1[p.size()-q[1]+i]=2;
        }else{
            break;
        }
    }

    for(int i=0;i<p.size();i++){
            if(ans1[i])aa[1].push_back(p[i]);
            else aa[2].push_back(p[i]);
    }
    for(int i=0;i<3;i++){
        // for(auto a:aa[i])cout<<a<<endl;
        // cout<<endl;
        // cout<<endl;
        int x = accumulate(aa[i].begin(),aa[i].end(), 0LL);
        if(x+x>=tot){
            cout<<"NO\n";
            return ;
        }
    }
    cout<<"YES\n";
    vector<bool> vis(3);
    for(int i=0;i<3;i++){
        if(aa[i].size()==na && vis[i]==0){
            vis[i]=1;
            for(auto x:aa[i]){
                cout<<x<<" ";
            }
            cout<<"\n";
            break;
        }
    }
    for(int i=0;i<3;i++){
        if(aa[i].size()==nb && vis[i]==0){
            vis[i]=1;
            for(auto x:aa[i]){
                cout<<x<<" ";
            }
            cout<<"\n";
            break;
        }
    }
    for(int i=0;i<3;i++){
        if(aa[i].size()==nc && vis[i]==0){
            vis[i]=1;
            for(auto x:aa[i]){
                cout<<x<<" ";
            }
            cout<<"\n";
            break;
        }
    }
}
signed main(){
    ios_base::sync_with_stdio();
    cin.tie(0);
    int T;
    cin>>T;
    while(T--){
        solve();
    }
}