HDU 2389 Rain on your Parade
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
#include<vector>
#include<queue>
#include<cmath>
using namespace std;
#define INF 0x3fffffff
#define maxn 3005
bool vis[maxn];
int Px[maxn], Py[maxn];///Px[i]表示在X集合的i 所匹配Y的值的编号是 P[i]
int dx[maxn], dy[maxn], Head[maxn];///dx记录X集合每个点所在的层, dy同dx
int n, m, depth, k;
struct node
{
int x, y, s;///这个点的坐标集,人移动的速度
}Peo[maxn], Umb[maxn];///分别表示人的坐标点集, 和伞所在位置的点集 struct Edge
{
int v, next;
}e[maxn*maxn]; void Add(int a,int b)
{
e[k].v = b;
e[k].next = Head[a];
Head[a] = k;
k ++;
} void Init()
{
k = ;
memset(e, , sizeof(e));
memset(Px, -, sizeof(Px));
memset(Py, -, sizeof(Py));
memset(Head, -, sizeof(Head));
} bool BFS()
{
queue<int> Q;
depth = INF;///记录深度 memset(dx, -, sizeof(dx));
memset(dy, -, sizeof(dy)); for(int i=; i<n; i++)
{
if(Px[i] == -)
{///以x方的点为源点,进行广搜, 并且是没有加入匹配的点
Q.push(i);
dx[i] = ;
}
} while( Q.size() )
{
int u = Q.front();
Q.pop(); if(dx[u] > depth)///?????
break; for(int i=Head[u]; i != -; i = e[i].next)
{
int v = e[i].v;
if(dy[v] == -)
{
dy[v] = dx[v] + ;
/**说明这个点是没有加入到匹配内的, 也就是说我们找到了一条增广路
(因为我们的源点是没进入匹配内的点,而源点到达的这个点也是没有进入匹配的点,
两个未进入匹配的点,相连了肯定是一对匹配,也就是增广路)*/
if(Py[v] == -)
depth = dy[v];
else///否则源点所连接的就是匹配过的点。
{
dx[Py[v]] = dy[v] + ;
Q.push( Py[v] );
}
}
}
}
return depth != INF;
} bool Find(int u)
{
for(int i=Head[u]; i != -; i = e[i].next)
{
int v = e[i].v;
if(!vis[v] && dx[u] == dy[v] - )///由dx[i]到 dy[v] 可以寻得一条增广路
{
vis[v] = true; /**增广路不在这里,因为我们在dy[v]这个深度的时候已经找到增广路了*/
if(Py[v] != - && dy[v] == depth)
continue; if(Py[v] == - || Find(Py[v]))
{
Py[v] = u;
Px[u] = v;
return true;
}
}
}
return false;
} int solve()
{
int ans = ;
while( BFS() )///确定是否存在增广路
{
memset(vis, false, sizeof(vis));
for(int i=; i<n; i++)
{
if(Px[i] == - && Find(i) )
ans ++;
}
}
return ans;
} int main()
{
int T, time, cas = ;///n 人的个数 m伞的个数
scanf("%d", &T); while(T--)
{
scanf("%d %d", &time, &n);
for(int i=; i<n; i++)
scanf("%d %d %d",&Peo[i].x, &Peo[i].y, &Peo[i].s);
scanf("%d", &m);
for(int i=; i<m; i++)
scanf("%d %d", &Umb[i].x, &Umb[i].y); Init(); for(int i=; i<n; i++)///构图
{
for(int j=; j<m; j++)
{
double len = sqrt( 1.0*(Peo[i].x-Umb[j].x)*(Peo[i].x-Umb[j].x) + (Peo[i].y-Umb[j].y)*(Peo[i].y-Umb[j].y) );
if(len <= time*Peo[i].s)
Add(i,j);
}
}
int ans = solve(); printf("Scenario #%d:\n%d\n\n", cas ++, ans);
}
return ;
} /*
3
1 0 1
1 0 1
0 1 0
*/
HDU 2389 Rain on your Parade的更多相关文章
- 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 最大匹配(模板题)【HK算法】
<题目链接> 题目大意:有m个宾客,n把雨伞,预计时间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 ...
- 【HDOJ】2389 Rain on your Parade
读题显然是二分图匹配,看成guest与umbrella的匹配.匈牙利果断TLE了,其实时间卡的相当紧.HK过的,750ms. /* 2389 */ #include <iostream> ...
- Rain on your Parade
Rain on your Parade Time Limit: 6000/3000 MS (Java/Others) Memory Limit: 655350/165535 K (Java/Ot ...
- Hdu2389 Rain on your Parade (HK二分图最大匹配)
Rain on your Parade Problem Description You’re giving a party in the garden of your villa by the sea ...
随机推荐
- php做站点购物车 你搞懂了吗?
网上购物现已成为时尚,客户选择一个商品将其放入到购物车,然后返回继续购物或者去收银台,这个功能怎样实现呢?今天capucivar就将使用PHP来实现这个购物车的功能. 首先,做一个简单的首页,从数据库 ...
- [Javascript] Proper use of console.assert in JavaScript
Learn about console.assert, which is syntactic sugar for logging an error the console when a given c ...
- 边走边学Nodejs (基础入门篇)
1.什么是Node.js Nodejs ,或者node, 是一个基于ChromeJavaScript执行时建立的平台.用于方便地搭建响应速度快.易于扩展的网络应用.Node.js 使用事件驱动, 非堵 ...
- sublime text3 插件配置
(转) sublme text 全程指引:http://www.cnblogs.com/figure9/p/sublime-text-complete-guide.html 使用Package Con ...
- HDU 5120 Intersection(几何模板题)
题意:给定两个圆环,求两个圆环相交的面积. 思路:由于圆心和半径不一样,分了好多种情况,后来发现只要把两个圆相交的函数写好之后就不需要那么复杂了.两个圆相交的面积的模板如下: double area_ ...
- 鼠标滑动判断与y轴的距离
$(window).scroll(function(){ var y = window.scrollY;//754到达 ...
- linux File Handling commands 'ls'.
ref:Linux / Unix Command: ls NAME ls - list directory contents SYNOPSIS ls [OPTION]... [FILE]... DES ...
- 开始学习<p>标签,添加段落
如果想在网页上显示文章,这时就需要<p>标签了,把文章的段落放到<p>标签中. 语法: <p>段落文本</p> 注意一段文字一个<p>标签, ...
- Direct2D 加载位图
说明: 通过WIC从文件加载位图. 可缩放后加载到内存. 源码: HRESULT LoadImageFormFile( IWICImagingFactory *pWicFactory, ID2D1Re ...
- python challenge 16
前情回顾:上一篇 第16关地址 打开16关,又是一张奇奇怪怪很多点点的图片,应该又是与PIL库有关的. 页面的标题是:let me get this straight.这是英语中的一句俚语,意思是让我 ...