【HDU】2389 Rain on your Parade 二分匹配 Hopcroft-Krap算法

2024-09-05 15:32

本文主要是介绍【HDU】2389 Rain on your Parade 二分匹配 Hopcroft-Krap算法,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

传送门:【HDU】2389 Rain on your Parade


题目分析:

这题目非要我学Hopcroft-Krap= =||。。普通的DFS版的二分匹配不行,最大流又爆内存。。不得不学更好的算法了。

二分匹配的其他性质我也不多说了,不会的自行搜索,网上很多的。

现在我主要对该算法的实现发表一下自己的见解。(算法复杂度的证明不会,论文没看太懂)

该算法的核心思想是通过bfs寻找多条相同长度的最短增广路实现多路增广,那么怎么进行bfs?

首先,将还未匹配的X集合的顶点加入队列。然后用这个队列中的元素挨个查找对应的Y集合的元素。

如果查找到的Y集合的元素是没被查找过的,如果其已经被覆盖,那么将与他匹配的X集合的元素入队;如果是没有被覆盖的,那么就确定了这次找到的最短增广路的长度了,之后所有大于等于该长度的X集合的节点都不用再用来继续查找了,因为不可能找到长度等于最短增广路长度的增广路了。当然找增广路的时候需要给所有的顶点距离标号,表示到起点的距离。

通过bfs与处理以后,在dfs中只走标号差等于1的路,可以节省很多时间。算法证明的复杂度是V^0.5*E。


代码如下:


#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;#define REP( i , n ) for ( int i = 0 ; i < n ; ++ i )
#define REPF( i , a , b ) for ( int i = a ; i <= b ; ++ i )
#define REPV( i , a , b ) for ( int i = a ; i >= b ; -- i )
#define clear( a , x ) memset ( a , x , sizeof a )const int MAXN = 3005 ;
const int MAXE = 10000000 ;
const int INF = 0x3f3f3f3f ;struct Edge {int v , n ;Edge () {}Edge ( int var , int next ) :v ( var ) , n ( next ) {}
} ;struct Node {int x , y , v ;void input () {scanf ( "%d%d%d" , &x , &y , &v ) ;}
} ;Edge E[MAXE] ;
Node A[MAXN] ;
int H[MAXN] , cntE ;
int Lx[MAXN] , Ly[MAXN] ;
int dx[MAXN] , dy[MAXN] ;
bool vis[MAXN] ;
int x[MAXN] , y[MAXN] ;
int Q[MAXN << 1] , head , tail ;
int t , n , m ;
int dis ;void addedge ( int u , int v ) {E[cntE] = Edge ( v , H[u] ) ;H[u] = cntE ++ ;
}int Hopcroft_Krap () {dis = INF ;head = tail = 0 ;clear ( dx , -1 ) ;clear ( dy , -1 ) ;REP ( i , n )if ( Lx[i] == -1 ) {dx[i] = 0 ;Q[tail ++] = i ;}while ( head != tail ) {int u = Q[head ++] ;if ( dx[u] >= dis )continue ;for ( int i = H[u] ; ~i ; i = E[i].n ) {int v = E[i].v ;if ( dy[v] == -1 ) {dy[v] = dx[u] + 1 ;if ( Ly[v] == -1 )dis = dy[v] ;else {dx[Ly[v]] = dy[v] + 1 ;Q[tail ++] = Ly[v] ;}}}}return dis != INF ;
}int find ( int u ) {for ( int i = H[u] ; ~i ; i = E[i].n ) {int v = E[i].v ;if ( !vis[v] && dy[v] == dx[u] + 1 ) {vis[v] = 1 ;if ( ~Ly[v] && dy[v] == dis )continue ;else if ( Ly[v] == -1 || find ( Ly[v] ) ) {Lx[u] = v ;Ly[v] = u ;return 1 ;}}}return 0 ;
}int match () {int ans = 0 ;clear ( Lx , -1 ) ;clear ( Ly , -1 ) ;while ( Hopcroft_Krap () ) {clear ( vis , 0 ) ;REP ( i , n )if ( Lx[i] == -1 )ans += find ( i ) ;}return ans ;
}int dist ( int i , int j ) {int X = A[i].x - x[j] ;int Y = A[i].y - y[j] ;return X * X + Y * Y ;
}void solve () {cntE = 0 ;clear ( H , -1 ) ;scanf ( "%d%d" , &t , &n ) ;REP ( i , n )A[i].input () ;scanf ( "%d" , &m ) ;REP ( i , m )scanf ( "%d%d" , &x[i] , &y[i] ) ;REP ( i , n ) {int tmp = t * t * A[i].v * A[i].v ;REP ( j , m )if ( dist ( i , j ) <= tmp )addedge ( i , j ) ;}printf ( "%d\n\n" , match () ) ;
}int main () {int T , cas ;for ( scanf ( "%d" , &T ) , cas = 1 ; cas <= T ; ++ cas ) {printf ( "Scenario #%d:\n" , cas ) ;solve () ;}return 0 ;
}


