# 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];
set<pair<int,int>>v[1000000+10];
vector<int>ans;
int cnt;
int du[1000000+10];
vector<pair<int,pair<int,int>>>temp;
void dfs(int now)
{
   for(auto it:v[now])
   {

       if(book[it.second])
        continue;
       book[it.second]=1;
       temp.push_back(make_pair(it.first,make_pair(now,it.second)));
       dfs(it.first);
       cnt++;
       ans.push_back(it.second);
       temp.push_back(make_pair(now,it));


   }
}
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]].insert(make_pair(r[i],i));
            v[r[i]].insert(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;

            temp.clear();
          dfs(l[i]);

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


          for(auto it:temp)

          {
              v[it.first].erase(it.second);
          }
           // 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;
}