from random import randint
inf=998244353

def calc(a,b):
    return ((star[a][0]-star[b][0])**2+(star[a][1]-star[b][1])**2)*star[a][2]*star[b][2]

def sell():
    pos=0
    go={0}
    ans=[]
    ansdis=0
    while len(go)<N:
        next=-1
        dis=inf
        for i in range(N):
            if i not in go:
                if dist[pos][i][0]<dis:
                    next=i
                    dis=dist[pos][i][0]
        ans+=dist[pos][next][1]
        pos=next
        ansdis+=dis
        go.add(pos)
    ans+=dist[pos][0][1]
    return ansdis,ans+[0]

N,M=map(int,input().split())
star=[tuple(map(int,input().split()+["5"])) for i in range(N)]+[[randint(0,1000),randint(0,1000),1] for i in range(M)]
prest=[]
preans=[]
predis=inf
t=7
for _ in range(t):
	dist=[[[calc(i,j),[j]] for i in range(N+M)]for j in range(N+M)]
	for i in range(N+M):
		for j in range(N+M):
			for k in range(N+M):
				if dist[j][i][0]+dist[i][k][0]<dist[j][k][0]:
					dist[j][k]=[dist[j][i][0]+dist[i][k][0],dist[j][i][1]+dist[i][k][1]]
	ansdis,ans=sell()
	if predis>ansdis:
		predis=ansdis
		preans=ans[:]
		prest.clear()
		for i in range(M):
			prest.append(star[N+i][:2])
for i in range(M):
	print(*prest[i])
print(len(preans))
for i in preans:
	if i>=N:
		print(2,i-N+1)
	else:
		print(1,i+1)