題目描述

Freda控制著N座可以發射導彈的防禦塔。每座塔都有足夠數量的導彈，但是每座塔每次只能發射一枚。在發射導彈時，導彈需要T1秒才能從防禦塔中射出，而在發射導彈後，發射這枚導彈的防禦塔需要T2分鐘來冷卻。
所有導彈都有相同的勻速飛行速度V，並且會沿著距離最短的路徑去打擊目標。計算防禦塔到目標的距離Distance時，你只需要計算水平距離，而忽略導彈飛行的高度。導彈在空中飛行的時間就是 (Distance/V) 分鐘，導彈到達目標後可以立即將它擊毀。
現在，給出N座導彈防禦塔的座標，M個入侵者的座標，T1、T2和V，你需要求出至少要多少分鐘才能擊退所有的入侵者。

輸入格式

第一行五個正整數N,M,T1,T2,V。
接下來M行每行兩個整數，代表入侵者的座標。
接下來N行每行兩個整數，代表防禦塔的座標。

輸出格式

輸出一個實數，表示最少需要多少分鐘才能擊中所有的入侵者，四捨五入保留六位小數。

提示

對於40%的資料，N,M<=20.
對於100%的資料， 1≤N≤50, 1≤M≤50，座標絕對值不超過10000，T1,T2,V不超過2000.

樣例資料

分析：

《lyd演算法進階》二分多重匹配例題。二分規定的時間內是否所有的入侵者能否被消滅即可，即是否存在最大二分匹配。

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<string>
#include<vector>
#include<queue>
#include<stack>
#include<cmath>
#include<set>
#include<map>
using namespace std;
#define rep(i,a,n) for (int i=a;i<n;i++)
#define per(i,a,n) for (int i=n-1;i>=a;i--)
#define oo cout<<"!!!"<<endl;
typedef long long ll;
typedef unsigned long long ull;
#define ms(s) memset(s, 0, sizeof(s))
const int inf = 0x3f3f3f3f;
//head

const int maxn = 1111;
std::vector<int> g[maxn];
double t1,t2,V;
double dis[maxn][maxn];
int n,m;

struct max_mactch
{
	vector<int>g[maxn];
	bool vis[maxn];
	int left[maxn];

	void init()
	{
		rep(i,0,m+11)g[i].clear();
		memset(left,-1,sizeof left);
	}

	bool match(int u)
	{
		rep(i,0,g[u].size())
		{
			int v = g[u][i];
			if(!vis[v])
			{
				vis[v] = true;
				if(left[v] == -1 || match(left[v]))
				{
					left[v] = u;
					return true;
				}
			}
		}
		return false;
	}

	bool C(double mid)
	{
		double T = mid;
		init();
		int k = 1;

		while(T>=t1)
		{
			T -= t1;
			k++;
			rep(i,1,m+1)
				rep(j,1,n+1)
				if(dis[j][i] <= T)
					g[i].push_back(m+(k-1)*n+j);
			T-=t2;
		}
		int ans = 0;
		rep(i,1,m+1)
		{
			ms(vis);
			if(match(i))
				ans++;
		}
		return ans == m;

	}
}MM;
struct node
{
	double x,y;
}a[maxn];

int main() 
{
	
	cin>>n>>m>>t1>>t2>>V;
	t1 /= 60;

	rep(i,1,m+1)
	{
		cin>>a[i].x>>a[i].y;
	}
	rep(i,1,n+1)
	{
		double x,y;
		cin>>x>>y;
		rep(j,1,m+1)
		{
			dis[i][j] = sqrt((x-a[j].x) * (x-a[j].x) + (y-a[j].y) * (y - a[j].y)) / V;
		}

	}
	double l = t1,r = 1111111;
	rep(i,1,111)
	{
		double mid = (l+r)/2;
		if(MM.C(mid))r = mid;
		else l = mid;
	}
	printf("%.6lf",r);
}