导弹防御塔

题目描述：

Freda的城堡——
“Freda，城堡外发现了一些入侵者！”
“喵...刚刚探究完了城堡建设的方案数，我要歇一会儿嘛lala~”
“可是入侵者已经接近城堡了呀！”
“别担心，rainbow，你看呢，这是我刚设计的导弹防御系统的说~”

“喂...别卖萌啊……”

Freda控制着N座可以发射导弹的防御塔。每座塔都有足够数量的导弹，但是每座塔每次只能发射一枚。在发射导弹时，导弹需要T1秒才能从防御塔中射出，而在发射导弹后，发射这枚导弹的防御塔需要T2分钟来冷却。
所有导弹都有相同的匀速飞行速度V，并且会沿着距离最短的路径去打击目标。计算防御塔到目标的距离Distance时，你只需要计算水平距离，而忽略导弹飞行的高度。导弹在空中飞行的时间就是 (Distance/V) 分钟，导弹到达目标后可以立即将它击毁。
现在，给出N座导弹防御塔的坐标，M个入侵者的坐标，T1、T2和V，你需要求出至少要多少分钟才能击退所有的入侵者。
 

输入描述:

第一行五个正整数N,M,T1,T2,V。 接下来M行每行两个整数，代表入侵者的坐标。 接下来N行每行两个整数，代表防御塔的坐标。 

输出描述:

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

输入

3 3 30 20 1
0 0
0 50
50 0
50 50
0 1000
1000 0

输出

91.500000

备注:

对于40%的数据，N,M≤20

对于100%的数据，1≤N≤50,1≤M≤50，坐标绝对值不超过10000，T1,T2,V不超过2000.

方法：

做这道题的时候，因为开了知识点提示，所以知道用二分图。二分图需要满足右边最多只能被选择一次，而本题一座防御塔可以被选择多次，所以用拆点的方式，用多个点表示一座防御塔的多个导弹，可是普通二分图只能解决最大匹配问题，无法选择最优方案，所以这道题要再用二分来做，mid在此题的定义是能否在mid的时间内消灭入侵者，如果 mid 时间内，第 i 个入侵者能被第 j 座塔的第 k 个导弹击中，那么就在第 i 个左部节点和第 (k - 1) *n + j 个右部节点之间连一条无向边。

#include <iostream>
#include <cstring>
#include <vector>
#include <cmath>

#define x first
#define y second
 
using namespace std;

const int N=55,M=2600;
const double eps=1e-9;
typedef pair<int,int> pii;
pii a[N],b[N];
double dist[N][N];
bool st[M];
int match[M];
vector<int> adj[N];
int n,m;
double t1,t2,v;

double get_dist(int x1,int y1,int x2,int y2)
{
    double x=x1-x2,y=y1-y2;
    return sqrt(x*x+y*y);
}

bool dfs(int u)
{
    for(int i=0;i<adj[u].size();i++)
    {
        int j=adj[u][i];
        if(st[j])
            continue;
        st[j]=1;
        if(!match[j]||dfs(match[j]))
        {
            match[j]=u;
            return true;
        }
    }
    return false;
}

bool check(double mid)
{
    int p = (mid + t2) / (t1 + t2);
    p=min(m,p);
    for(int i=1;i<=m;i++)
    {
        adj[i].clear();
        for(int j=1;j<=n;j++)
        {
            for(int k=0;k<p;k++)
            {
                if(k*t2+(k+1)*t1+dist[i][j]/v < mid)//第 i 个入侵者能被第 j 座塔的第 k 个导弹击中
                    adj[i].push_back(n*k+j);
            }
        }
    }
    memset(match, 0, sizeof match);
    for(int i=1;i<=m;i++)
    {
        memset(st,0,sizeof st);
        if(!dfs(i)) return false;
    }
    return true;
}

int main()
{
    cin>>n>>m>>t1>>t2>>v;
    t1/=60;
    for(int i=1;i<=m;i++)
    {
        int x,y;
        cin>>x>>y;
        a[i]={x,y};
    }
    for(int i=1;i<=n;i++)
    {
        int x,y;
        cin>>x>>y;
        b[i]={x,y};
    }
    for(int i=1;i<=m;i++)//建图
    {
        for(int j=1;j<=n;j++)
        {
            dist[i][j]=get_dist(a[i].first,a[i].second,b[j].first,b[j].second);
        }
    }
    //二分
    double l=0,r=1e9;
    while(r-l>eps)
    {
        double mid=(l+r)/2;
        if(check(mid))r=mid;
        else l=mid;
    }
    printf("%.6f",l);
    return 0;
}

我的思路到拆点就结束了，看了这篇博客才知道需要二分

374. 导弹防御塔_满地都是六便士他却抬头看向了月亮的博客-CSDN博客