bzoj3035: 导弹防御塔

Description

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

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

Input

第一行五个正整数N,M,T1,T2,V。
接下来M行每行两个整数，代表入侵者的坐标。
接下来N行每行两个整数，代表防御塔的坐标。

Output

输出一个实数，表示最少需要多少分钟才能击中所有的入侵者，四舍五入保留六位小数。

二分答案，网络流判定

#include<cstdio>
#include<cmath>
const int N=1000000;
int n,m,S,T;
int mx[64],my[64],nx[64],ny[64];
int es[N],enx[N],e0[N],ep,h[N],q[N],ev[N];
double t1,t2,v;
bool bfs(){
    for(int i=1;i<=T;i++)h[i]=0;
    int ql=0,qr=0;
    h[q[qr++]=S]=1;
    while(ql!=qr){
        int w=q[ql++];
        for(int i=e0[w];i;i=enx[i])if(ev[i]){
            int u=es[i];
            if(!h[u])h[q[qr++]=u]=h[w]+1;
        }
    }
    return h[T];
}
int dfs(int w,int f){
    if(w==T)return f;
    int c,used=0;
    for(int i=e0[w];i;i=enx[i])if(ev[i]){
        int u=es[i];
        if(h[u]!=h[w]+1)continue;
        c=f-used;
        if(c>ev[i])c=ev[i];
        c=dfs(u,c);
        ev[i]-=c;
        ev[i^1]+=c;
        used+=c;
        if(f==used)return f;
    }
    h[w]=0;
    return used;
}
inline void adde(int a,int b){
    es[ep]=b;enx[ep]=e0[a];ev[ep]=1;e0[a]=ep++;
    es[ep]=a;enx[ep]=e0[b];ev[ep]=0;e0[b]=ep++;
}
bool chk(double t){
    ep=2;
    S=n*m+m+1;T=S+1;
    for(int i=1;i<=T;i++)e0[i]=0;
    for(int i=1;i<=n;i++){
        for(int k=0;k<m;k++){
            adde(S,(i-1)*m+k+1);
        }
    }
    for(int i=1;i<=m;i++)adde(n*m+i,T);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            double x=mx[j]-nx[i],y=my[j]-ny[i];
            double tm=sqrt(x*x+y*y)/v;
            if(tm>t)continue;
            for(int k=0;k<m;k++)if(tm+t1+(t1+t2)*k<=t){
                adde((i-1)*m+k+1,n*m+j);
            }
        }
    }
    int ans=0;
    while(bfs())ans+=dfs(S,0x3f3f3f3f);
    return ans==m;
}
int main(){
    scanf("%d%d%lf%lf%lf",&n,&m,&t1,&t2,&v);
    for(int i=1;i<=m;i++)scanf("%d%d",mx+i,my+i);
    for(int i=1;i<=n;i++)scanf("%d%d",nx+i,ny+i);
    t1/=60.;
    double L=0,R=1e5;
    for(int i=0;i<44;i++){
        double M=(L+R)*.5;
        if(chk(M))R=M;
        else L=M;
    }
    printf("%.6f
",L);
    return 0;
}