[P2526][SHOI2001]小狗散步

#include <bits/stdc++.h>

using namespace std;
typedef pair<int,int> P;
#define X first
#define Y second
const int MAXN=105;
P a[MAXN],b[MAXN];
vector<int> G[MAXN];
int n,m,mat[MAXN],vis[MAXN],idx=1;

int dfs(int x)
{
    vis[x]=idx;
    for(int i=0;i<G[x].size();i++)
    {
        int v=G[x][i],m=mat[v];
        if(m==-1||vis[m]!=idx&&dfs(m))
        {mat[v]=x;return 1;}
    }
    return 0;
}

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,-1,sizeof(mat));
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&a[i].X,&a[i].Y);
    for(int i=1;i<=m;i++)
        scanf("%d%d",&b[i].X,&b[i].Y);
    for(int i=2;i<=n;i++)
    {
        double d=dist(a[i],a[i-1]);
        for(int j=1;j<=m;j++)
            if(2*d>=dist(a[i],b[j])+dist(a[i-1],b[j])) 
                G[j].push_back(i);
    }    
    
    int res=0;
    for(int i=1;i<=m;i++,idx++) res+=dfs(i);
    printf("%d\n",res+n);    
    for(int i=1;i<=n;i++)
    {
        if(mat[i]!=-1) 
            printf("%d %d ",b[mat[i]].X,b[mat[i]].Y);
        printf("%d %d ",a[i].X,a[i].Y);
    }
    return 0;
}