ZOJ3460 Missile（神奇的建图）

http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3460

感叹出题人的脑洞。。。

这里把第i个基地第j次发射导弹这一事件记作点(j-1)*N+i，这样就成功处理了分次发射导弹的问题。。每个点都和目标相连，边权为所花的时间，即：j*T1+（j-1）*T2+d/V，然后二分时间限制。求二分匹配时如果对应边的权大于时间限制就不选此边，然后看最大匹配数是否等于m，若否则增大时间限；否则减小时间限。直到这一次与上一次的时间限小于1e-6，既达到目标精度。。

代码：

#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
using namespace std;
#define N 55
#define M 2510
#define eps 1e-8
#define INF 200000000
#define min(x,y) ((x)<(y)?(x):(y))
#define max(x,y) ((x)>(y)?(x):(y))
#define d(x,y) sqrt((x)*(x)+(y)*(y))
double G[N][M],limit,A[N][2],B[M][2];
bool vis[M];
int Nx,Ny,Mx[N],My[M],dx[N],dy[M],dis;
bool bfs(){
    queue<int> q;dis=INF;int u;
    memset(dx,-1,sizeof(dx));
    memset(dy,-1,sizeof(dy));
    for(int i=1;i<=Nx;++i){
        if(Mx[i]==-1){q.push(i);dx[i]=0;}
    }
    while(!q.empty()){
        u=q.front();q.pop();
        if(dx[u]>dis) break;
        for(int i=1;i<=Ny;++i){
            if(G[u][i]<=limit&&dy[i]==-1){
                dy[i]=dx[u]+1;
                if(My[i]==-1) dis=dy[i];
                else{
                    dx[My[i]]=dy[i]+1;q.push(My[i]);
                }
            }
        }
    }
    return dis!=INF;
}
bool dfs(int u){
    for(int i=1;i<=Ny;++i){
        if(G[u][i]<=limit&&!vis[i]&&dy[i]==dx[u]+1){
            vis[i]=true;
            if(My[i]!=-1&&dy[i]==dis) continue;
            if(My[i]==-1||dfs(My[i])){
                My[i]=u;Mx[u]=i;return true;
            }
        }
    }
    return false;
}
int MaxMatch(){
    int ans=0;
    memset(Mx,-1,sizeof(Mx));
    memset(My,-1,sizeof(My));
    while(bfs()){
        memset(vis,false,sizeof(vis));
        for(int i=1;i<=Nx;++i){
            if(Mx[i]==-1&&dfs(i)) ++ans;
        }
    }
    return ans;
}
int main(){
    int n,m,ans;double L,R,T1,T2,v;
    while(scanf("%d%d%lf%lf%lf",&n,&m,&T1,&T2,&v)!=EOF){
        L=INF;R=-INF;Nx=m;Ny=n*m;T1=T1/60;
        for(int i=1;i<=m;++i) scanf("%lf%lf",&A[i][0],&A[i][1]);
        for(int i=1;i<=n;++i){
            scanf("%lf%lf",&B[i][0],&B[i][1]);
            for(int j=1;j<=m;++j){
                G[j][i]=d(A[j][0]-B[i][0],A[j][1]-B[i][1])/v+T1;
                L=min(L,G[j][i]);R=max(R,G[j][i]);
                for(int k=1;k<m;++k){
                    G[j][k*n+i]=G[j][i]+(T1+T2)*k;
                    L=min(L,G[j][k*n+i]);R=max(R,G[j][k*n+i]);
                }
            }
        }
        while(R-L>eps){
            limit=(L+R)/2;
            ans=MaxMatch();
            if(ans==m) R=limit;
            else L=limit;
        }
        printf("%.6lf\n",limit);
    }
    return 0;
}