#include <bits/stdc++.h>
using namespace std;
const int INF = 2e9;
const int MOD = 1e9 + 7;
const int MAX = 5e5 + 10;
const long long LNF = 2e18;
class Seg_t {
  int n;
  struct node {
    pair<int, int> mn;
    int lz;
  };
  const node E = {{INF, -1}, 0};
  vector<node> T;
  void busy(int v, int s, int e) {
    if (s == e) return;
    int z = T[v].lz;
    T[v].lz = 0;
    T[v * 2].lz += z;
    T[v * 2].mn.first += z;
    T[v * 2 + 1].lz += z;
    T[v * 2 + 1].mn.first += z;
  }
  void upt(int v, int s, int e, int l, int r, int x) {
    if (e < l || r < s) return;
    node &nd = T[v];
    if (l <= s && e <= r) {
      nd.lz += x;
      nd.mn.first += x;
      return;
    }
    busy(v, s, e);
    upt(v * 2, s, (s + e) / 2, l, r, x);
    upt(v * 2 + 1, (s + e) / 2 + 1, e, l, r, x);
    T[v].mn = min(T[v * 2].mn, T[v * 2 + 1].mn);
  }
  void init(int v, int s, int e) {
    if (s == e) {
      T[v] = {{0, -s}, 0};
      return;
    }
    init(v * 2, s, (s + e) / 2);
    init(v * 2 + 1, (s + e) / 2 + 1, e);
    T[v].mn = min(T[v * 2].mn, T[v * 2 + 1].mn);
  }

 public:
  void init(int sz) {
    T.clear();
    n = sz;
    int k = 1;
    while (k <= sz) k *= 2;
    T.resize(k * 2, E);
    init(1, 1, n);
  }
  void add(int l, int r, int x) {
    if (l > r) return;
    upt(1, 1, n, l, r, x);
  }
  void add(int p, int x) { add(p, p, x); }
  pair<int, int> get() {
    pair<int, int> res = T[1].mn;
    res.second *= -1;
    return res;
  }
} Seg;
vector<int> G[MAX];
int main() {
  ios::sync_with_stdio(0);
  cin.tie(0);
  int n, X[MAX], P[MAX];
  int t;
  cin >> t;
  for (int tt = 1; tt <= t; tt++) {
    cin >> n;
    Seg.init(n + 1);
    for (int i = 1; i <= n; i++) cin >> X[i];
    X[n + 1] = n + 1;
    for (int i = 1; i <= n + 1; i++) G[i].clear();
    for (int i = 1; i <= n; i++)
      if (X[i] > 0) G[X[i]].push_back(i);
    for (int i = 1; i <= n + 1; i++) {
      if (X[i] > 0) Seg.add(i, X[i] - 1, 1);
    }
    bool fail = false;
    for (int val = n + 1; val >= 1; val--) {
      pair<int, int> now = Seg.get();
      int deg, idx;
      tie(deg, idx) = now;
      if (deg != 0) {
        fail = true;
        break;
      }
      P[idx] = val;
      if (X[idx] > 0) Seg.add(idx + 1, X[idx] - 1, -1);
      Seg.add(idx, INF);
      for (int x : G[idx]) Seg.add(x, -1);
    }
    if (fail) {
      cout << "-1\n";
      continue;
    }
    for (int i = 1; i <= n; i++) cout << P[i] << ' ';
    cout << '\n';
  }
  return 0;
}
