# include<bits/stdc++.h>
using namespace std;
typedef long long int ll;

ll b[1000000+10],lisan[1000000+10],l[1000000+10],r[1000000+10],a[1000000+10];
int book[1000000+10];
vector<pair<int,int>>v[1000000+10];
vector<int>ans;
int cnt;
int du[1000000+10];
void dfs(int now)
{


    for(int i=v[now].size()-1;i>=0;i--)
    {
        if(book[v[now][i].second])
            continue;
        book[v[now][i].second]=1;
        ans.push_back(v[now][i].second);
        cnt++;
        dfs(v[now][i].first);
        v[now].pop_back();

    }
}
int main ()
{

    int t;
    cin>>t;

    while(t--)
    {
        int n;
        scanf("%d",&n);
        int len=0;
        ans.clear();
        for(int i=1; i<=n; i++)
        {
            scanf("%lld",&a[i]);
            len++;
            lisan[len]=a[i]+i;
            len++;
            lisan[len]=a[i]-i;
            l[i]=a[i]+i;
            r[i]=a[i]-i;
        }


        sort(lisan+1,lisan+1+2*n);
        len=unique(lisan+1,lisan+1+2*n)-lisan-1;

        for(int i=1; i<=len; i++)
        {
            v[i].clear();
            book[i]=0;
        }

      //  cout<<endl;

        for(int i=1; i<=n; i++)
        {
            l[i]=lower_bound(lisan+1,lisan+1+len,l[i])-lisan;
            r[i]=lower_bound(lisan+1,lisan+1+len,r[i])-lisan;

            v[l[i]].push_back(make_pair(r[i],i));
            v[r[i]].push_back(make_pair(l[i],i));

           // cout<<l[i]<<" "<<r[i]<<endl;
        }


      //  cout<<endl;

        int flag=0;
        for(int i=1; i<=n; i++)
        {
            cnt=0;

          dfs(l[i]);

          if(cnt%2==1)
          {
              flag=1;
              break;
          }


           // cout<<cnt<<endl;

        }

        if(flag)
        {
            cout<<"No"<<'\n';
        }
        else
        {

              cout<<"Yes"<<'\n';
             for(int i=0;i<ans.size();i++)
             {
                 cout<<ans[i]<<" "<<ans[i+1]<<'\n';
                 i++;
             }

             cout<<'\n';
        }
    }
    return 0;
}