#include<stdio.h>
#include<cstring>
#include<algorithm>
#define inf 1000000000
using namespace std;
int i,j,k,l,s,n,m,E[405],e[405],Val[100005],gt,v[405],q[205],Q[205],qq[205],r,R;
int W[405],tot,D[405],P[405],next[100005],last[205],to[100005],val[100005],S,T,dis[205];
int ans,test,x;
inline void doit() {
    for (int i=2;i<=400;i++) if (!v[i]) for (int j=2;j*i<=400;j++) v[j*i]=1;
}
inline bool spfa() {
    int l=0,r=1; for (int i=S;i<=T;i++) E[i]=dis[i]=0;
    dis[S]=1; e[S]=1; qq[1]=S;
    while (l<r) {
        int k=qq[++l];
        for (int i=last[k];i;i=next[i]) {
            if (!val[i]) continue;
            if (dis[to[i]]<dis[k]+Val[i]) {
                E[to[i]]=i^1; dis[to[i]]=dis[k]+Val[i];
                if (!e[to[i]]) e[qq[++r]=to[i]]=1;
            }
        }
        e[k]=0;
    }
    if (dis[T]) return true;
    return false;
}    
inline int Dfs(int o,int p) {
    if (o==T) return p;
    if (!p) return 0;
    int sum=0;
    for (int i=last[o];i;i=next[i]) {
        if (to[E[to[i]]]!=o||!val[i]||dis[to[i]]!=dis[o]+Val[i]) continue;
        int temp=Dfs(to[i],min(p-sum,val[i]));
        if (temp) {
            val[i]-=temp;
            sum+=temp;
            val[i^1]+=temp;
        }
        else dis[to[i]]=0;
    }
    return sum;
}            
inline void add(int o,int p,int w,int e) {
    next[++tot]=last[o]; last[o]=tot; to[tot]=p; val[tot]=w; Val[tot]=e;
}
inline void check() {
    tot=1;
    T=r+s+1;
    for (i=S;i<=T;i++) last[i]=0;
    int gtw=0;
    for (i=1;i<=r;i++) D[i]=val[P[i]];
    for (i=1;i<=r;i++) {
        W[i]=0;
        if (val[P[i]]) {
            W[i]=1;
            if (v[q[i]+1]) {
                printf("NO\n");
                return ;
            }
            if (D[i]==1) add(i+s,T,1,0),add(T,i+s,0,0),gtw++;
            else if (D[i]==2) add(i+s,T,2,0),add(T,i+s,0,0),gtw+=2;
        }
    }
    if (!gtw&&s<=2&&s) {
        printf("NO\n");
        return ;
    }
    for (i=1;i<=s;i++) add(S,i,2,0),add(i,S,0,0);
    for (i=1;i<=s;i++) for (j=1;j<=r;j++) if (W[j]) add(i,j+s,1,0),add(j+s,i,0,0);
    ans=0;
    while (spfa()) ans+=Dfs(S,inf);
    if (ans==gtw) printf("YES\n");
    else printf("NO\n");
}
inline void lianbian() {
    memset(last,0,sizeof(last));
    S=0,T=r+R+1; tot=1;
    for (i=1;i<=R;i++) add(S,i,2,0),add(i,S,0,0);
    for (i=1;i<=r;i++) P[i]=tot+1,add(R+i,T,2,0),add(T,i+R,0,0);
    for (i=1;i<=R;i++) for (j=1;j<=r;j++) if (!v[Q[i]+q[j]]) {
        if (v[q[j]+1]) 
        add(i,j+R,1,1),add(j+R,i,0,-1);
        else add(i,j+R,1,0),add(j+R,i,0,0);
    }
}
int main() {
    scanf("%d",&test);
    doit();
    while (test--) {
        scanf("%d",&n);
        r=R=s=0;
        for (i=1;i<=n;i++) {
            scanf("%d",&x);
            if (x==1) s++;
            else if (!(x%2)) q[++r]=x;
            else Q[++R]=x;
        }
            if (s==1) {
                Q[++R]=1;
                lianbian();
                ans=0;
                while (spfa()) ans+=Dfs(S,inf);
                if (ans==R*2&&ans==r*2) printf("YES\n");
                else printf("NO\n");
            } 
            else {     
            lianbian();
            gt=0;    
            while (spfa()) gt+=Dfs(S,inf);
            if (gt!=R*2) printf("NO\n");
            else 
            check();
        }
    }
}