#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

const int N = 1e5 + 10, M = 3e6 + 10;

int n;
int a[N];
int h[N], e[M], ne[M], idx;
int match[N];
bool st[N];

void add(int a, int b)
{
	e[idx] = b, ne[idx] = h[a], h[a] = idx ++ ;
}

bool find(int u)
{
	for(int i=h[u]; i!=-1; i=ne[i])
	{
		int j = e[i];
		if(!st[j])
		{
			st[j] = true;
			if(match[j] == 0 || find(match[j]))
			{
				match[j] = u;
				return true;
			}
		}
	}
	
	return false;
}

int main()
{
	int T;
	scanf("%d", &T);
	while(T -- )
	{
		memset(h, -1, sizeof h);
		memset(match, 0, sizeof match);
		idx = 0;
		
		scanf("%d", &n);
		for(int i=1; i<=n; i++)
			scanf("%d", &a[i]);

		for(int i=1; i<=n; i++)
			for(int j=i+1; j<=n; j++)
			{
				if(abs(i-j)==abs(a[i]-a[j]))
					add(i, j), add(j, i);
			}

		int cnt = 0;
		for(int i=1; i<=n; i++)
		{
			memset(st, 0, sizeof st);
			if(find(i)) cnt ++ ;
		}
		
		if(cnt < n/2)
		{
			puts("No");
			continue;
		}
		puts("Yes");
		
		memset(st, 0, sizeof st);
		for(int i=1; i<=n; i++)
			if(!st[i])
			{
				st[i] = st[match[i]] = true;
				printf("%d %d\n", i, match[i]);
			}
	}
	
	return 0;
}