#include<bits/stdc++.h>
using namespace std;
using ll=long long;
template<typename T>
ostream& operator << (ostream &out,const vector<T>&x){
	if(x.empty())return out<<"[]";
	out<<'['<<x[0];
	for(int len=x.size(),i=1;i<len;i++)out<<','<<x[i];
	return out<<']';
}
template<typename T>
vector<T> ary(const T *a,int l,int r){
	return vector<T>{a+l,a+1+r};
}
template<typename T>
void debug(T x){
	cerr<<x<<'\n';
}
template<typename T,typename ...S>
void debug(T x,S ...y){
	cerr<<x<<' ',debug(y...);
}
const int N=1e3+10;
int T,n,fa[N],a[N],vis[N];
int find(int x){
	return fa[x]==x?x:fa[x]=find(fa[x]);
}
vector<tuple<int,int,int> >ans;
void merge(int x,int y){
	ans.push_back({2,x,y});
	x=find(x),y=find(y);
	if(x^y)fa[x]=y;
}
vector<int>to[N];
void dfs(int u){
	for(int v:to[u]){
		dfs(v);
	}
	if(fa[u]==u)for(int v:to[u]){
		merge(v,u);
	}
	if(fa[u]==u)for(int v:to[u]){
		find(v);
		ans.push_back({1,v,0});
	}
	// if(u==1){
		// debug(u,ary(fa,1,n));
	// }
}
bool get(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++)scanf("%d",&fa[i]);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		if(a[i]^i)to[a[i]].push_back(i);
	}
	fill(vis,vis+1+n,0);
	for(int i=1;i<=n;i++)if(a[i]^i){
		vis[i]=1;
		for(int j=a[i];;j=a[j]){
			if(j==i)return 0;
			if(vis[j])break;
			vis[j]=1;
			if(a[j]==j)break;
		}
	}
	ans.clear();
	for(int i=1;i<=n;i++)if(a[i]==i)dfs(i);
	// debug(ary(fa,1,n),ary(a,1,n));
	for(int i=1;i<=n;i++)if(fa[i]^a[i])return 0;
	puts("YES");
	printf("%d\n",(int)ans.size());
	for(auto x:ans){
		if(get<0>(x)==1)printf("%d %d\n",get<0>(x),get<1>(x));
		else printf("%d %d %d\n",get<0>(x),get<1>(x),get<2>(x));
	}
	return 1;
}
void clr(){
	for(int i=1;i<=n;i++)to[i].clear();
}
int main(){
	for(scanf("%d",&T);T--;clr())if(!get())puts("NO");
	return 0;
}