/*

								_/                                      _/
	   _/_/        _/_/_/      _/_/_/        _/_/_/      _/  _/_/               _/_/_/
	_/    _/    _/            _/    _/    _/    _/      _/_/          _/       _/    _/
   _/    _/    _/            _/    _/    _/    _/      _/            _/       _/    _/
	_/_/        _/_/_/      _/    _/      _/_/_/      _/            _/       _/    _/

*/
#include<bits/stdc++.h>
#define int long long

using namespace std;

void solve(){
	int n;
	cin>>n;
	vector<int>a(n);
	for(int i=0;i<n;++i) cin>>a[i];
	unordered_map<int,int>mp1,mp2;
	int c1=0,c2=n;
	vector<vector<array<int,2>>>e(2*n);
	for(int i=0;i<n;++i){
		int u=i-a[i],v=i+a[i];
		if(!mp1.count(u)) mp1[u]=c1++;
		if(!mp2.count(v)) mp2[v]=c2++;
		u=mp1[u],v=mp2[v];
		e[u].push_back({v,i+1});
		e[v].push_back({u,i+1});
	}
	vector<array<int,2>>res;
	vector<int>dep(2*n),vis(2*n);
	auto dfs=[&](auto dfs,int u)->int{
		vis[u]=1;
		int cur=0;
		for(auto [v,w]:e[u]){
			if(vis[v]){
				if(dep[v]<=dep[u]) continue;
				if(!cur){
					cur=w;
				}
				else{
					res.push_back({cur,w});
					cur=0;
				}
			}
			else{
				dep[v]=dep[u]+1;
				int x=dfs(dfs,v);
				if(x){
					res.push_back({x,w});
				}
				else{
					if(cur){
						res.push_back({w,cur});
						cur=0;
					}
					else cur=w;
				}
			}
		}
		return cur;
	};
	for(int i=0;i<2*n;++i){
		if(!vis[i]) dfs(dfs,i);
	}
	if(res.size()*2==n){
		cout<<"Yes\n";
		for(auto [x,y]:res) cout<<x<<" "<<y<<"\n";
	}
	else cout<<"No\n";
}

signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0);
	int T;cin>>T;
	while(T--) solve();
}