from random import randint,random
import time
import sys
input=sys.stdin.readline
inf=pow(2,61)-1

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)
    st_temp=10**7
    c=0.99
    if strict:
    	c+=0.009
    en_temp=10**-2
    ans.append(0)
    temp=st_temp
    loop=0
    while True:
        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.72,(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
    if strict:
        return ans
    return ansdis

def change(x,y,p):
    global star,dist
    star[N+p][0]+=x
    star[N+p][1]+=y
    for i in range(N+M):
        dist[i][N+p]=[calc(i,N+p),[i]]
        dist[N+p][i]=[calc(i,N+p),[N+p]]
    for i in range(N+M):
        for j in range(N+M):
            if dist[i][N+p][0]+dist[N+p][j][0]<dist[i][j][0]:
                dist[j][k]=[dist[i][N+p][0]+dist[N+p][j][0],dist[i][N+p][1]+dist[N+p][j][1]]

t1=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)]
ser=randint(0,N-1)
star[N]=[star[ser][0],star[ser][1],1]
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]]
temp=10**2
c=0.99
while time.time()-t1<0.8:
    x=randint(-10,10)
    y=randint(-10,10)
    p=randint(0,M-1)
    if 0<=star[N+p][0]+x<=1000 and 0<=star[N+p][1]+y<=1000:
        pre=sell(False)
        change(x,y,p)
        now=sell(False)
        if (pre-now)<0 and random()>pow(2.72,(pre-now)/temp):
            change(-x,-y,p)
    temp=temp*c
for i in range(M):
    print(*star[N+i][:2])
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)