Hlg 1481 二分图匹配+二分.cpp

题意：

一个船队 VS 一只章鱼

给出 n 个定点<章鱼头+章鱼手> 再给出 p 个可移动点<船长+船员>..

对于 n 个定点 给出坐标 x y

对于 p 个可移动点 给出坐标 x y 和速度 v

先输入 T 组样例..

对于定点 n 第一个给出章鱼头坐标 接下来n-1行是每只章鱼手的坐标

对于移动点 p 第一个给出船长的信息 接下来p-1行是每个船员的信息

思路：

二分图匹配..

所谓最快时间内让船长击打到章鱼头 就是找出某个船员到章鱼手的最长时间

所以二分查找每个船员到对应章鱼手的花费时间..

如果以这个时间作为最长时间..那么就把比船员-章鱼手花费它时间少的连线..求最大匹配..匹配结果==章鱼手数的时候..

就继续二分..确保精确值..

答案就是这个时间+船长到章鱼头花费的时间..

Tips：

类型一定要明确好..已经不是第一次出现类型隐式转换导致数值改变的问题了.

这样的二分可以确定精确值..

二分的时候..什么时候左移..什么时候右移..

Code：

View Code 
  1 #include <stdio.h>
  2 #include <cstring>
  3 #include <cmath>
  4 using namespace std;
  5 #define clr(x) memset(x, 0, sizeof(x))
  6 const int INF = 0x1f1f1f1f;
  7 
  8 struct Node
  9 {
 10     double x;
 11     double y;
 12     double v;
 13 }pir[110], ten[110], ca, he;
 14 
 15 struct Edge
 16 {
 17     int next;
 18     int to;
 19 }edge[1000000];
 20 int tot;
 21 int head[110];
 22 
 23 void add(int s, int u)
 24 {
 25     edge[tot].to = u;
 26     edge[tot].next = head[s];
 27     head[s] = tot++;
 28 }
 29 
 30 int link[110], vis[110];
 31 double v[110][110];
 32 int sum;
 33 int n, p;
 34 double l, r, mid;
 35 
 36 int dfs(int x)
 37 {
 38     for(int i = head[x]; i != 0; i = edge[i].next) {
 39         int y = edge[i].to;
 40         if(!vis[y]) {
 41             vis[y] = true;
 42             if(link[y] == 0 || dfs(link[y])) {
 43                 link[y] = x;
 44                 return true;
 45             }
 46         }
 47     }
 48     return false;
 49 }
 50 
 51 void solve()
 52 {
 53     clr(link);
 54     sum = 0;
 55     for(int i = 1; i <= n; ++i) {
 56         clr(vis);
 57         if(dfs(i))
 58             sum++;
 59     }
 60 }
 61 
 62 double dis(Node a, Node b)
 63 {
 64     return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));
 65 }
 66 
 67 bool check(double t)
 68 {
 69     tot = 1;
 70     clr(head);
 71     for(int i = 1; i <= n; ++i)
 72         for(int j = 1; j <= p; ++j) {
 73             if(v[i][j] < t)// ?<=
 74                 add(i, j);
 75     }
 76     solve();
 77     if(sum == n)
 78         return true;
 79     else return false;
 80 }
 81 
 82 double solvee()
 83 {
 84     l = 0;
 85     while(l < r) {
 86         mid = (l+r)/2;
 87         if(check(mid)) r = mid-0.0000000001;
 88         else l = mid+0.0000000001;
 89     }
 90     return l;
 91 }
 92 
 93 int main()
 94 {
 95     int i, j, k;
 96     int T;
 97     //freopen("E:\\ACM\\Mess\\stdin-stdout\\in.txt", "r", stdin);
 98     scanf("%d", &T);
 99     while(T--)
100     {
101         clr(v);
102         r = 0;
103 
104         scanf("%d %d", &n, &p);
105         scanf("%lf %lf %lf", &ca.x, &ca.y, &ca.v);
106         for(i = 1; i <= p; ++i)
107             scanf("%lf %lf %lf", &pir[i].x, &pir[i].y, &pir[i].v);
108         scanf("%lf %lf", &he.x, &he.y);
109         for(i = 1; i <= n; ++i) {
110             scanf("%lf %lf", &ten[i].x, &ten[i].y);
111         }
112 
113         for(i = 1; i <= n; ++i)
114         for(j = 1; j <= p; ++j) {
115             v[i][j] = dis(pir[j], ten[i])/pir[j].v;
116             if(v[i][j] > r) r = v[i][j];
117         }
118         double ans = solvee()+dis(ca, he)/ca.v;
119         printf("%.9lf\n", ans);
120     }
121     return 0;
122 }

 

 

 

posted @ 2012-09-14 12:29  Griselda.  阅读(196)  评论(0)  编辑  收藏  举报