#include<bits/stdc++.h>
using namespace std;
const int N=200005;
int t,n,a[N],dfn[N];
vector<int>vx,vy;
vector<pair<int,int> >ans,e[N];
bool dfs(int k,int id){
    vector<int>v;
    dfn[k]=++dfn[0];
    for(pair<int,int>i:e[k]){
        int x=i.first,y=i.second;
        if (y==id)continue;
        if (!dfn[x]){
            if (dfs(x,y))v.push_back(y);
        }
        else{
            if (dfn[k]<dfn[x])v.push_back(y);
        }
    }
    bool flag=1;
    if (v.size()&1)flag=0,v.push_back(id);
    for(int i=0;i<v.size();i+=2)ans.push_back(make_pair(v[i],v[i+1]));
    return flag;
}
int main(){
    scanf("%d",&t);
    while (t--){
        scanf("%d",&n);
        vx.clear(),vy.clear();
        for(int i=1;i<=n;i++){
            scanf("%d",&a[i]);
            vx.push_back(a[i]-i);
            vy.push_back(a[i]+i);
        }
        sort(vx.begin(),vx.end());
        sort(vy.begin(),vy.end());
        for(int i=0;i<=(n<<1);i++)dfn[i]=0,e[i].clear();
        for(int i=1;i<=n;i++){
            int x=lower_bound(vx.begin(),vx.end(),a[i]-i)-vx.begin()+1;
            int y=lower_bound(vy.begin(),vy.end(),a[i]+i)-vy.begin()+1;
            e[x].push_back(make_pair(y+n,i));
            e[y+n].push_back(make_pair(x,i));
        }
        bool flag=0;
        ans.clear();
        for(int i=1;i<=(n<<1);i++)
            if (!dfs(i,0)){
                flag=1,puts("No");
                break;
            }
        if (!flag){
            puts("Yes");
            for(pair<int,int>i:ans)printf("%d %d\n",i.first,i.second);
        }
    }
    return 0;
}