这篇关于【HDU】2389 Rain on your Parade 二分匹配 Hopcroft-Krap算法的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1139352

相关文章

Nginx location匹配模式与规则详解

《Nginxlocation匹配模式与规则详解》:本文主要介绍Nginxlocation匹配模式与规则,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、环境二、匹配模式1. 精准模式2. 前缀模式(不继续匹配正则)3. 前缀模式(继续匹配正则)4. 正则模式(大

Java 正则表达式URL 匹配与源码全解析

《Java正则表达式URL匹配与源码全解析》在Web应用开发中,我们经常需要对URL进行格式验证,今天我们结合Java的Pattern和Matcher类,深入理解正则表达式在实际应用中... 目录1.正则表达式分解:2. 添加域名匹配 (2)3. 添加路径和查询参数匹配 (3) 4. 最终优化版本5.设计思

openCV中KNN算法的实现

《openCV中KNN算法的实现》KNN算法是一种简单且常用的分类算法,本文主要介绍了openCV中KNN算法的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的... 目录KNN算法流程使用OpenCV实现KNNOpenCV 是一个开源的跨平台计算机视觉库,它提供了各

Python中使用正则表达式精准匹配IP地址的案例

《Python中使用正则表达式精准匹配IP地址的案例》Python的正则表达式(re模块)是完成这个任务的利器,但你知道怎么写才能准确匹配各种合法的IP地址吗,今天我们就来详细探讨这个问题,感兴趣的朋... 目录为什么需要IP正则表达式?IP地址的基本结构基础正则表达式写法精确匹配0-255的数字验证IP地

浅谈配置MMCV环境,解决报错,版本不匹配问题

《浅谈配置MMCV环境,解决报错,版本不匹配问题》:本文主要介绍浅谈配置MMCV环境,解决报错,版本不匹配问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录配置MMCV环境,解决报错,版本不匹配错误示例正确示例总结配置MMCV环境,解决报错,版本不匹配在col

springboot+dubbo实现时间轮算法

《springboot+dubbo实现时间轮算法》时间轮是一种高效利用线程资源进行批量化调度的算法,本文主要介绍了springboot+dubbo实现时间轮算法,文中通过示例代码介绍的非常详细,对大家... 目录前言一、参数说明二、具体实现1、HashedwheelTimer2、createWheel3、n

详解nginx 中location和 proxy_pass的匹配规则

《详解nginx中location和proxy_pass的匹配规则》location是Nginx中用来匹配客户端请求URI的指令,决定如何处理特定路径的请求,它定义了请求的路由规则,后续的配置(如... 目录location 的作用语法示例:location /www.chinasem.cntestproxy

SpringBoot实现MD5加盐算法的示例代码

《SpringBoot实现MD5加盐算法的示例代码》加盐算法是一种用于增强密码安全性的技术,本文主要介绍了SpringBoot实现MD5加盐算法的示例代码,文中通过示例代码介绍的非常详细,对大家的学习... 目录一、什么是加盐算法二、如何实现加盐算法2.1 加盐算法代码实现2.2 注册页面中进行密码加盐2.

Java时间轮调度算法的代码实现

《Java时间轮调度算法的代码实现》时间轮是一种高效的定时调度算法,主要用于管理延时任务或周期性任务,它通过一个环形数组(时间轮)和指针来实现,将大量定时任务分摊到固定的时间槽中,极大地降低了时间复杂... 目录1、简述2、时间轮的原理3. 时间轮的实现步骤3.1 定义时间槽3.2 定义时间轮3.3 使用时

如何通过Golang的container/list实现LRU缓存算法

《如何通过Golang的container/list实现LRU缓存算法》文章介绍了Go语言中container/list包实现的双向链表,并探讨了如何使用链表实现LRU缓存,LRU缓存通过维护一个双向... 目录力扣:146. LRU 缓存主要结构 List 和 Element常用方法1. 初始化链表2.