#include<cstdio>
#include<map>
using namespace std;
const int o=2e5+10;
int n,T,a[o],tot,h[o],cnt,d[o],f[o],ansx[o],ansy[o],ans,st[o],tp;bool flg;map<int,int> mp1,mp2;
struct Edge{int v,p;}e[o];
inline void ad(int U,int V){e[++cnt].v=V;e[cnt].p=h[U];h[U]=cnt;}
void dfs(int nw){
	for(int i=h[nw];i;i=e[i].p) if(!d[e[i].v]) d[e[i].v]=d[nw]+1,dfs(e[i].v);
	tp=0;
	for(int i=h[nw];i;i=e[i].p)
		if(d[e[i].v]==d[nw]+1&&f[e[i].v]) ansx[++ans]=i/2,ansy[ans]=f[e[i].v];
		else if(d[e[i].v]>d[nw]) st[++tp]=i/2;
	for(int i=2;i<=tp;i+=2) ansx[++ans]=st[i-1],ansy[ans]=st[i];https://qoj.ac/contest/1093/problem/5423/statistics
	if(tp&1) f[nw]=st[tp];
}
int main(){
	for(scanf("%d",&T);T--;mp1.clear(),mp2.clear(),tot=ans=flg=0){
		scanf("%d",&n);cnt=1;
		for(int i=1,x,y;i<=n;++i){
			scanf("%d",&a[i]);
			if(mp1.find(a[i]-i)==mp1.end()) mp1[a[i]-i]=++tot;
			if(mp2.find(a[i]+i)==mp2.end()) mp2[a[i]+i]=++tot;
			ad(x=mp1[a[i]-i],y=mp2[a[i]+i]);ad(y,x);
		}
		for(int i=1;i<=tot;++i) if(!d[i]){d[i]=1;dfs(i);if(f[i]){flg=1;break;}}
		if(flg) printf("No\n");
		else{
			printf("Yes\n");
			for(int i=1;i<=n/2;++i) printf("%d %d\n",ansx[i],ansy[i]);
		}
		for(int i=1;i<=tot;++i) h[i]=d[i]=f[i]=0;
	}
	return 0;
}