from random import randint,random
import time
import sys
input=sys.stdin.readline
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(strict):
    pos=0
    go={0}
    ans=[0]
    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.append(next)
        pos=next
        ansdis+=dis
        go.add(pos)
    ans.append(0)
    temp=st_temp
    loop=0
    while strict:
        loop+=1
        l=randint(1,N-1)
        r=randint(1,N-1)
        while l==r:
            r=randint(1,N-1)
        if l>r:
            l,r=r,l
        pre=(dist[ans[l-1]][ans[l]][0]+dist[ans[r]][ans[r+1]][0])-(dist[ans[l-1]][ans[r]][0]+dist[ans[l]][ans[r+1]][0])
        if pre>0 or random()<pow(2.71,(pre/temp)):
            ansdis-=pre
            for i in range((r-l+1)//2):
                ans[l+i],ans[r-i]=ans[r-i],ans[l+i]
        temp*=c
        if temp<en_temp:
            break
    return ansdis,ans

st=time.time()
N,M=map(int,input().split())
star=[tuple(map(int,input().split()+["5"])) for i in range(N)]+[[0,0,1] for i in range(M)]
star[N]=[star[0][0],star[0][1],1]
st_temp=10**7
c=0.9995
en_temp=10**-6
for i in range(M-2):
    nea=0
    pos=-1
    for j in range(N):
        for k in range(i+1):
            pre=(star[N+k][0]-star[j][0])**2+(star[N+k][1]-star[j][1])**2
            if pre>nea:
                nea=pre
                pos=j
    star[N+i+1]=[star[pos][0],star[pos][1],1]
while True:
    gr=[[0,0,0] for j in range(M-1)]
    for j in range(N):
        nea=inf
        pos=-1
        for k in range(M-1):
            pre=(star[N+k][0]-star[j][0])**2+(star[N+k][1]-star[j][1])**2
            if pre<nea:
                nea=pre
                pos=k
        gr[pos][0]+=star[j][0]
        gr[pos][1]+=star[j][1]
        gr[pos][2]+=1
    ok=0
    for j in range(M-1):
        if gr[j][2]!=0:
            if star[N+j][0]!=gr[j][0]//gr[j][2]:
                ok=1
                star[N+j][0]=gr[j][0]//gr[j][2]
            if star[N+j][1]!=gr[j][1]//gr[j][2]:
                ok=1
                star[N+j][1]=gr[j][1]//gr[j][2]
    if ok==0:
        break
sumx=0
sumy=0
for i in range(M-1):
    sumx+=star[N+i][0]
    sumy+=star[N+i][1]
star[N+M-1]=[sumx//(M-1),sumy//(M-1),1]
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]]
for i in range(M):
    print(*star[N+i][:2])
ansdis,ans=sell(True)
finalans=[]
for i in range(len(ans)-1):
    for j in dist[ans[i]][ans[i+1]][1]:
        finalans.append(j)
finalans.append(0)
print(len(finalans))
for i in finalans:
    if i>=N:
        print(2,i-N+1)
    else:
        print(1,i+1)