# HDU 2389 二分图最大匹配之Hopcroft-Karp优化 O(n^0.5*m)

T个测试用例， 时间time, n个人，接下来n行是坐标 和 每个人的速度 v。 m,接下来m行是 雨伞的坐标。

#include <bits/stdc++.h>
using namespace std;
const int MAXN = 3010;
const int MAXM = 3010*3010;
const int INF = 0x3f3f3f3f;
struct Edge {
int v;
int next;
} edge[MAXM];
struct node {
double x, y;
double v;
} a[MAXN], b[MAXN];
int nx, ny;
int cnt;
int t;
int dis;
int first[MAXN];
int dx[MAXN], dy[MAXN]; /*dx[i]表示左集合i顶点的距离编号，dy[i]表示右集合i顶点的距离编号*/
int vis[MAXN]; //寻找增广路的标记数组
void init() {
cnt = 0;
memset(first, -1, sizeof(first));
}
void read_graph(int u, int v) {
edge[cnt].v = v;
edge[cnt].next = first[u], first[u] = cnt++;
}
int bfs() {
queue<int> q;
dis = INF;
memset(dx, -1, sizeof(dx));
memset(dy, -1, sizeof(dy));
for(int i = 0; i < nx; i++)	{
q.push(i);
dx[i] = 0;
}
}
while(!q.empty())	{
int u = q.front();
q.pop();
if(dx[u] > dis) break;
for(int e = first[u]; e != -1; e = edge[e].next)		{
int v = edge[e].v;
if(dy[v] == -1)			{
dy[v] = dx[u] + 1;
if(ylink[v] == -1) dis = dy[v];
else				{
}
}
}
}
return dis != INF;
}
int find(int u) {
for(int e = first[u]; e != -1; e = edge[e].next){
int v = edge[e].v;
if(!vis[v] && dy[v] == dx[u]+1)	{
vis[v] = 1;
if(ylink[v] != -1 && dy[v] == dis) continue;
return 1;
}
}
}
return 0;
}
int MaxMatch() {
int ans = 0;
while(bfs())	{
memset(vis, 0, sizeof(vis));
for(int i = 0; i < nx; i++) if(xlink[i] == -1)	{
ans += find(i);
}
}
return ans;
}
double dist(const node a, const node b) {
return sqrt((a.x-b.x)*(a.x-b.x) + (a.y-b.y)*(a.y-b.y));
}
init();
int Time;
scanf("%d", &Time);
scanf("%d", &nx);
for(int i = 0; i < nx; i++)	{
scanf("%lf%lf%lf", &a[i].x, &a[i].y, &a[i].v);
}
scanf("%d", &ny);
for(int i = 0; i < ny; i++)	{
scanf("%lf%lf", &b[i].x, &b[i].y);
}
for(int i = 0; i < nx; i++)	{
for(int j = 0; j < ny; j++)		{
double limit = a[i].v*Time;
double s = dist(a[i], b[j]);
}
}
}
void solve() {
int ans = MaxMatch();
printf("%d\n\n", ans);
}
int main() {
int T, times = 0;
scanf("%d", &T);
while(T--)	{
printf("Scenario #%d:\n", ++times);
solve();
}
return 0;
}

0 条评论

• ### ICPC Asia Shenyang 2019 Dudu's maze

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

• ### Educational Codeforces Round 67 (Rated for Div. 2) A~E 贪心，构造，线段树，树的子树

Educational Codeforces Round 67 (Rated for Div. 2)

• ### [LightOJ-1356] Prime Independence 二分图+素数分解

数据大，需要用优化的二分图，对每个数求出素因数，不独立的两个数之间就差一个素因数，若 a 去掉这个素因数得到b

• ### C#访问非托管DLL

前段时间写了一个读卡器程序，使用到了一个DLL。在.NET 中调用DLL还是相当简单的。我也是转别人的代码，自己的不写。下面就是代码： C#中要使用动态链接库...

• ### HDU1878 欧拉回路

Problem Description 欧拉回路是指不令笔离开纸面，可画过图中每条边仅一次，且可以回到起点的一条回路。现给定一个图，问是否存在欧拉回路？ Inp...

• ### 20:话题焦点人物

20:话题焦点人物 查看 提交 统计 提问 总时间限制: 1000ms 内存限制: 65536kB描述 微博提供了一种便捷的交流平台。一条微博中，可以提及其它...

• ### DES算法-C/C++实现

闲来无事，车一下轮子，折腾了大半天才搞懂 DESDES 是干毛子的，看了好多博客才弄清楚这个算法的具体原理，真是心累。 只是简单的实现，功能比较简陋，因为参考的...

• ### win2d 画出好看的图形

本文告诉大家，win2d 不需要从零开始做，以前做出来的很多库其实只需要做很小修改就可以做出好看的效果，而且用在 UWP 上。本文修改原先 大神写的 GDI 图...

• ### 洛谷P3763 [TJOI2017]DNA(后缀数组 RMQ)

具体来说，每次比较当前后缀和\(S_0\)的lcp，如果长度\(< N\)的话就从不合法的位置继续匹配