[luogu P2526][SHOI2001]小狗散步（二分图最大匹配）

Description

---

P2526 [SHOI2001]小狗散步
Grant喜欢带着他的小狗Pandog散步。Grant以一定的速度沿着固定路线走，该路线可能自交。Pandog喜欢游览沿途的景点，不过会在给定的N个点和主人相遇。小狗和主人同时从(X1，Y1)点出发，并同时在(Xn，Yn)点汇合。小狗的速度最快是Grant的两倍。当主人从一个点以直线走向另一个点时，Pandog跑向一个它感兴趣的景点。Pandog每次与主人相遇之前最多只去一个景点。

你现在的任务是：为Pandog寻找一条路线(有可能与主人的路线部分相同)，使它能够游览最多的景点，并能够准时与主人在给定地点相遇或者汇合。

Input

---

输入文件第一行是两个整数N和M( 1≤N，M≤100 )；
输入文件第二行的N个坐标给出了Grant的散步路线，即Pandog和主人相遇地点；
输入文件第三行的M个坐标给出了所有Pandog感兴趣的景点。
所有输入的坐标均不相同，且绝对值不超过1000。

Output

---

输出小狗的移动路线。
第一行是经过的点数，第二行依次为经过的点的坐标(直角坐标系)

Sample Input

---

4 5
1 4 5 7 5 2 -2 4
-4 -2 3 9 1 2 -1 3 8 -3

Sample Output

---

6
1 4 3 9 5 7 5 2 1 2 -2 4

Solution

---

Bravo!居然是洛谷上第一个A掉这道题的人…虽然本来也没几个人做
因为每次小狗与主人相遇前只能去一个景点，所以我们可以把这个情景看做二分图，尽可能地给每一段主人要走的路匹配一个景点

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<string>
#include<cmath>
#define eps 1e-6
using namespace std;
struct Point
{
    int x,y;
    Point(int x,int y):x(x),y(y){}
    Point():x(0),y(0){}
}p[105];
struct line
{
    Point p1,p2;
    line(Point p1,Point p2):p1(p1),p2(p2){}
    line():p1(Point()),p2(Point()){}
}l[105];
int n,m,visited[105],link[105],res=0;
double dis(Point u,Point v)
{
    return sqrt(double((u.x-v.x)*(u.x-v.x)+(u.y-v.y)*(u.y-v.y)));
}
bool check(Point P,line L)
{
    double d1=dis(P,L.p1);
    double d2=dis(P,L.p2);
    if(2*dis(L.p1,L.p2)-d1-d2>=-eps)return true;
    return false;
}
int dfs(int i)
{
    for(int j=1;j<n;j++)
    {
        if(check(p[i],l[j])&&!visited[j])
        {
            visited[j]=1; 
            if(!link[j]||dfs(link[j]))
            {
                link[j]=i;
                return true;
            }
        }
    }
    return false;
}
int main()
{
    scanf("%d%d",&n,&m);
    int x1,y1,x2,y2;
    Point p1,p2;
    scanf("%d%d",&x2,&y2);
    p2=Point(x2,y2);
    for(int i=1;i<n;i++)
    {
        scanf("%d%d",&x1,&y1);
        p1=Point(x1,y1);
        l[i]=line(p2,p1);
        p2=p1;
    }
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&x1,&y1);
        p[i]=Point(x1,y1);
    }
    for(int i=1;i<=m;i++)
    {
        memset(visited,0,sizeof(visited));
        if(dfs(i))res++;
    }
    printf("%d\n",res+n);
    printf("%d %d ",l[1].p1.x,l[1].p1.y);
    for(int i=1;i<n;i++)
    {
        if(link[i])
        {
            printf("%d %d ",p[link[i]].x,p[link[i]].y);
        }
        printf("%d %d ",l[i].p2.x,l[i].p2.y);
    }
    return 0;
}