#include<bits/stdc++.h>
using namespace std;
#define rep(i,x,y) for (int i=(x);i<=(y);i++)
#define drep(i,y,x) for (int i=(y);i>=(x);i--)
#define pii pair<int,int>
#define fir first
#define sec second
#define MP make_pair
template<typename T> bool chkmin(T &x,T y){return x>y?x=y,1:0;}
template<typename T> bool chkmax(T &x,T y){return x<y?x=y,1:0;}
void file() {
    #ifdef zqj
    freopen("a.in","r",stdin);
    #endif
}
typedef long long ll;
#define sz 233333

int n;
int a[sz];

void work() {
    cin>>n;
    vector<int>V;
    static int ans[sz];
    rep(i,1,n) {
        cin>>a[i];
        if (a[i]==i) V.push_back(i);
        else ans[i]=i;
    }
    if (V.size()==1u) {
        int x=V[0];
        int ok=-1;
        rep(i,1,n) if (a[i]!=x) ok=i;
        if (ok==-1) return cout<<"NO\n",void();
        ans[x]=ok,ans[ok]=x;
    }
    else {
        rep(i,0,(int)V.size()-1) ans[V[i]]=V[(i+1)%V.size()];
    }
    cout<<"YES\n";
    rep(i,1,n) cout<<ans[i]<<" \n"[i==n];
}

int main(){
    file();
    ios::sync_with_stdio(false),cin.tie(0);
    int T; cin>>T;
    while (T--) work();
    return 0;
}