【ACWing】374. 导弹防御塔（二分图多重匹配+二分）

Freda的城堡遭受了 M 个入侵者的攻击！

Freda控制着 N 座导弹防御塔，每座塔都有足够数量的导弹，但是每次只能发射一枚。

在发射导弹时，导弹需要 T1T1 秒才能从防御塔中射出，而在发射导弹后，发射这枚导弹的防御塔需要 T2T2 分钟来冷却。

所有导弹都有相同的匀速飞行速度 V，并且会沿着距离最短的路径去打击目标。

计算防御塔到目标的距离Distance时，你只需要计算水平距离，而忽略导弹飞行的高度。

导弹在空中飞行的时间就是 (Distance/V) 分钟，导弹到达目标后可以立即将它击毁。

现在，给出 N 座导弹防御塔的坐标，M 个入侵者的坐标，T1,T2T1,T2 和 V。因为Freda的小伙伴Rainbow就要来拜访城堡了，你需要求出至少多少分钟才能击退所有的入侵者。

输入格式

第一行五个正整数N,M,T1T1,T2T2,V。

接下来 M 行每行两个整数，代表入侵者的坐标。

接下来 N 行每行两个整数，代表防御塔的坐标。

输出格式

输出一个实数，表示最少需要多少分钟才能击中所有的入侵者，四舍五入保留六位小数。

思路：这个题目的话，明显是个多重的二分图多重匹配，因为明显这个题目的导弹是一组，防御塔是一组，并且是没有任何交集的，但是我们的导弹是可以多次发射的，这是很关键的一点，所以我们选择多重匹配解决方案中的拆点

还有就是这个题目的答案，也就是每一次我们选择某个导弹发射多次的时候，答案是满足单调性的，这时候可以选择二分，二分答案，这样既可以确定出有一个防御塔有多少个导弹去打，

然后就是按照防御塔建图，检出图来就可以跑一遍匈牙利算法，看看是不是每一个塔都有导弹去打。

代码：

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<vector>
#include<cmath>
using namespace std;
const int maxn=50+10;
struct Node
{
    int x,y;
    Node() {}
    Node(int a,int b):x(a),y(b) {}
} def[maxn],ta[maxn];
int vis[maxn*maxn],pre[maxn*maxn];
double times[maxn][maxn];
double t1,t2,V;
vector<int >G[maxn];
int n,m;
bool DFS(int u)
{
    for(int i=0; i<G[u].size(); i++)
    {
        int v=G[u][i];
        if(vis[v])
            continue;
        vis[v]=1;
        if(pre[v]==0||DFS(pre[v]))
        {
            pre[v]=u;
            return true;
        }
    }
    return false;
}
bool judge(double mid)
{
    for(int i=1; i<=m; i++)
        G[i].clear();
    for(int i=1; i<=m; i++)
        for(int j=1; j<=n; j++)
        {
            double now=0;
            for(int k=0; k<=m; k++)
            {
                now=k*(t1+t2);
                if(mid-(now+t1+times[i][j])>1e-7)
                {
                    G[i].push_back(j*(m+1)+k-m-1);
                }
                else
                    break;
            }
        }
    memset(pre,0,sizeof(pre));
    for(int i=1; i<=m; i++)
    {
        memset(vis,0,sizeof(vis));
        if(!DFS(i))
            return false;
    }
    return true;
}
int main()
{
    cin>>n>>m>>t1>>t2>>V;
    t1/=60.0;
    for(int i=1; i<=m; i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        ta[i]=Node(x,y);
    }
    for(int i=1; i<=n; i++)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        def[i]=Node(x,y);
    }
    for(int i=1; i<=m; i++)
        for(int j=1; j<=n; j++)
            times[i][j]=sqrt((ta[i].x-def[j].x)*(ta[i].x-def[j].x)+(ta[i].y-def[j].y)*(ta[i].y-def[j].y))/V;
    /*for(int i=1;i<=m;i++)
    {
        for(int j=1;j<=n;j++)
            cout<<times[i][j]<<" ";
        cout<<endl;
    }*/
    double l=0,r=1000000,ans=r;
    while(r-l>1e-9)
    {
        double mid=(l+r)/2.0;
        if(judge(mid))
            ans=min(ans,mid),r=mid;
        else
            l=mid;
    }
    printf("%.6lf\n",ans);
}