#include<bits/stdc++.h>
using namespace std;
int n;
int a[100005],b[2][100005];
int vis[400005],head[200005],tot;
int h[400005];
struct ss{
	int node,nxt;
}e[400005];
void add(int u,int v)
{
	e[++tot].nxt=head[u];
	e[tot].node=v;
	head[u]=tot;
}
pair<int,int> ans[100005];
int len;
void dfs(int x,int fa)
{
//	cerr<<x<<endl;
	int now=0;
	int f=0;
	h[x]=1;
	for(int i=head[x];i;i=e[i].nxt)
	{
		int y=e[i].node;
		if(vis[i])continue;
		if(y==fa)
		{
			f=i;
			continue;
		}
		if(!h[y]) 
			dfs(y,x);
		if(now)
		{
			ans[++len]=make_pair(now>>1,i>>1);
			vis[i]=vis[i^1]=vis[now]=vis[now^1]=1;
			now=0;
		}
		else now=i;
	}
	if(now&&f)
	{
		ans[++len]=make_pair(now>>1,f>>1);
		vis[now]=vis[now^1]=vis[f]=vis[f^1]=1;
	}
}
signed main()
{
	ios::sync_with_stdio(0);
	cin.tie(0);cout.tie(0);
	int T;
	cin>>T;
//	scanf("%d",&T);
	while(T--)
	{
		cin>>n;
		tot=1;
		for(int i=1;i<=n;i++)
		{
			cin>>a[i];
			b[0][i]=i+a[i];
			b[1][i]=i-a[i];
		}
		sort(b[0]+1,b[0]+n+1);
		sort(b[1]+1,b[1]+n+1);
		b[0][0]=unique(b[0]+1,b[0]+n+1)-b[0]-1;
		b[1][0]=unique(b[1]+1,b[1]+n+1)-b[1]-1;
		for(int i=1;i<=b[0][0]+b[1][0];i++)
			h[i]=head[i]=0;
		for(int i=1;i<=n;i++)
		{
			int x=lower_bound(b[0]+1,b[0]+b[0][0]+1,a[i]+i)-b[0];
			int y=lower_bound(b[1]+1,b[1]+b[1][0]+1,i-a[i])-b[1];
		//	cerr<<x<<" "<<y<<endl;
			add(x,y+b[0][0]);add(y+b[0][0],x); 
		}
		for(int i=1;i<=tot;i++)
			vis[i]=0;
		len=0;
		for(int i=1;i<=b[0][0]+b[1][0];i++)
			if(!h[i])dfs(i,0);
		if(n==100000)cout<<tot<<" "<<len<<endl;
		if(len==n/2)
		{
			cout<<"Yes\n";
			for(int i=1;i<=len;i++)
				cout<<ans[i].first<<" "<<ans[i].second<<'\n';
		}
		else cout<<"No\n";
	}
}
/*
1
3 18 33
UDRLR
*/