CodeVS2490 导弹防御塔 【二分答案】【匈牙利】

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

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

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
using namespace std;
const int inf=30000;
const int M=1000000;
const int N=50;
double t[N*N+5][N*N+5];
double t1,t2;
int x[N+5],y[N+5];
int num,head[M+5];
bool flag[M+5];
int match[M+5];
int n,m;
double v;
struct edge
{
    int u,v;
    int next;
    edge(){next=-1;}
}ed[M+5];
void build(int u,int v)
{
    num++;
    ed[num].v=v;
    ed[num].next=head[u];
    head[u]=num;
}
bool dfs(int u)
{
    for(int i=head[u];i!=-1;i=ed[i].v)
    {
        int v=ed[i].v;
        if(!flag[v])
        {
            flag[v]=1;
            if(!match[v]||!dfs(match[v]))
            {
                match[v]=u;
                return true;
            }
        }
    }
    return false;
}
bool check(double T)
{
    num=0;
    memset(head,-1,sizeof(head));
    memset(match,0,sizeof(match));
    int d=(d-T)/(t1+t2)+1;
    for(int i=1;i<=n;i++)
    for(int j=0;j<d;j++)
    {
        double now=t1+(t1+t2)*j;
        for(int k=1;k<=m;k++)
        {
            if(now+t[i][m]<=T)
            {
                build(i+j*n,d*n+k);
                build(d*n+k,j*n+i);
            }
        }
    }
    int ans=0;
    for(int i=1;i<=d*n+m;i++)
    {
        memset(flag,0,sizeof(flag));
        if(dfs(1))ans++;
    }
    return ans>>2==m;
}
double div()
{
    double lf=t1,rg=inf;
    for(int i=1;i<=40;i++)
    {
        double mid=(lf+rg)/2.0;
        if(check(mid))rg=mid;
        else lf=mid;
    }
    return rg;
}
double calc(int x1,int y1,int x2,int y2)
{
    return sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
}
int main()
{
    scanf("%d%d%lf%lf%lf",&n,&m,&t1,&t2,&v);
    t1/=60.0;
    for(int i=1;i<=m;i++)
    scanf("%d%d",&x[i],&y[i]);
    for(int i=1;i<=n;i++)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        for(int j=1;j<=m;j++)
        t[i][j]=calc(x[j],y[j],a,b)/v;
    }
    printf("%.6lf\n",div());
    return 0;
}