The reachability problem asks to decide if there exists a path from one vertex to another in a digraph. In a grid digraph, the vertices are the points of a two-dimensional square grid, and an edge can occur between a vertex and its immediate horizontal and vertical neighbors only.
Asano and Doerr (CCCG'11) presented the first simultaneous time-space bound for reachability in grid digraphs by solving the problem in polynomial time and $O(n^{1/2 + ε})$ space. In 2018, the space complexity was improved to $\tilde{O}(n^{1/3})$ by Ashida and Nakagawa (SoCG'18).
In this paper, we show that there exists a polynomial-time algorithm that uses $O(n^{1/4 + ε})$ space to solve the reachability problem in a grid digraph containing $n$ vertices. We define and construct a new separator-like device called pseudoseparator to develop a divide-and-conquer algorithm. This algorithm works in a space-efficient manner to solve reachability.
论文ID : 1902.00488标题 : On Solving Reachability in Grid Digraphs using a Pseudoseparator作者 : Rahul Jain, Raghunath Tewari分类 : cs.CC (Computational Complexity), cs.DS (Data Structures and Algorithms)发表时间/会议 : Theory of Computing, Volume 19 (2), 2023 (会议版本发表于FSTTCS 2019)论文链接 : https://arxiv.org/abs/1902.00488 有向图的可达性问题要求判断是否存在从一个顶点到另一个顶点的路径。在网格有向图中,顶点是二维方格网格上的点,边只能存在于顶点与其相邻的水平和垂直邻居之间。Asano和Doerr (CCCG'11)首次给出了网格有向图可达性问题的同时时间-空间界限,在多项式时间和O ( n 1 / 2 + ε ) O(n^{1/2+ε}) O ( n 1/2 + ε ) 空间内解决该问题。2018年,Ashida和Nakagawa (SoCG'18)将空间复杂度改进到O ~ ( n 1 / 3 ) \tilde{O}(n^{1/3}) O ~ ( n 1/3 ) 。本文证明存在一个多项式时间算法,使用O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) 空间解决包含n n n 个顶点的网格有向图的可达性问题。我们定义并构造了一种称为伪分离器的新分离器类设备,开发了一种分治算法,以空间高效的方式解决可达性问题。
理论意义 :有向图可达性问题是NL-complete问题,捕获了非确定性对数空间的复杂性,对计算复杂性理论具有根本重要性应用价值 :许多网络相关问题的算法将其作为子程序使用算法挑战 :标准遍历算法(DFS、BFS)需要线性空间,而Savitch算法虽然只需O ( log 2 n ) O(\log^2 n) O ( log 2 n ) 空间但时间复杂度为n O ( log n ) n^{O(\log n)} n O ( l o g n ) 一般有向图 :Barnes等人的算法达到n / 2 Θ ( log n ) n/2^{\Theta(\sqrt{\log n})} n / 2 Θ ( l o g n ) 空间和多项式时间,但仍不能回答Wigderson提出的是否存在O ( n 1 − ε ) O(n^{1-ε}) O ( n 1 − ε ) 空间多项式时间算法的问题平面有向图 :Imai等人给出O ( n 1 / 2 + ε ) O(n^{1/2+ε}) O ( n 1/2 + ε ) 空间算法,Asano等人改进到O ~ ( n 1 / 2 ) \tilde{O}(n^{1/2}) O ~ ( n 1/2 ) 空间网格有向图 :作为平面有向图的子类,Asano-Doerr算法达到O ( n 1 / 2 + ε ) O(n^{1/2+ε}) O ( n 1/2 + ε ) 空间,Ashida-Nakagawa改进到O ( n 1 / 3 ) O(n^{1/3}) O ( n 1/3 ) 空间本文旨在进一步改进网格有向图可达性问题的空间复杂度,通过引入新的伪分离器概念,实现O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) 空间复杂度的突破。
主要理论结果 :证明了存在多项式时间算法使用O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) 空间解决n n n 顶点网格有向图的可达性问题新概念引入 :定义并构造了伪分离器(pseudoseparator)概念,这是一种新的分离器类设备创新算法设计 :开发了基于伪分离器的分治算法,有效利用辅助图的交叉性质技术突破 :显著改进了现有最好结果,从O ( n 1 / 3 ) O(n^{1/3}) O ( n 1/3 ) 提升到O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) 输入 :m × m m×m m × m 网格有向图G G G ,其中顶点集为[ m ] × [ m ] = { 0 , 1 , . . . , m } × { 0 , 1 , . . . , m } [m]×[m] = \{0,1,...,m\}×\{0,1,...,m\} [ m ] × [ m ] = { 0 , 1 , ... , m } × { 0 , 1 , ... , m } ,边( u , v ) (u,v) ( u , v ) 存在当且仅当∣ u 1 − v 1 ∣ + ∣ u 2 − v 2 ∣ = 1 |u_1-v_1|+|u_2-v_2|=1 ∣ u 1 − v 1 ∣ + ∣ u 2 − v 2 ∣ = 1 ,以及两个查询顶点s , t s,t s , t
输出 :判断是否存在从s s s 到t t t 的有向路径
约束 :算法必须在多项式时间内完成,空间使用量为O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) ,其中n = ( m + 1 ) 2 n=(m+1)^2 n = ( m + 1 ) 2
将m × m m×m m × m 网格有向图G G G 分割成m 2 α m^{2α} m 2 α 个子网格,每个子网格大小为m 1 − α × m 1 − α m^{1-α}×m^{1-α} m 1 − α × m 1 − α :
对于子网格G [ i , j ] G[i,j] G [ i , j ] ,构造辅助图A u x α ( G ) [ i , j ] Aux_α(G)[i,j] A u x α ( G ) [ i , j ] ,顶点集为边界顶点 如果子网格内存在从u u u 到v v v 的路径,则在辅助图中添加边( u , v ) (u,v) ( u , v ) 最终辅助图A u x α ( G ) Aux_α(G) A u x α ( G ) 包含所有子网格的辅助图 定义4.1 :对于辅助图A u x α ( G ) Aux_α(G) A u x α ( G ) 的导出子图H H H ,子图Q Q Q 是f ( h ) f(h) f ( h ) -伪分离器,当且仅当H ⋄ Q H⋄Q H ⋄ Q 中每个连通分量的大小至多为f ( h ) f(h) f ( h ) ,其中H ⋄ Q H⋄Q H ⋄ Q 表示从H H H 中移除Q Q Q 的顶点和与Q Q Q 中边相交的边后得到的图。
构造过程 :
构造p l a n a r ( H ) planar(H) pl ana r ( H ) :选择H H H 中不相交的边的最大子集 使用算法1完成边界,然后三角化得到H ~ \tilde{H} H ~ 应用Imai等人的平面分离器算法得到顶点集S S S 构造伪分离器p s e p ( H ) psep(H) p se p ( H ) ,包含S S S 中的顶点和相关的屏蔽边 引理3.5 :如果辅助图中两条边e 1 = ( u 1 , v 1 ) e_1=(u_1,v_1) e 1 = ( u 1 , v 1 ) 和e 2 = ( u 2 , v 2 ) e_2=(u_2,v_2) e 2 = ( u 2 , v 2 ) 相交,则辅助图也包含边( u 1 , v 2 ) (u_1,v_2) ( u 1 , v 2 ) 和( u 2 , v 1 ) (u_2,v_1) ( u 2 , v 1 ) 。
引理3.6 :如果e 1 e_1 e 1 和e 2 e_2 e 2 都与边f = ( x , y ) f=(x,y) f = ( x , y ) 相交,且e 1 e_1 e 1 比e 2 e_2 e 2 更接近x x x ,则边( u 1 , v 2 ) (u_1,v_2) ( u 1 , v 2 ) 也在辅助图中。
Input: 辅助图的导出子图H,查询顶点x,y
1. 如果|H| ≤ m^{1/8},使用DFS直接解决
2. 否则:
a. 构造h^{1-β}-伪分离器Q
b. 维护visited数组标记Q中的顶点和边
c. 初始化visited[x] := 1
d. 进行h次迭代,每次更新visited数组
e. 返回visited[y]是否为1
Input: 网格有向图G,查询顶点s,t
1. 如果G小于m^{1/8}×m^{1/8},使用DFS
2. 否则调用ImplicitAuxReach(Aux_α(G),s,t)
3. 当查询辅助图中的边时,递归调用GridReach
伪分离器概念 :扩展了传统分离器,允许使用边来"分离"图,利用辅助图的交叉性质交叉性质利用 :巧妙利用引理3.5和3.6,使得路径跨越不同组件时只需存储少量顶点递归结构设计 :通过适当选择参数α α α 和β β β ,实现空间复杂度的优化隐式图表示 :不显式存储辅助图,而是按需递归计算边的存在性本文主要进行理论分析,通过数学证明建立算法的正确性和复杂度界限。
空间复杂度 :S ( m ) = S ( m 1 − α ) + O ( ( m 1 + α ) 1 / 2 + β / 2 ) S(m) = S(m^{1-α}) + O((m^{1+α})^{1/2+β/2}) S ( m ) = S ( m 1 − α ) + O (( m 1 + α ) 1/2 + β /2 ) 时间复杂度 :T ( m ) = q ( m ) ⋅ T ( m 1 − α ) + p ( m ) T(m) = q(m)·T(m^{1-α}) + p(m) T ( m ) = q ( m ) ⋅ T ( m 1 − α ) + p ( m ) ,其中p ( m ) p(m) p ( m ) 和q ( m ) q(m) q ( m ) 是多项式参数选择 :对于任意常数ε > 0 ε > 0 ε > 0 ,选择适当的α α α 和β β β 使得空间复杂度为O ( m 1 / 2 + ε ) O(m^{1/2+ε}) O ( m 1/2 + ε ) 通过归纳法证明AuxReach算法的正确性,关键在于证明路径从一个组件到另一个组件时,算法能正确标记相应的顶点或边。
定理1.1 :对于每个ε > 0 ε > 0 ε > 0 ,存在多项式时间算法使用O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) 空间解决n n n 顶点网格有向图的可达性问题。
空间复杂度 :从之前的O ( n 1 / 3 ) O(n^{1/3}) O ( n 1/3 ) 改进到O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) 时间复杂度 :保持多项式时间递归深度 :常数深度,确保空间的有效重用引理4.8 :证明了构造的p s e p ( H ) psep(H) p se p ( H ) 确实是h 1 − β h^{1-β} h 1 − β -伪分离器引理5.1 :证明了AuxReach算法的正确性引理5.2 :建立了算法的空间和时间复杂度界限Savitch算法:O ( log 2 n ) O(\log^2 n) O ( log 2 n ) 空间,n O ( log n ) n^{O(\log n)} n O ( l o g n ) 时间 Barnes等人:n / 2 Θ ( log n ) n/2^{\Theta(\sqrt{\log n})} n / 2 Θ ( l o g n ) 空间,多项式时间 平面有向图 :Imai等人:O ( n 1 / 2 + ε ) O(n^{1/2+ε}) O ( n 1/2 + ε ) 空间 Asano等人:O ~ ( n 1 / 2 ) \tilde{O}(n^{1/2}) O ~ ( n 1/2 ) 空间 其他图类 :高亏格图:O ~ ( n 2 / 3 g 1 / 3 ) \tilde{O}(n^{2/3}g^{1/3}) O ~ ( n 2/3 g 1/3 ) 空间 H H H -minor-free图:O ~ ( n 2 / 3 ) \tilde{O}(n^{2/3}) O ~ ( n 2/3 ) 空间层次平面图:O ( n ε ) O(n^ε) O ( n ε ) 空间 Asano-Doerr (2011):O ( n 1 / 2 + ε ) O(n^{1/2+ε}) O ( n 1/2 + ε ) 空间 Ashida-Nakagawa (2018):O ( n 1 / 3 ) O(n^{1/3}) O ( n 1/3 ) 空间 本文 (2023):O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) 空间 本文成功将网格有向图可达性问题的空间复杂度从O ( n 1 / 3 ) O(n^{1/3}) O ( n 1/3 ) 改进到O ( n 1 / 4 + ε ) O(n^{1/4+ε}) O ( n 1/4 + ε ) ,这是该问题空间复杂度的显著改进。
伪分离器 :提供了一种新的图分解工具,适用于非平面图交叉性质利用 :巧妙利用辅助图的结构性质算法设计 :展示了如何在保持多项式时间的同时大幅减少空间使用常数因子 :算法中涉及多个常数,实际常数可能较大适用范围 :仅适用于网格有向图,不能直接推广到一般平面图下界缺失 :未提供相应的下界结果平面图推广 :虽然网格图可达性能归约到平面图,但由于二次膨胀,直接改进平面图算法可能更有效下界研究 :建立网格图可达性问题的空间下界实用算法 :开发具有更好常数因子的实用算法理论突破 :在重要问题上取得显著的复杂度改进技术创新 :伪分离器概念具有独创性,为相关问题提供新思路严谨证明 :数学证明完整严谨,技术细节处理得当写作清晰 :论文结构清晰,概念定义准确实用性限制 :算法较为复杂,常数因子可能很大,实用价值有限推广困难 :方法难以直接推广到更一般的图类缺乏实验 :纯理论工作,缺乏实际性能评估学术价值 :为计算复杂性理论做出重要贡献技术影响 :伪分离器概念可能启发相关研究后续工作 :为进一步改进奠定基础主要适用于理论计算机科学研究,特别是:
空间复杂度理论研究 图算法设计 几何算法分析 论文引用了22篇重要文献,涵盖了可达性问题、平面图算法、分离器理论等相关领域的关键工作,为研究提供了坚实的理论基础。
这篇论文在网格有向图可达性问题上取得了重要的理论突破,虽然实用价值有限,但其技术创新和理论贡献对计算复杂性理论具有重要意义。伪分离器的概念和相关技术可能会启发更多相关研究。