#include<bits/stdc++.h>
#define pb push_back
using namespace std;
const int N=2e5+5;
int n,A[N],Flg,D[N];
vector<int>L,R;
unordered_map<int,int>Ma[2];
struct Edge{int v,Id;};
vector<Edge>E[N],Ans;
int Dfs(int u,int f,int Lst)
{
	D[u]=D[f]+1;
	vector<int>e;
	for(Edge i:E[u])
	{
		if(i.v==f)continue;
		if(D[i.v])
		{
			if(D[i.v]<D[u])e.pb(i.Id);
		}
		else
		{
			int K=Dfs(i.v,u,i.Id);
			if(K!=-1)e.pb(K);
		}
	}
	while(e.size()>1)
	{
		int x=e.back();e.pop_back();
		int y=e.back();e.pop_back();
		Ans.pb({x,y});
	}
	if(e.size())
	{
		if(Lst==-1)Flg=0;
		Ans.pb({Lst,e.back()});
		return -1;
	}
	return Lst;
}
inline void Work()
{
	L.clear(),R.clear();
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&A[i]);
		L.pb(i-A[i]),R.pb(i+A[i]);
	}
	sort(L.begin(),L.end()),sort(R.begin(),R.end());
	L.erase(unique(L.begin(),L.end()),L.end());
	R.erase(unique(R.begin(),R.end()),R.end());
	int m1=L.size(),m2=R.size();
	for(int i=0;i<m1;i++)Ma[0][L[i]]=i;
	for(int i=0;i<m2;i++)Ma[1][R[i]]=i;
	for(int i=0;i<m1+m2;i++)E[i].clear();
	for(int i=1;i<=n;i++)
	{
		int l=Ma[0][i-A[i]],r=m1+Ma[1][i+A[i]];
		E[l].pb({r,i}),E[r].pb({l,i});
	}
	Flg=1,Ans.clear();
	memset(D,0,sizeof(D));
	for(int i=0;i<m1+m2;i++)if(!D[i])Dfs(i,0,-1);
	if(!Flg)return puts("No"),void();
	puts("Yes");
	for(Edge i:Ans)printf("%d %d\n",i.v,i.Id);
}
int main()
{
	int T;scanf("%d",&T);
	while(T--)Work();
}