题目链接：戳我

我们先做一个DP，就能求出来到前i位的最长上升子序列的长度（maxx[i]数组）。

然后我们考虑求最小割——给每个点拆点，如果要割掉这个点，就相当于把in[i]--out[i]这条边给割掉了。

然后如果在最长上升子序列中，该位下面可以接很多，那么就从该位代表的out向下面可以接的in连边。

然后跑最小割就行了QAQ

但是......之后还需要求字典序最小的最小割方案。
那么我们先把边从小到大排个序，然后把所有可能割边都捞出来，贪心地选择。
但是有一个问题，就是我们选择了一个边之后，它会对我们之后的决策产生影响（比如说，本来应该选的边可以不选了）
所以我们要强行退流（具体操作可以看代码）。我们知道选了这条边之后，这条边的可替代品就不需要选了，这些可替代它的边显然是存在于S-u,v-T之间的。所以我们直接退掉他们满流边的流量就行啦！

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#define S 0
#define T 2*n+1
#define MAXN 1510
#define INF 0x3f3f3f3f
using namespace std;
inline int read()
{int x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch==-1)f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=(x<<1)+(x<<3)+(ch^48);ch=getchar();}return x*f;
}
int TT,n,t=1,cnt,kkk;
int num[MAXN];
int ans[MAXN],in[MAXN],out[MAXN],done[MAXN];
int head[MAXN],cur[MAXN],dis[MAXN],dp[MAXN][MAXN],maxx[MAXN];
struct Node{int id,a,b,c;}node[MAXN];
struct Edge{int nxt,to,dis;}edge[MAXN*MAXN*2];
inline void add(int from,int to,int dis)
{
//  printf("[%d %d] %d\n",from,to,dis);edge[++t].nxt=head[from],edge[t].to=to,edge[t].dis=dis,head[from]=t;edge[++t].nxt=head[to],edge[t].to=from,edge[t].dis=0,head[to]=t;
}
inline bool check(int s,int t)
{queue<int>q;memset(done,0,sizeof(done));q.push(s);done[s]=1;while(!q.empty()){int u=q.front();q.pop();if(u==t) return true;for(int i=head[u];i;i=edge[i].nxt){int v=edge[i].to;if(edge[i].dis&&done[v]==0) q.push(v),done[v]=1;}}return false;
}
inline bool bfs(int s,int t)
{queue<int>q;memset(dis,0x3f,sizeof(dis));memcpy(cur,head,sizeof(head));q.push(s);dis[s]=0;while(!q.empty()){int u=q.front(); q.pop();for(int i=head[u];i;i=edge[i].nxt){int v=edge[i].to;if(edge[i].dis&&dis[v]==0x3f3f3f3f){dis[v]=dis[u]+1;q.push(v);}}}if(dis[t]==0x3f3f3f3f) return false;return true;
}
inline int dfs(int x,int t,int f)
{if(!f||x==t) return f;int used=0,w;for(int i=cur[x];i;i=edge[i].nxt){int v=edge[i].to;cur[x]=i;if(dis[v]==dis[x]+1&&(w=dfs(v,t,min(edge[i].dis,f)))){edge[i].dis-=w;edge[i^1].dis+=w;used+=w,f-=w;if(!f) break;}}return used;
}
inline int dinic(int s,int t)
{int cur_ans=0;while(bfs(s,t)) cur_ans+=dfs(s,t,INF);return cur_ans;
}
inline void pre()
{for(int i=1;i<=n;i++){maxx[i]=1;for(int j=1;j<i;j++)if(node[i].a>node[j].a)maxx[i]=max(maxx[i],maxx[j]+1);}   for(int i=1;i<=n;i++) kkk=max(kkk,maxx[i]);
}
inline void add_edge()
{for(int i=1;i<=n;i++) in[i]=i;for(int i=n+1;i<=2*n;i++) out[i-n]=i;for(int i=1;i<=n;i++) add(in[i],out[i],node[i].b),num[i]=t-1;for(int i=1;i<=n;i++){for(int j=1;j<i;j++){if(maxx[j]+1==maxx[i]&&node[j].a<node[i].a)add(out[j],in[i],INF);}}for(int i=1;i<=n;i++) if(maxx[i]==1) add(S,in[i],INF);for(int i=1;i<=n;i++) if(maxx[i]==kkk) add(out[i],T,INF);
}
inline bool cmp(struct Node x,struct Node y){return x.c<y.c;}
inline void solve()
{for(int i=1;i<=n;i++){int x=node[i].id;if(check(in[x],out[x])) continue;ans[++cnt]=x;while(bfs(T,out[x])) dfs(T,out[x],INF);while(bfs(in[x],S)) dfs(in[x],S,INF);edge[num[x]].dis=0,edge[num[x]^1].dis=0;}
}
int main()
{#ifndef ONLINE_JUDGEfreopen("ce.in","r",stdin);#endifTT=read();while(TT--){memset(head,0,sizeof(head));cnt=0,kkk=0;t=1;n=read();for(int i=1;i<=n;i++) node[i].a=read();for(int i=1;i<=n;i++) node[i].b=read();for(int i=1;i<=n;i++) node[i].c=read();pre();add_edge();printf("%d ",dinic(S,T));for(int i=1;i<=n;i++) node[i].id=i;sort(&node[1],&node[n+1],cmp);solve();printf("%d\n",cnt);sort(&ans[1],&ans[cnt+1]);for(int i=1;i<=cnt;i++) printf("%d ",ans[i]); puts("");}return 0;
}