Hlg 1481 二分图匹配+二分.cpp

题意：

一个船队 VS 一只章鱼

给出 n 个定点<章鱼头+章鱼手> 再给出 p 个可移动点<船长+船员>..

对于 n 个定点 给出坐标 x y

对于 p 个可移动点 给出坐标 x y 和速度 v

先输入 T 组样例..

对于定点 n 第一个给出章鱼头坐标 接下来n-1行是每只章鱼手的坐标

对于移动点 p 第一个给出船长的信息 接下来p-1行是每个船员的信息

思路：

二分图匹配..

所谓最快时间内让船长击打到章鱼头 就是找出某个船员到章鱼手的最长时间

所以二分查找每个船员到对应章鱼手的花费时间..

如果以这个时间作为最长时间..那么就把比船员-章鱼手花费它时间少的连线..求最大匹配..匹配结果==章鱼手数的时候..

就继续二分..确保精确值..

答案就是这个时间+船长到章鱼头花费的时间..

Tips：

类型一定要明确好..已经不是第一次出现类型隐式转换导致数值改变的问题了.

这样的二分可以确定精确值..

二分的时候..什么时候左移..什么时候右移..

Code：

#include <stdio.h>
#include <cstring>
#include <cmath>
using namespace std;
#define clr(x) memset(x, 0, sizeof(x))
const int INF = 0x1f1f1f1f;

struct Node
{
    double x;
    double y;
    double v;
}pir[110], ten[110], ca, he;

struct Edge
{
    int next;
    int to;
}edge[1000000];
int tot;
int head[110];

void add(int s, int u)
{
    edge[tot].to = u;
    edge[tot].next = head[s];
    head[s] = tot++;
}

int link[110], vis[110];
double v[110][110];
int sum;
int n, p;
double l, r, mid;

int dfs(int x)
{
    for(int i = head[x]; i != 0; i = edge[i].next) {
        int y = edge[i].to;
        if(!vis[y]) {
            vis[y] = true;
            if(link[y] == 0 || dfs(link[y])) {
                link[y] = x;
                return true;
            }
        }
    }
    return false;
}

void solve()
{
    clr(link);
    sum = 0;
    for(int i = 1; i <= n; ++i) {
        clr(vis);
        if(dfs(i))
            sum++;
    }
}

double dis(Node a, Node b)
{
    return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
}

bool check(double t)
{
    tot = 1;
    clr(head);
    for(int i = 1; i <= n; ++i)
        for(int j = 1; j <= p; ++j) {
            if(v[i][j] < t)// ?<=
                add(i, j);
    }
    solve();
    if(sum == n)
        return true;
    else return false;
}

double solvee()
{
    l = 0;
    while(l < r) {
        mid = (l+r)/2;
        if(check(mid)) r = mid-0.0000000001;
        else l = mid+0.0000000001;
    }
    return l;
}

int main()
{
    int i, j, k;
    int T;
    //freopen("E:\\ACM\\Mess\\stdin-stdout\\in.txt", "r", stdin);
    scanf("%d", &T);
    while(T--)
    {
        clr(v);
        r = 0;

        scanf("%d %d", &n, &p);
        scanf("%lf %lf %lf", &ca.x, &ca.y, &ca.v);
        for(i = 1; i <= p; ++i)
            scanf("%lf %lf %lf", &pir[i].x, &pir[i].y, &pir[i].v);
        scanf("%lf %lf", &he.x, &he.y);
        for(i = 1; i <= n; ++i) {
            scanf("%lf %lf", &ten[i].x, &ten[i].y);
        }

        for(i = 1; i <= n; ++i)
        for(j = 1; j <= p; ++j) {
            v[i][j] = dis(pir[j], ten[i])/pir[j].v;
            if(v[i][j] > r) r = v[i][j];
        }
        double ans = solvee()+dis(ca, he)/ca.v;
        printf("%.9lf\n", ans);
    }
    return 0;
}