#include <bits/stdc++.h>

using namespace std;

typedef pair<int,int> P;

#define X first

#define Y second

const int MAXN=;

P a[MAXN],b[MAXN];

vector<int> G[MAXN];

int n,m,mat[MAXN],vis[MAXN],idx=;

int dfs(int x)

{

    vis[x]=idx;

    for(int i=;i<G[x].size();i++)

    {

        int v=G[x][i],m=mat[v];

        if(m==-||vis[m]!=idx&&dfs(m))

        {mat[v]=x;return ;}

    }

    return ;

}

double dist(P a,P b)

{return sqrt((a.X-b.X)*(a.X-b.X)+(a.Y-b.Y)*(a.Y-b.Y));}

int main()

{

    memset(mat,-,sizeof(mat));

    scanf("%d%d",&n,&m);

    for(int i=;i<=n;i++)

        scanf("%d%d",&a[i].X,&a[i].Y);

    for(int i=;i<=m;i++)

        scanf("%d%d",&b[i].X,&b[i].Y);

    for(int i=;i<=n;i++)

    {

        double d=dist(a[i],a[i-]);

        for(int j=;j<=m;j++)

            if(*d>=dist(a[i],b[j])+dist(a[i-],b[j]))

                G[j].push_back(i);

    }    

    int res=;

    for(int i=;i<=m;i++,idx++) res+=dfs(i);

    printf("%d\n",res+n);

    for(int i=;i<=n;i++)

    {

        if(mat[i]!=-)

            printf("%d %d ",b[mat[i]].X,b[mat[i]].Y);

        printf("%d %d ",a[i].X,a[i].Y);

    }

    return ;

}