HDU 2389 Rain on your Parade 最大匹配(模板题)【HK算法】
<题目链接>
题目大意:
有m个宾客,n把雨伞,预计时间t后将会下大雨,告诉你每个宾客的位置和速度,每把雨伞的位置,问你最多几个宾客能够拿到伞。
解题分析:
本题就是要我们求人与伞之间的最大匹配,但是数据量较大,匈牙利算法复杂度为$O(n \times m)$,会超时,所以这里用的是复杂度为$O(\sqrt{n} \times m)$的HK算法。
当然,本题也可以用Dinic最大流来求解最大匹配,并且这种方法求得到最大匹配的时间复杂度也为$O(sqrt(n)*m)$
#include <bits/stdc++.h>
using namespace std;
#define INF 0x3f3f3f3f
const int M = +;
struct point{
double x,y,s;
}a[M],b[M];
bool map[M][M];
bool vis[M];
int dx[M],dy[M];/*dx[i]表示左集合i顶点的距离编号,dy[i]表示右集合i顶点的距离编号,还有就是这里就是bfs找增广路用的数组 对于u-->v可达就有dy[v] = dx[u] + 1 */
int cx[M],cy[M];
int dis;
int n,T,m;
bool searchpath()
{
queue<int>q;
dis=INF;
memset(dx,-,sizeof(dx));
memset(dy,-,sizeof(dy));
for(int i=;i<=n;i++)
{
//cx[i]表示左集合i顶点所匹配的右集合的顶点序号
if(cx[i]==-)
{
//将未遍历的节点 入队 并初始化次节点距离为0
q.push(i);
dx[i]=;
}
}
//广度搜索增广路径
while(!q.empty())
{
int u=q.front();
q.pop();
if(dx[u]>dis) break;
//取右侧节点,所以用m。
for(int v=;v<=m;v++)
{
//右侧节点的增广路径的距离
if(map[u][v]&&dy[v]==-)
{
dy[v]=dx[u]+;//v对应的距离 为u对应距离加1
if(cy[v]==-) dis=dy[v];
else
{
dx[cy[v]]=dy[v]+;
q.push(cy[v]);
}
}
}
}
return dis!=INF;
}
//寻找路径 深度搜索
int dfs(int u)
{
for(int v=;v<=m;v++)
{
//如果该点没有被遍历过 并且距离为上一节点+1
if(map[u][v]&&!vis[v]&&dy[v]==dx[u]+)
{
//对该点染色,即是进行标记
vis[v]=true;
if(cy[v]!=-&&dy[v]==dis)
continue;
if(cy[v]==-||dfs(cy[v]))
{
cx[u]=v;
cy[v]=u;
return ;
}
}
}
return ;
}
//得到最大匹配的数目
int MaxMatch()
{
memset(cx,-,sizeof(cx));
memset(cy,-,sizeof(cy));
int res=;
while(searchpath())
{
memset(vis,false,sizeof(vis));
for(int i=;i<=n;i++)
{
if(cx[i]==-)
res+=dfs(i);
}
}
return res;
}
int main()
{
int t,k=;
scanf("%d",&t);
while(t--)
{
memset(map,false,sizeof(map));
scanf("%d",&T);
scanf("%d",&n);
for(int i=;i<=n;i++)
scanf("%lf%lf%lf",&a[i].x,&a[i].y,&a[i].s);
scanf("%d",&m);
for(int i=;i<=m;i++)
scanf("%lf%lf",&b[i].x,&b[i].y);
for(int i=;i<=n;i++)
for(int j=;j<=m;j++)
{
if(a[i].s*T>=sqrt((a[i].x-b[j].x)*(a[i].x-b[j].x)+(a[i].y-b[j].y)*(a[i].y-b[j].y)))
map[i][j]=true;
}
int ans=MaxMatch();
printf("Scenario #%d:\n%d\n\n",k++,ans);
}
}
HK
HDU 2389 Rain on your Parade 最大匹配(模板题)【HK算法】的更多相关文章
- HDU 2389 Rain on your Parade / HUST 1164 4 Rain on your Parade(二分图的最大匹配)
HDU 2389 Rain on your Parade / HUST 1164 4 Rain on your Parade(二分图的最大匹配) Description You're giving a ...
- HDU 2389 ——Rain on your Parade——————【Hopcroft-Karp求最大匹配、sqrt(n)*e复杂度】
Rain on your Parade Time Limit:3000MS Memory Limit:165535KB 64bit IO Format:%I64d & %I64 ...
- HDU 2389 Rain on your Parade(二分匹配,Hopcroft-Carp算法)
Rain on your Parade Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 655350/165535 K (Java/Ot ...
- HDU 2389 Rain on your Parade
大意:在一个二维坐标系上有nx个人和ny把伞,每个人都有自己的移动速度,问有多少人可以再 time 时间内移动到不同的雨伞处(不允许两个人共用一把伞). 输入数据: 第一行是一个T代表T组测试数据 ...
- Hdu 3289 Rain on your Parade (二分图匹配 Hopcroft-Karp)
题目链接: Hdu 3289 Rain on your Parade 题目描述: 有n个客人,m把雨伞,在t秒之后将会下雨,给出每个客人的坐标和每秒行走的距离,以及雨伞的位置,问t秒后最多有几个客人可 ...
- HDOJ 2389 Rain on your Parade
HK.... Rain on your Parade Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 655350/165535 K ...
- HDU 4280:Island Transport(ISAP模板题)
http://acm.hdu.edu.cn/showproblem.php?pid=4280 题意:在最西边的点走到最东边的点最大容量. 思路:ISAP模板题,Dinic过不了. #include & ...
- HDU 4347 - The Closest M Points - [KDTree模板题]
本文参考: https://www.cnblogs.com/GerynOhenz/p/8727415.html kuangbin的ACM模板(新) 题目链接:http://acm.hdu.edu.cn ...
- HDU 2089 不要62(数位dp模板题)
http://acm.hdu.edu.cn/showproblem.php?pid=2089 题意:求区间内不包含4和连续62的数的个数. 思路: 简单的数位dp模板题.给大家推荐一个好的讲解博客.h ...
随机推荐
- 《Oracle DBA工作笔记:运维、数据迁移与性能调优》 PDF 下载
一:下载途径 二:本书图样 三:本书目录 第1篇 数据库运维篇第1章 数据库安装配置1.1 安装前的准备 11.2 安装数据库软件 51.2.1 方法1:OUI安装 61.2.2 方法2:静默安装 8 ...
- Confluence 6 使用 JMX 界面实时监控
使用 JMX 界面(Java Management Extensions API ),你可以实时的查看你 Confluence 运行实例的状态. JMX 使用的对象被称 MBeans (Managed ...
- jquery 中多选和全选
- SpringCloud注册中心(Eureka)
注册中心的作用就是用来方便接口暴露出来的一个管理工具,如果所有的接口只是暴露出来,没有一个同一的管理工具,又有谁知道你暴露了哪些接口了,Eureka是一个springcloud的服务注册和发现模块 搭 ...
- BeautifulSoup的基本操作
>>> from bs4 import BeautifulSoup #导入 >>> soup = BeautifulSoup(url.content," ...
- Django框架之Form组件
一.初探Form组件 在介绍Form组件之前,让大家先看看它强大的功能吧!Go... 下面我们来看看代码吧! 1.创建Form类 from django.forms import Form from ...
- PDF怎么去除页眉页脚,PDF页眉页脚编辑方法
我们在使用文件的时候需要编辑页眉页脚的时候,这个时候我们应该怎么做呢,相信别的文件大家都知道怎么编辑了,PDF文件大家都知道吗,最开始接触这个文件的时候小编觉得很难,之后找到技巧之后也并没有很难,今天 ...
- less 写关键帧动画
@keyframes 关键帧动画写在less里的时候,务必要写在所有的{}之外,不能被{}包裹 甚至务必写在代码的最后 不然报错 Compilation resulted in incorrect C ...
- 论文阅读笔记四:CTPN: Detecting Text in Natural Image with Connectionist Text Proposal Network(ECCV2016)
前面曾提到过CTPN,这里就学习一下,首先还是老套路,从论文学起吧.这里给出英文原文论文网址供大家阅读:https://arxiv.org/abs/1609.03605. CTPN,以前一直认为缩写一 ...
- ***小程序wx.getUserInfo不能弹出授权窗口后的解决方案
微信更新api后,wx.getUserInfo在开发和体验版本都不能弹出授权窗口.微信文档说明: 注意:此接口有调整,使用该接口将不再出现授权弹窗,请使用 <button open-type=& ...