#include<cstdio>
#include<iostream>
#include<vector>
#include<cstring>
#include<queue>
#include<algorithm>
#include<cassert>
#define RI register int
#define CI const int&
using namespace std;
const int N=505,INF=1e9;
int T,n,m,b[N],valid[N],mat[N],rb[N]; vector <int> a[N];
namespace Network_Flow
{
    const int NN=1005,MM=5e6+5;
    struct edge
    {
        int to,nxt,v;
    }e[MM<<1]; int cnt=1,head[NN],cur[NN],dep[NN],s,t;
    inline void addedge(CI x,CI y,CI z)
    {
        e[++cnt]=(edge){y,head[x],z}; head[x]=cnt;
        e[++cnt]=(edge){x,head[y],0}; head[y]=cnt;
    }
    #define to e[i].to
    inline bool BFS(void)
    {
        memset(dep,0,(t+1)*sizeof(int)); dep[s]=1;
        queue <int> q; q.push(s);
        while (!q.empty())
        {
            int now=q.front(); q.pop();
            for (RI i=head[now];i;i=e[i].nxt)
            if (e[i].v&&!dep[to]) dep[to]=dep[now]+1,q.push(to);
        }
        return dep[t];
    }
    inline int DFS(CI now,CI tar,int dis)
    {
        if (now==tar) return dis; int ret=0;
        for (RI& i=cur[now];i&&dis;i=e[i].nxt)
        if (e[i].v&&dep[to]==dep[now]+1)
        {
            int dist=DFS(to,tar,min(dis,e[i].v));
            if (!dist) dep[to]=INF;
            dis-=dist; ret+=dist; e[i].v-=dist; e[i^1].v+=dist;
            if (!dis) return ret;
        }
        if (!ret) dep[now]=INF; return ret;
    }
    #undef to
    inline int Dinic(int ret=0)
    {
        while (BFS()) memcpy(cur,head,(t+1)*sizeof(int)),ret+=DFS(s,t,INF); return ret;
    }
    inline void build(void)
    {
        cnt=1; for (RI i=0;i<=t;++i) head[i]=0;
        for (RI i=1;i<=n;++i) addedge(s,i,1);
        for (RI i=1;i<=m;++i) if (b[i]>0) addedge(n+i,t,b[i]);
        for (RI i=1;i<=n;++i)
        for (RI j=(int)a[i].size()-1;j>=0;--j) addedge(i,n+a[i][j],1);
    }
};
using namespace Network_Flow;
int main()
{
   	//freopen("1.in","r",stdin);
    for (scanf("%d",&T);T;--T)
    {
        scanf("%d%d",&n,&m); s=n+m+1; t=s+1;
        for (RI i=1;i<=n;++i) valid[i]=1;
        for (RI i=1;i<=m;++i) scanf("%d",&b[i]),rb[i]=0;
        for (RI i=1;i<=n;++i)
        {
            int x; scanf("%d",&x); a[i].resize(x);
            for (RI j=0;j<a[i].size();++j) scanf("%d",&a[i][j]);
        }
        build(); int flow=Dinic();
        printf("%d\n",flow);
        for (RI i=1;i<=n;++i)
        {
            mat[i]=0;
            for (RI j=head[i];j;j=e[j].nxt)
            if (1<=e[j].to-n&&e[j].to-n<=m&&e[j].v==0)
            {
                mat[i]=e[j].to-n; ++rb[mat[i]]; break;
            }
        }
        vector <int> ans;
        for (RI k=1;k<=flow;++k)
        {
            bool find=0;
            for (RI i=1;i<=n;++i)
            {
                if (!mat[i]||!valid[i]) continue;
                bool flag=1;
                for (auto x:a[i])
                {
                    if (x==mat[i]) break;
                    if (rb[x]>0) { flag=0; break; }
                }
                if (flag)
                {
                    ans.push_back(i); valid[i]=0;
                    --rb[mat[i]]; find=1; break;
                }
            }
            assert(find==1);
        }
        for (RI i=1;i<=n;++i) if (!mat[i]) ans.push_back(i);
        for (auto x:ans) printf("%d ",x); putchar('\n');
    }
    return 0;
}