POJ 1034 The dog task（二分图匹配）

http://poj.org/problem?id=1034

题意：

猎人和狗一起出去，狗的速度是猎人的两倍，给出猎人的路径坐标，除了这些坐标外，地图上还有一些有趣的点，而我们的狗，就是要尽量去多的有趣的点。前提是在猎人到达一个点的时候，狗必须也正好到。

思路：

二分图匹配，猎人的路径在一边，有趣的点在另一边，我们就是求一个最大匹配。

#include<iostream>
#include<algorithm>
#include<string>
#include<cstring>
#include<cmath>
using namespace std;

struct node
{
    int x, y;
};

node Bob[105], Dog[105];
int g[105][105];
int vis[105];
int match[105];
int path[105];
int n, m;


double dis(int x1, int y1, int x2, int y2)
{
    return sqrt((double)(x2 - x1)*(x2 - x1) + (y2 - y1)*(y2 - y1));
}


int dfs(int x)
{
    for (int i = 0; i < m; i++)
    {
        if (g[x][i] && !vis[i])
        {
            vis[i] = 1;
            if (match[i]==-1 || dfs(match[i]))
            {
                match[i] = x;
                return 1;
            }
        }
    }
    return 0;
}

int main()
{
    //freopen("D:\txt.txt", "r", stdin);
    while (cin >> n >> m)
    {
        for (int i = 0; i < n; i++)
            cin >> Bob[i].x >> Bob[i].y;
        for (int i = 0; i < m; i++)
            cin >> Dog[i].x >> Dog[i].y;

        memset(g, 0, sizeof(g));
        memset(match, -1, sizeof(match));

        for (int i = 0; i < n - 1; i++)
        for (int j = 0; j < m; j++)
        {
            double d1 = dis(Bob[i].x, Bob[i].y, Bob[i + 1].x, Bob[i + 1].y);
            double d2 = dis(Bob[i].x, Bob[i].y, Dog[j].x, Dog[j].y);
            double d3 = dis(Bob[i + 1].x, Bob[i + 1].y, Dog[j].x, Dog[j].y);
            if (d2 + d3 - 2 * d1 <= 0)
                g[i][j] = 1;
        }
        int ans = 0;
        for (int i = 0; i < n; i++)
        {
            memset(vis, 0, sizeof(vis));
            if(dfs(i)) ans++;
        }
        cout << ans + n << endl;

        memset(path, -1, sizeof(path));
        for (int i = 0; i < m; i++)
        {
            if (match[i] != -1)
                path[match[i]] = i;
        }
        for (int i = 0; i < n - 1; i++)
        {
            cout << Bob[i].x << " " << Bob[i].y << " ";
            if (path[i] != -1)
                cout << Dog[path[i]].x << " " << Dog[path[i]].y << " ";
        }
        cout << Bob[n - 1].x << " " << Bob[n - 1].y << endl;
    }
    return 0;
}