hdu2389 Rain on your Parade--HK算法 & 最大匹配数
来源:互联网 发布:台湾知见旅行社 编辑:程序博客网 时间:2024/06/04 19:51
原题链接:http://acm.hdu.edu.cn/showproblem.php?pid=2389
题意:pNum个人,uNum把伞,给定这些人和伞的二维坐标,和每个人的行进速度。问在t时间内,最多有多少人能拿到伞?
建二分图,左边是人,与右边是伞,如果这个人到伞的时间不大于t就连线,最后求最大匹配就行,不过这题用匈牙利算法会超时,所以用HK算法。
#define _CRT_SECURE_NO_DEPRECATE #include<cstring>#include<iostream>#include<algorithm>#include<cmath>#include<queue>#include<vector>#include<string>#define INF 99999999;using namespace std;struct People{double x;double y;double t;};struct Umbrella{double x;double y;};int T;int time;//约束时间int pNum;//人的数目int uNum;//伞的数目People people[3005];Umbrella umbrella[3005];int dis;int cx[3005];int cy[3005];int dx[3005];int dy[3005];bool g[3005][3005];bool vis[3005];double distance(People p, Umbrella u){return sqrt((p.x - u.x)*(p.x - u.x) + (p.y - u.y)*(p.y - u.y));}bool searchPath(){dis = INF;queue<int> Q;memset(dx, -1, sizeof(dx));memset(dy, -1, sizeof(dy));for (int i = 0; i < pNum; i++){if (cx[i] == -1){Q.push(i);dx[i] = 0;}}while (!Q.empty()){int u = Q.front();Q.pop();if (dx[u] > dis)break;for (int v = 0; v < uNum; v++){if (g[u][v] && dy[v] == -1){dy[v] = dx[u] + 1;if (cy[v] == -1)dis = dy[v];else{dx[cy[v]] = dy[v] + 1;Q.push(cy[v]);}}}}return dis != INF;}int dfs(int u){for (int v = 0; v < uNum; v++){if (g[u][v] && dy[v] == dx[u] + 1 && !vis[v]){vis[v] = 1;if (cy[v] != -1 && dy[v] == dis)continue;if (cy[v] == -1 || dfs(cy[v])){cx[u] = v;cy[v] = u;return 1;}}}return 0;}int HK(){int ans = 0;memset(cx, -1, sizeof(cx));memset(cy, -1, sizeof(cy));while (searchPath()){memset(vis, 0, sizeof(vis));for (int u = 0; u < pNum; u++){if (cx[u] == -1)ans += dfs(u);}}return ans;}int main(){scanf("%d", &T);int cas = 1;while (T--){scanf("%d%d", &time, &pNum);for (int i = 0; i < pNum; i++)scanf("%lf%lf%lf", &people[i].x, &people[i].y, &people[i].t);scanf("%d", &uNum);for (int i = 0; i < uNum; i++)scanf("%lf%lf", &umbrella[i].x, &umbrella[i].y);memset(g, 0, sizeof(g));for (int i = 0; i < pNum; i++)for (int j = 0; j < uNum; j++)if (distance(people[i], umbrella[j]) / people[i].t <= time)g[i][j] = 1;printf("Scenario #%d:\n%d\n\n", cas++, HK());}return 0;}
1 0
- hdu2389 Rain on your Parade--HK算法 & 最大匹配数
- hdu2389 Rain on your Parade【二分图最大匹配-HK算法】
- HDU2389 Rain on your Parade 【二分图最大匹配+HK算法】
- Hdu2389 Rain on your Parade (HK二分图最大匹配)
- HDU2389 Rain on your Parade(HK算法)
- HDU2389 Rain on your Parade(二分图匹配 ,HK算法 )
- HDU2389---Rain on your Parade (二分图匹配-HK算法)
- hdu2389 Rain on your Parade (最大匹配Hopcroft--Karp)
- hdu2389 Rain on your Parade(二分匹配hopcroft-carp算法)
- hdu 2389 Rain on your Parade 二分匹配 HK算法
- 【二分匹配】 HDU 2389 Rain on your Parade HK算法
- hdu 2389 Rain on your Parade(二分最大匹配HK算法)
- HDU2389 Rain on your Parade 二分图匹配——Hopcroft_Karp算法
- HDU2389 Rain on your Parade 二分匹配 Hopcroft-Carp的算法+模版
- HDU2389-Rain on your Parade(二分图匹配Hopcroft-Karp算法)
- hdu 2389 Rain on your Parade 二分图匹配+HK算法
- HDU 2389 Rain on your Parade(HK算法)
- hdoj 2389 Rain on your Parade 【HK】
- 事件冒泡
- shape的使用
- Python——我的第一个Python程序cmd直接输入
- Myeclipse-Java head space
- POJ_3627_贪心
- hdu2389 Rain on your Parade--HK算法 & 最大匹配数
- mysqldump: dumper has insufficent privileges to SHOW CREATE FUNCTION
- 干货:机器学习领域的几种主要学习方式,如深度学习
- 如何打开Mac OSX原生的读写NTFS功能(支持OS X EI,秒杀NTFS+)
- uva 11270 插头dp入门
- eclipse如何搭载Tomcat
- Json解析
- unity JavaScript 与 c# 脚本之间的通信
- Android Studio (AS)快捷键