题意很简单，显然是对客人与伞做匹配，求最大匹配数。
一开始用了匈牙利，果断TLE。。。
好吧，原来二分匹配还有这个奇葩算法。。然而敲完还是不明白这个算法啥意思。。。
附两个算法的时间复杂度：
匈牙利：O(VE)
H-K ：O(V^0.5 E)
1 #include<cstdio> 2 #include<cstring> 3 #include<vector> 4 #include<queue> 5 using namespace std; 6 int nx,ny,mx[3300],my[3300],dx[3300],dy[3300],vis[3300]; 7 vector<int> v[3300]; 8 struct node1 9 { 10 int x,y,sp; 11 }per[3300]; 12 struct node2 13 { 14 int x,y; 15 }umb[3300]; 16 bool dfs(int u) 17 { 18 int v2,i; 19 for(i=0;i<v[u].size();i++) 20 { 21 v2=v[u][i]; 22 if(!vis[v2]&&dy[v2]==dx[u]+1) 23 { 24 vis[v2]=1; 25 if(my[v2]==-1||dfs(my[v2])) 26 { 27 my[v2]=u; 28 mx[u]=v2; 29 return true; 30 } 31 } 32 } 33 return false; 34 } 35 int hkmatch() 36 { 37 int ret=0; 38 memset(mx,-1,sizeof(mx)); 39 memset(my,-1,sizeof(my)); 40 while(true) 41 { 42 int i; 43 bool flag=false; 44 queue<int> q; 45 memset(dx,0,sizeof(dx)); 46 memset(dy,0,sizeof(dy)); 47 for(i=1;i<=nx;i++) 48 { 49 if(mx[i]==-1) 50 { 51 q.push(i); 52 } 53 } 54 while(!q.empty()) 55 { 56 int z=q.front(); 57 q.pop(); 58 for(i=0;i<v[z].size();i++) 59 { 60 int vv=v[z][i]; 61 if(!dy[vv]) 62 { 63 dy[vv]=dx[z]+1; 64 if(my[vv]==-1) flag=true; 65 else 66 { 67 dx[my[vv]]=dy[vv]+1; 68 q.push(my[vv]); 69 } 70 } 71 } 72 } 73 if(flag==false) break; 74 memset(vis,0,sizeof(vis)); 75 for(i=1;i<=nx;i++) 76 { 77 if(mx[i]==-1&&dfs(i)) 78 ret++; 79 } 80 } 81 return ret; 82 } 83 84 int main() 85 { 86 int cas,i,j,numcas=0; 87 scanf("%d",&cas); 88 while(cas--) 89 { 90 int t; 91 scanf("%d",&t); 92 scanf("%d",&nx); 93 for(i=1;i<=nx;i++) 94 v[i].clear(); 95 for(i=1;i<=nx;i++) 96 scanf("%d%d%d",&per[i].x,&per[i].y,&per[i].sp); 97 scanf("%d",&ny); 98 for(i=1;i<=ny;i++) 99 scanf("%d%d",&umb[i].x,&umb[i].y); 100 for(i=1;i<=nx;i++) 101 { 102 for(j=1;j<=ny;j++) 103 { 104 if(per[i].sp*t*per[i].sp*t>=(per[i].x-umb[j].x)*(per[i].x-umb[j].x)+(per[i].y-umb[j].y)*(per[i].y-umb[j].y)) 105 v[i].push_back(j); 106 } 107 } 108 printf("Scenario #%d: ",++numcas); 109 printf("%d ",hkmatch()); 110 } 111 return 0; 112 }