#include <bits/stdc++.h>
using namespace std;
const int MAX = 5e5;
int SPF[MAX + 5];
void update(int i, int val, vector<int> &BIT)
{
for (; i < BIT.size(); i += (i & (-i)))
BIT[i] += val;
}
int query(int i, vector<int>& BIT)
{
int res = 0;
for (; i > 0; i -= (i & (-i)))
res += BIT[i];
return res;
}
int solve(vector< pair<int, int> > &X, vector<int> &F)
{
vector<int> BIT(X.size() + 1);
for (int i = 1; i <= X.size(); i++)
update(i, 1, BIT);
for (auto &pr: X)
{
int indx = pr.first, target = pr.second + 1;
int left = 1, right = X.size(), res = -1;
while (left <= right)
{
int mid = (left + right) / 2;
if (query(mid, BIT) >= target)
res = mid, right = mid - 1;
else
left = mid + 1;
}
if (res == -1)
return 1;
F[indx] = res;
update(res, -1, BIT);
}
return 0;
}
int main() {
ios::sync_with_stdio(0);
cin.tie(NULL), cout.tie(NULL);
int t;
cin >> t;
SPF[1] = 1;
for (int i = 2; i <= MAX; i++)
if (SPF[i] == 0)
for (int j = i; j <= MAX; j += i)
if (SPF[j] == 0)
SPF[j] = i;
while (t--)
{
int n;
cin >> n;
vector< vector< pair<int, int> > > X(n + 1);
for (int i = 1; i <= n; i++)
{
int c;
cin >> c;
X[ SPF[i] ].push_back({i, c});
}
int f = 0;
vector<int> F(n + 1);
for (int i = 1; i <= n; i++)
if (!X[i].empty())
f |= solve(X[i], F);
if (f)
{
cout << "No\n";
continue;
}
cout << "Yes\n";
for (int i = 1; i <= n; i++)
cout << F[i] << " ";
cout << "\n";
}
return 0;
}