#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define MP make_pair
const int MAXN=1e6+5;
map<ll,int> M;int tot=0;
int ls[MAXN],rs[MAXN];
vector<pair<int,int> > ve[2*MAXN];// 到达节点，边编号 
vector<int> ans[MAXN];
int n;int a[MAXN];bool vis[MAXN];
bool dfs(int u,int fa){// 是否用完 
	vis[u]=1;
	for(auto i:ve[u]){
		if(i.first==fa) continue;
		if(vis[i.first]) ans[u].push_back(i.second);
	}
	for(auto i:ve[u]){
		if(vis[i.first]) continue;
		if(dfs(i.first,u)) ans[i.first].push_back(i.second);
		else ans[u].push_back(i.second);
	}
	return ans[u].size()&1; 
}
void solve(){
	cin>>n;
	for(int i=1;i<=n;i++){
		cin>>a[i];
	}
	M.clear();tot=0;
	for(int i=1;i<=n;i++){
		if(!M[a[i]+i])
			M[a[i]+i]=++tot;
		ls[i]=M[a[i]+i];
	}
	M.clear();
	for(int i=1;i<=n;i++){
		if(!M[a[i]-i])
			M[a[i]-i]=++tot;
		rs[i]=M[a[i]-i];
	}
	for(int i=1;i<=tot;i++) ve[i].clear();
	for(int i=1;i<=n;i++){
		ve[ls[i]].push_back(MP(rs[i],i));
		ve[rs[i]].push_back(MP(ls[i],i));
	}
	for(int i=1;i<=tot;i++) ans[i].clear(),vis[i]=0;
	for(int i=1;i<=tot;i++)
		if(!vis[i]) dfs(i,0);
	for(int i=1;i<=tot;i++)
		if(ans[i].size()&1){
			cout<<"No"<<endl;
			return;
		}
	cout<<"Yes"<<endl;
	for(int i=1;i<=tot;i++)
		for(int j=0;j<ans[i].size();j+=2)
			cout<<ans[i][j]<<" "<<ans[i][j+1]<<endl;
	return; 
}
int main(){
//	freopen("input.in","r",stdin);
//	freopen("answer.out","w",stdout);
	ios::sync_with_stdio(false);
	int _;cin>>_;
	while(_--) solve();
}