#include <bits/stdc++.h>
using namespace std;

#define int long long

const int N=1e5+5;

int T,n,a[N],dep[N];
map<int,int> x,y;
vector<pair<int,int>> ans,G[N];

int dfs(int u,int pa)
{
    dep[u]=dep[pa]+1;
    vector<int> vec;
    for(auto i:G[u]){
        int v=i.first,id=i.second;
        if(v==pa) continue;
        if(dep[v]){
            if(dep[v]>dep[u])
                vec.emplace_back(id);
            continue;
        }
        int res=dfs(v,u);
        if(res) ans.emplace_back(id,res);
        else vec.emplace_back(id);
    }
    int res=0;
    if(vec.size()&1) res=vec.back(),vec.pop_back();
    for(int i=0;i<vec.size();i+=2)
        ans.emplace_back(vec[i],vec[i+1]);
    return res;
}

void solve()
{
    x.clear(),y.clear(),ans.clear();
    memset(dep,0,sizeof(dep));
    cin>>n;
    for(int i=1;i<=n;i++)
        cin>>a[i];
    for(int i=1;i<=n;i++)
        x[a[i]+i]=1,y[a[i]-i]=1;
    int cnt=0;
    for(auto &i:x) i.second=++cnt;
    for(auto &i:y) i.second=++cnt;
    for(int i=1;i<=cnt;i++)
        G[i].clear();
    for(int i=1;i<=n;i++){
        int u=x[a[i]+i],v=y[a[i]-i];
        G[u].emplace_back(v,i);
        G[v].emplace_back(u,i);
    }
    for(int i=1;i<=cnt;i++){
        if(dep[i]) continue;
        if(dfs(i,0))
            return cout<<"No\n",void();
    }
    cout<<"Yes\n";
    for(auto i:ans)
        cout<<i.first<<' '<<i.second<<'\n';
}

signed main()
{
    ios::sync_with_stdio(0),cin.tie(0),cout.tie(0);
    cin>>T;
    while(T--)
        solve();
    return 0;
}