#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=3005;
int T,n,m;
struct node{
	int x,y,id;
}A[N];
bool cmp(node a,node b){
	if(a.x==b.x) return a.y<b.y;
	return a.x<b.x;
}
bool vis[N];
int stk[N],tot;
bool check(int i,int j,int k){
	return 1ll*(A[j].y-A[i].y)*(A[k].x-A[j].x)<1ll*(A[k].y-A[j].y)*(A[j].x-A[i].x);
}
double get1(int i,int j,int k){
	double x1=A[i].x-A[j].x,y1=A[i].y-A[j].y;
	double x2=A[k].x-A[j].x,y2=A[k].x-A[j].y; 
	return (x1*x2+y1*y2)/sqrt(x1*x1+y1*y1)/sqrt(x2*x2+y2*y2);
}
int main(){
	scanf("%d",&T);
	while(T--){
		scanf("%d",&n);m=3*n;
		for(int i=1;i<=m;i++){
			scanf("%d%d",&A[i].x,&A[i].y);
			vis[i]=false;A[i].id=i;
		}
		sort(A+1,A+m+1,cmp);
		while(n--){
			tot=0;
			for(int i=1;i<=m;i++){
				if(vis[i]) continue;
				if(tot and A[i].x==A[stk[tot]].x) continue;
				while(tot>1 and !check(stk[tot-1],stk[tot],i)) tot--;
				stk[++tot]=i;
			}
			vis[stk[1]]=vis[stk[2]]=true;
			printf("%d %d ",A[stk[1]].id,A[stk[2]].id);
			int mi;double mm=-1;
			for(int i=1;i<=m;i++){
				if(vis[i]) continue;
				double mk=get1(stk[1],stk[2],i);
				if(mk>mm){
					mi=i;
					mm=mk;
				}
			}
			vis[mi]=true;
			printf("%d\n",A[mi].id);
		}
	}
	return 0;
}