题目描述请移步 https://www.luogu.com.cn/problem/P10936
- 题目简述
有n个防御塔，每个防御塔都有充足的导弹
导弹需要一定时间发出，又需要一定时间冷却
导弹有确定的速度，发出后会沿最短路径攻击任意一个入侵者
有m个入侵者，给定防御塔和入侵者的坐标，求至少多久才能击退所有入侵者
- 解题思路
经过分析，发现直接求最少多久似乎不太可行，于是考虑二分，将问题转化为验证一个时长是否可行
在一定时间内，每个防御塔可以发射固定数量个导弹，这些导弹必须可以击中所有的入侵者，否则不可行
我们可以处理出每个入侵者能被哪些导弹击中（第i个防御塔第1枚，第i个防御塔第2枚...），他必须被其中一枚击中
这有些类似二分图，将入侵者看作一类点，防御塔看作一类点，入侵者与能击中其的导弹连边
如果最大匹配小于入侵者数量，则不可行；反之，则可行
这样的复杂度或许有些高，但"N,M≤50"的数据还是可以通过的
- 易错点
-
题目中t1的单位是秒，需先转成分钟
-
本题许多变量涉及小数，如t1，二分l_r_mid，还有check函数上传的量，注意使用double
-
先输入的是m个入侵者的坐标，而非n个防御塔
-
代码实现
#include <iostream>
#include <utility>
#include <vector>
#include <cmath>
#include <cstring>
#include <iomanip>
using namespace std;
const int MAX=55;
const int MAXK=2505;
const double eps=1e-7;//二分精度
int n,m;
double t1,t2,v;
pair <double,double> a[MAX],b[MAX];
vector <int> edge[MAX];//哪些防御塔能打到这个敌人
int match[MAXK];//二分图匹配
bool vis[MAXK];
double dis(int i,int j){//欧几里得距离
return sqrt((a[i].first-b[j].first)(a[i].first-b[j].first)+(a[i].second-b[j].second)(a[i].second-b[j].second));
}
bool find(int x){//二分图
for(int j=0;j<edge[x].size();j++){
int now=edge[x][j];
if(vis[now]==false){
vis[now]=true;
if(match[now]==0||find(match[now])==true){
match[now]=x;
return true;
}
}
}
return false;
}
bool check(double x){
int tot=(x+t2)/(t1+t2);
tot=min(tot,m);//总共能发射多少枚导弹
for(int j=1;j<=m;j++){//敌人
edge[j].clear();
for(int i=1;i<=n;i++){//防御塔
for(int k=1;k<=tot;k++){
if(dis(i,j)/v+k*t1+(k-1)*t2<=x){//可以打到
edge[j].push_back((i-1)*tot+k);
}
}
}
}
memset(match,0,sizeof(match));
for(int i=1;i<=m;i++){
memset(vis,0,sizeof(vis));
if(find(i)==false){//如果有一个入侵者无法被击退则不可行
return false;
}
}
return true;
}
int main(){
cin>>n>>m>>t1>>t2>>v;
t1=t1/60.0;//注意t1单位是秒
for(int i=1;i<=m;i++){
cin>>b[i].first>>b[i].second;//入侵者
}
for(int i=1;i<=n;i++){
cin>>a[i].first>>a[i].second;//防御塔
}
double l=0,r=1e9;
while(r-l>eps){//小数二分
double mid=(l+r)/2.0;
if(check(mid)==true){
r=mid;
}else{
l=mid;
}
}
cout<<fixed<<setprecision(6)<<l;//依据题意保留小数
}