#include <bits/stdc++.h>
using namespace std;
struct UnionFind{
    vector<int> parent;
    vector<int> numberOfElements;
    UnionFind(int n){
        parent.resize(n);
        numberOfElements.resize(n);
        for(int i=0;i<n;++i){
            parent[i]=-1;
            numberOfElements[i]=1;
        }
    }
    int find(int a){
        vector<int> path(0);
        int pointer=a;
        while(parent[pointer]!=-1){
            path.push_back(pointer);
            pointer=parent[pointer];
        }
        for(int i:path){
            parent[i]=pointer;
        }
        return pointer;
    }
    void unite(int a,int b){
        int p=find(a);
        int q=find(b);
        if(p==q){
            return;
        }
        if(numberOfElements[p]<numberOfElements[q]){
            parent[p]=q;
            numberOfElements[q]+=numberOfElements[p];
        }else{
            parent[q]=p;
            numberOfElements[p]+=numberOfElements[q];
        }
    }
};

void solve(){
    int N,input;
    cin>>N;
    map<list<int>,int> m;
    UnionFind U(N);
    for(int j=2;j<=N;++j){
        for(int k=j*2;k<=N;k+=j){
            U.unite(j-1,k-1);
        }
    }
    for(int j=0;j<N;++j){
        cin>>input;
        ++m[{U.find(j),input}];
    }
    for(int j=0;j<N;++j){
        cin>>input;
        --m[{U.find(j),input}];
    }
    for(auto j=m.begin();j!=m.end();++j){
        if(j->second){
            cout<<"No"<<endl;
            return;
        }
    }
    cout<<"Yes"<<endl;
}

int main(){
    int T;
    cin>>T;
    for(int i=0;i<T;++i){
        solve();
    }
}