#include <bits/stdc++.h>
using namespace std;
 
typedef long long ll;
const ll mod = 998244353;
const int maxn = 2e5 + 5;
const int N = 4e5 + 10;

int a[maxn], hd1[maxn], tl1[maxn], hd2[maxn], tl2[maxn], nxt[N], n;
int p1[maxn], p2[maxn];
vector<int> g[maxn];

void dfs(int u)
{
    sort(g[u].begin(), g[u].end());
    hd1[u] = tl1[u] = u, hd2[u] = tl2[u] = u + n;
    int tmph = 0, tmpt = 0;
    for (int v : g[u])
    {
        dfs(v);
        nxt[tl2[u]] = hd2[v], tl2[u] = v;
        if (tmph && nxt[v]) nxt[tmpt] = nxt[v], tmpt = tl2[v];
        else if (nxt[v]) tmph = nxt[v], tmpt = tl2[v];
        if (nxt[hd1[v]]) nxt[tl1[u]] = nxt[hd1[v]], tl1[u] = tl1[v];
    }
    if (tmph) nxt[tl2[u]] = tmph, tl2[u] = tmpt;
    else nxt[tl2[u]] = 0;
    reverse(g[u].begin(), g[u].end());
    for (int v : g[u]) nxt[hd1[v]] = hd1[u], hd1[u] = hd1[v];
    swap(hd1[u], hd2[u]), swap(tl1[u], tl2[u]);
}

// bool c[maxn];

void work()
{
    cin >> n;
    for (int i = 0; i <= n; i++) g[i].clear(), hd1[i] = hd2[i] = tl1[i] = tl2[i] = 0;
    for (int i = 1; i <= n; i++) cin >> a[i];
    for (int i = 1; i <= n; i++)
        if (a[i] > i) return cout << "No\n", void();
        else if (a[i] < i) g[a[i]].emplace_back(i);
    for (int i = 1; i <= 2 * n; i++) nxt[i] = 0;

    int t = 0;
    for (int i = 1; i <= n; i++)
        if (a[i] == i)
        {
            dfs(i);
            for (int j = hd1[i], k = hd2[i]; j; j = nxt[j], k = nxt[k]) p1[++t] = (j > n ? j - n : j), p2[t] = (k > n ? k - n : k);
        }
    cout << "Yes\n";
    for (int i = 1; i <= n; i++) cout << p1[i] << " ";
    for (int i = 1; i <= n; i++) cout << p2[i] << " ";
    cout << "\n";
}

int main()
{
    freopen("c.in", "r", stdin);
    ios::sync_with_stdio(0), cin.tie(0);

    int T; cin >> T;
    while (T--) work();
    return 0;
}