题目描述:
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;
}
我的思路到拆点就结束了,看了这篇博客才知道需要二分