zoj 3460 二分+二分图匹配

  1 /*
  2 大致题意：
  3 
  4 用n个导弹发射塔攻击m个目标。每个发射架在某个时刻只能为
  5 一颗导弹服务，发射一颗导弹需要准备t1的时间，一颗导弹从发
  6 射到击中目标的时间与目标到发射架的距离有关。每颗导弹发
  7 射完成之后发射架需要t2的时间进入下个发射流程。现在问
  8 最少需要多少时间可以击毁所有m个目标。
  9 
 10 
 11 大致思路：
 12 二分枚举这个最大时间的最小值，每次按照这个枚举的时间构出
 13 二分图，求最大匹配来判定枚举值是否符合要求。
 14 
 15 注意单位，T1要除于60转化成分的
 16 
 17 */
 18 
 19 #include<stdio.h>
 20 #include<string.h>
 21 #include<iostream>
 22 #include<algorithm>
 23 #include<math.h>
 24 #include<vector>
 25 using namespace std;
 26 const double eps=1e-8;
 27 const int MAXN=2550;
 28 int linker[MAXN];
 29 bool used[MAXN];
 30 vector<int>g[60];
 31 int uN;
 32 bool dfs(int u)
 33 {
 34     for(int i=0;i<g[u].size();i++)
 35     {
 36         if(!used[g[u][i]])
 37         {
 38             used[g[u][i]]=true;
 39             if(linker[g[u][i]]==-1||dfs(linker[g[u][i]]))
 40             {
 41                 linker[g[u][i]]=u;
 42                 return true;
 43             }
 44         }
 45     }
 46     return false;
 47 }
 48 int hungary()
 49 {
 50     int u;
 51     int res=0;
 52     memset(linker,-1,sizeof(linker));
 53     for(u=0;u<uN;u++)
 54     {
 55         memset(used,false,sizeof(used));
 56         if(dfs(u))res++;
 57     }
 58     return res;
 59 }
 60 
 61 int N,M;
 62 double  T1,T2,V;
 63 struct Node
 64 {
 65     int x,y;
 66 };
 67 Node node1[60],node2[60];
 68 double d[60][60];
 69 double tt[MAXN][60];
 70 
 71 double dis(Node a,Node b)
 72 {
 73     return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
 74 }
 75 
 76 void init()
 77 {
 78     for(int i=0;i<N;i++)
 79       for(int j=0;j<M;j++)
 80         d[i][j]=dis(node1[i],node2[j]);
 81     for(int k=0;k<M;k++)
 82       for(int i=0;i<N;i++)
 83         for(int j=0;j<M;j++)
 84         {
 85             tt[i*M+k][j]=k*T2+(k+1)*T1+d[i][j]/V;
 86         }
 87     uN=M;
 88 }
 89 
 90 double solve()
 91 {
 92     double l=0;
 93     double r=200000000000.0;
 94     double mid;
 95     while(r-l>=eps)
 96     {
 97         mid=(l+r)/2;
 98         for(int i=0;i<M;i++)g[i].clear();
 99         for(int i=0;i<M*N;i++)
100           for(int j=0;j<M;j++)
101           {
102               if(tt[i][j]<=mid)g[j].push_back(i);
103           }
104         if(hungary()==M)
105         {
106            r=mid;
107         }
108         else l=mid;
109     }
110     printf("%.6lf\n",r);
111 }
112 int main()
113 {
114     //freopen("in.txt","r",stdin);
115     //freopen("out.txt","w",stdout);
116     while(scanf("%d%d%lf%lf%lf",&N,&M,&T1,&T2,&V)!=EOF)
117     {
118         T1/=60;//这个注意，没有除一直得不到答案，纠结
119         for(int i=0;i<M;i++)scanf("%d%d",&node2[i].x,&node2[i].y);
120         for(int i=0;i<N;i++)scanf("%d%d",&node1[i].x,&node1[i].y);
121         init();
122         solve();
123     }
124     return 0;
125 }