#include<bits/stdc++.h>
using namespace std;
constexpr int maxn=2e5+5;
int n;vector<pair<int,int>> e[maxn];
int dep[maxn],tmp[maxn];
vector<pair<int,int>> ans;
void dfs(int p,int f){
	dep[p]=dep[f]+1;
	int z=0;
	for(auto [x,y]:e[p]){
		if(!dep[x])dfs(x,p);
		if(dep[x]>dep[p]){
			if(tmp[x]) ans.emplace_back(tmp[x],y),tmp[x]=0;
			else{
				if(z) ans.emplace_back(z,y),z=0;
				else z=y;
			}
		}
	}
	tmp[p]=z;
}
inline void solve(){
	unordered_map<int,int> p[2];
	int tot=0;
	auto f=[&](int op,int x){if(!p[op].count(x))p[op][x]=++tot;return p[op][x];};
	ans.clear();
	cin>>n;
	for(int i=1;i<=n;i++){
		int x;cin>>x;
		int a=f(0,x+i),b=f(1,x-i);
		// cout<<a<<' '<<b<<'\n';
		e[a].emplace_back(b,i),e[b].emplace_back(a,i);
	}
	for(int i=1;i<=tot;i++) if(!dep[i]) dfs(i,0);
	if(ans.size()==n/2){
		// cout<<"!!!";for(int i=1;i<=tot;i++) cout<<tmp[i]<<' ';cout<<'\n';
		cout<<"Yes\n";
		for(auto [x,y]:ans)cout<<x<<' '<<y<<'\n';
	}
	else cout<<"No\n";
	for(int i=1;i<=tot;i++) e[i].clear(),dep[i]=tmp[i]=0;
}
int main(){
	int t;cin>>t;
	while(t--) solve();
}