[二分+匈牙利]zoj 3460:Missile
大致题意:
? ? 用n个导弹发射塔攻击m个目标。每个发射架在某个时刻只能为一颗导弹服务,发射一颗导弹需要准备t1的时间,一颗导弹从发射到击中目标的时间与目标到发射架的距离有关。每颗导弹发射完成之后发射架需要t2的时间进入下个发射流程。现在问最少需要多少时间可以击毁所有m个目标。
?
大致思路:
? ? 二分枚举这个最大时间的最小值,每次按照这个枚举的时间构出二分图,求最大匹配来判定枚举值是否符合要求。
?
#include<iostream>#include<cstring>#include<cmath>#include<cstdio>using namespace std;const int inf=1<<30;const double eps=1e-8;int N,M;double tt1,tt2,v;int mis[100][2],tar[100][2];double dis[3000][60];double getdis(int a,int b){ double tmp=(mis[a][0]-tar[b][0])*(mis[a][0]-tar[b][0])+(mis[a][1]-tar[b][1])*(mis[a][1]-tar[b][1]); double res=sqrt(tmp); return res;}void getinit(){ int i,j,k; for(i=1;i<=N;i++) // 枚举第i个发射塔 { for(j=1;j<=M;j++) //第i个发射塔的第j次发射 { for(k=1;k<=M;k++) //第j次发射,攻击第k个目标 { dis[(i-1)*M+j][k]=(j-1)*tt2+getdis(i,k)/v+j*tt1; } } }}int map[3000][60];bool vis[3000];int linkk[3000];int dfs(int s){ for(int i=1;i<=M;i++){ if(!vis[i]&&map[s][i]){ vis[i]=1; if(linkk[i]==-1||dfs(linkk[i])){ linkk[i]=s; return 1; } } } return 0;}bool check(double mid){// cout<<"check"<<mid<<endl; int i,j,res=0; memset(map,0,sizeof(map)); for(i=1;i<=N*M;i++) { for(j=1;j<=M;j++) { if(dis[i][j]<=mid) { map[i][j]=1; } } } memset(linkk,-1,sizeof(linkk)); for(i=1;i<=N*M;i++){ memset(vis,0,sizeof(vis)); if(dfs(i))res++; } if(res==M)return 1; return 0;}int main(){ int i; while(scanf("%d%d%lf%lf%lf",&N,&M,&tt1,&tt2,&v)!=EOF) { tt1/=60.0; for(i=1;i<=M;i++) { scanf("%d%d",&tar[i][0],&tar[i][1]); } for(i=1;i<=N;i++) { scanf("%d%d",&mis[i][0],&mis[i][1]); } getinit(); double left=0,right=200000000,mid; while(right-left>=eps) { mid=(left+right)/2.0; if(check(mid))right=mid; else left= mid; } printf("%.6f\n",left); } return 0;}?