#include <iostream>
#include <cstdio>
#include <vector>
#include <map>
using namespace std;
int t,n,a[200005];
int tot1,tot2;
map<int,int> id1,id2;
vector<int> e[400005],e_id[400005];
int depth[400005];
int ans[200005][2],tot;
void add(int x,int y){
	++tot;
	ans[tot][0]=x,ans[tot][1]=y;
	return;
}
int dfs(int now,int fa,int d){
	depth[now]=d;
	int qwq=-1;
	for (int i=0;i<(int)e[now].size();i++){
		int v=e[now][i],w=e_id[now][i];
		if (v==fa)continue;
		if (depth[v]==0){
			int t=dfs(v,now,d+1);
			if (t!=-1)add(t,w);
			else if (qwq!=-1)add(qwq,w),qwq=-1;
			else qwq=w;
		}
		else if (depth[v]<depth[now]){
			if (qwq!=-1)add(qwq,w),qwq=-1;
			else qwq=w;
		}
	}
	return qwq;
}
int main(){
	scanf("%d",&t);
	while(t--){
		scanf("%d",&n);
		tot1=tot2=0;
		id1.clear(),id2.clear();
		for (int i=1;i<=2*n;i++){
			e[i].clear();
			e_id[i].clear();
			depth[i]=0;
		}
		for (int i=1;i<=n;i++){
			scanf("%d",&a[i]);
			if (id1[a[i]-i]==0)id1[a[i]-i]=++tot1;
			if (id2[a[i]+i]==0)id2[a[i]+i]=++tot2;
			int x=id1[a[i]-i],y=id2[a[i]+i]+n;
			e[x].push_back(y);
			e_id[x].push_back(i);
			e[y].push_back(x);
			e_id[y].push_back(i);
		}
		tot=0;
		for (int i=1;i<=2*n;i++)
			if (depth[i]==0)dfs(i,0,1);
		if (tot==n/2){
			printf("Yes\n");
			for (int i=1;i<=tot;i++)printf("%d %d\n",ans[i][0],ans[i][1]);
		}
		else printf("No\n");
	}
	return 0;
}