【NOIP2013模拟】导弹防御塔

题目

Freda的城堡——
“Freda，城堡外发现了一些入侵者！”
“喵...刚刚探究完了城堡建设的方案数，我要歇一会儿嘛lala~”
“可是入侵者已经接近城堡了呀！”
“别担心，rainbow，你看呢，这是我刚设计的导弹防御系统的说~”
“喂...别卖萌啊……”
Freda控制着N座可以发射导弹的防御塔。每座塔都有足够数量的导弹，但是每座塔每次只能发射一枚。在发射导弹时，导弹需要T1秒才能从防御塔中射出，而在发射导弹后，发射这枚导弹的防御塔需要T2分钟来冷却。
所有导弹都有相同的匀速飞行速度V，并且会沿着距离最短的路径去打击目标。计算防御塔到目标的距离Distance时，你只需要计算水平距离，而忽略导弹飞行的高度。导弹在空中飞行的时间就是 (Distance/V) 分钟，导弹到达目标后可以立即将它击毁。
现在，给出N座导弹防御塔的坐标，M个入侵者的坐标，T1、T2和V，你需要求出至少要多少分钟才能击退所有的入侵者。

分析

首先二分时间，接着，就可以吧每个防御塔分成g个点，g表示在二分出的时间内可以发射多少枚导弹，然后把每个分出来的点和每个入侵者比较一下，如果这个点在二分出的时间内可以咔擦掉这个入侵者，就将它们连一条边。
然后。。。
就最大匹配一下就可以啦！
最大匹配可以用网络流或者匈牙利，如果不嫌麻烦的话就用网络流吧。

#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const int maxlongint=2147483647;
using namespace std;
double t1,t2,v,a[60][60],lf[60][3],xhm[60][3];
int n,m,tot,next[800000],to[800000],last[1000000],xyl[800000];
bool used[800000];
double dis(double x,double y,double x1,double y1)
{
    return sqrt((x1-x)*(x1-x)+(y1-y)*(y1-y));
}
int bj(int x,int y)
{
    next[++tot]=last[x];
    last[x]=tot;
    to[tot]=y;
}
bool find(int x)
{
    int i,j,k,l;
    for(i=last[x];i;i=next[i])
    {
        j=to[i];
        if(!used[j])
        {
            used[j]=true;
            if(xyl[j]==0 || find(xyl[j]))
            {
                xyl[j]=x;
                return true;
            }
        }
    }
    return false;
}
bool ddx(double limit)
{
    int i,j,k,l;
    double cs=(limit-t1)/(t1+t2)+1;
    tot=0;
    memset(to,0,sizeof(to));
    memset(last,0,sizeof(last));
    memset(next,0,sizeof(next));
    for(i=1;i<=n;i++)
        for(j=0;j<int(cs);j++)
        {
            double time=j*(t1+t2)+t1;
            for(k=1;k<=m;k++)
            {
                if(time+a[i][k]<=limit)
                {
                    bj(int(n*j+i),int(cs*n+k));
                    bj(int(cs*n+k),int(n*j+i));
                }
            }
        }
    int ans=0;
    memset(xyl,0,sizeof(xyl));
    for(i=1;i<=int(n*cs+m);i++)
    {
        memset(used,0,sizeof(used));
        if(find(i)) ans++;
    }
    return ans==m*2;
}
double rf()
{
    double l=0,r=50000,mid;
    int i,j,k;
    for(i=1;i<=50;i++)
    {
        mid=(l+r)/2;
        if(ddx(mid)) r=mid; 
        else l=mid;
    }
    return r;
}
int main()
{
    scanf("%d%d%lf%lf%lf",&n,&m,&t1,&t2,&v);
    t1=t1/60;
    int i,j,k,l,x,y;
    for(i=1;i<=m;i++)
    {
        scanf("%lf%lf",&xhm[i][1],&xhm[i][2]);
    }
    for(i=1;i<=n;i++)
    {
        scanf("%lf%lf",&lf[i][1],&lf[i][2]);
    }
    for(i=1;i<=n;i++)
        for(j=1;j<=m;j++)
        {
            a[i][j]=dis(lf[i][1],lf[i][2],xhm[j][1],xhm[j][2])/v;
        }
    double ans;
    ans=rf();
    printf("%.6lf",ans);
}