#include <bits/stdc++.h>
using namespace std;
map<int,int> mp;
vector< pair<int,int> > ans;
int T,n,a[100010],b[200010],Link[200010],dfn[200010],t,tot,sum;
bool vis[200010];
struct edge
{
	int v,nex,id;
}e[200010];
void Insert(int x,int y,int z)
{
	e[++t].nex=Link[x];e[t].v=y;Link[x]=t;e[t].id=z;return ;
}
void dfs(int now,int las)
{
	dfn[now]=++tot;vis[now]=false;
	vector <int> V;
	for (int i=Link[now];i;i=e[i].nex) {
		if (i==(las^1)) continue;
		if (!dfn[e[i].v]) {
			dfs(e[i].v,i);
			sum++;
			if (!vis[e[i].v]) 
			  V.push_back(e[i].id);
		}
		else if (dfn[e[i].v]<dfn[now]) {
			V.push_back(e[i].id);
			sum++;
		}
	}
	if ((int)V.size()&1) {
		for (int i=1;i<(int)V.size()-1;i++)
		  ans.push_back(make_pair(V[i],V[i+1]));
		ans.push_back(make_pair(V[0],e[las].id));
		vis[now]=true;
	}
	else 
	  for (int i=0;i<(int)V.size()-1;i++)
		ans.push_back(make_pair(V[i],V[i+1]));
	return ;
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	cin>>T;
	for (int i1=1;i1<=T;i1++) {
		cin>>n;int nu=0;
		for (int i=1;i<=n;i++) {
		    cin>>a[i];
		    b[++nu]=a[i]-i;
		    b[++nu]=a[i]+i;
		}
		sort(b+1,b+nu+1);
		int cnt=0;
		mp.clear();
		for (int i=1;i<=nu;i++)
		  if (mp.find(b[i])==mp.end()) mp[b[i]]=++cnt;
		t=1;
		for (int i=1;i<=cnt;i++) Link[i]=0;
		for (int i=1;i<=n;i++)
		  Insert(mp[a[i]-i],mp[i+a[i]],i),Insert(mp[a[i]+i],mp[a[i]-i],i);
		bool fla=true;ans.clear();
		for (int i=1;i<=cnt;i++) dfn[i]=0;
		for (int i=1;i<=cnt;i++)
		  if (!dfn[i]) {
		  	  tot=0;sum=0;dfs(i,0);
		  	  if (sum&1) {fla=false;break;}
		  }
		if (!fla) cout<<"No"<<'\n';
		  else {
		  	  cout<<"Yes"<<'\n';
		  	  for (int i=0;i<(int)ans.size();i++)
		  	    cout<<ans[i].first<<' '<<ans[i].second<<'\n';
		  }
	}
	return 0;
}