算法合集之动态树问题及其应用.pptx

Dynamic Trees Problem,and its applications湖南省长郡中学 袁昕颢xinhaoyuanatgmaildotcomOverview1.动态树问题n给出动态树问题的基本形式.2.解决动态树问题n提出新的Rake&Compress方法.3.动态树问题的应用n用最大流算法来说明动态树问题的应用.Part I.Dynamic Trees ProblemDynamic Trees Problemn动态树问题(Dynamic Trees Problem)是图论中一类非常重要的经典问题.许多图论算法,尤其是在线动态算法都将其作为瓶颈问题.n研究和解决该问题具有很高的理论价值和实际价值.n什么是动态树问题呢?Dynamic Trees Problemn维护一个包含N个点的森林,并且支持形态和权值信息的操作.n形态信息Dynamic Trees Problemn维护一个包含N个点的森林,并且支持形态和权值信息的操作.n形态信息nLink(u,v)添加边(u,v)Dynamic Trees Problemn维护一个包含N个点的森林,并且支持形态和权值信息的操作.n形态信息nLink(u,v)添加边(u,v)nCut(u,v)删除边(u,v)Dynamic Trees Problemn维护一个包含N个点的森林,并且支持形态和权值信息的操作.n形态信息nLink(u,v)添加边(u,v)nCut(u,v)删除边(u,v)nFind(u)找到u所在的树.Dynamic Trees Problemn维护一个包含N个点的森林,并且支持形态和权值信息的操作.n权值信息Dynamic Trees Problemn维护一个包含N个点的森林,并且支持形态和权值信息的操作.n权值信息n路径操作:对一条简单路径上的所有对象进行操作Dynamic Trees Problemn维护一个包含N个点的森林,并且支持形态和权值信息的操作.n权值信息n路径操作:对一条简单路径上的所有对象进行操作n树操作:对一棵树内的所有对象进行操作现有结果基本原理Euler TourPath DecomposingDivide and Conquer相关实现Euler Tour TreesST-Trees1,2(Self Adjust)Top-Trees3,5局限性不支持路径操作不支持树权操作常数过大理论补充n对于一个完整的动态树问题,目前公认的下界是O(log2N)per operation,并已经被上述方法达到.n但是由于巨大的常数因子,动态树在实践中并没有发挥应有的作用.n动态树问题仍然没有完美解决,并且仍然处在热烈讨论中.Part II.Solving Dynamic Trees ProblemNew Idean在这里,我向大家介绍一种新的解决动态树问题的思路.这种思路简单,而且,可以得到效率非常高的具体实现.I.树,与其平面刻画.n一棵树的平面刻画,直观地说就是将一棵树的点和边画在平面上.边与边仅在顶点处相交.I.树,与其平面刻画.n一棵树的平面刻画,直观地说就是将一棵树的点和边画在平面上.边与边仅在顶点处相交.n确定一棵树的平面刻画,等价于确定这棵树的Euler Tour.II.等价映射n事实上,所有解决动态树问题的方法,归根结底都使用等价映射的基本思想.n即,将任意形态的树(原树)映射到度限制,深度平均的新树(像树).III.Rake&Compressn这里介绍一种Rake&Compress5,6方法.n即将原树映射到一棵Rake&Compress Trees(Abbr.R&C Trees).nR&C Trees由Rake节点和Compress节点组成.1.Rake NodesnRake节点i是原树中以某节点为根的有根子树的映射.Root1.Rake NodesnRake节点i是原树中以某节点为根的有根子树的映射.n特别地,如果该节点仅包含根本身,那么该Rake节点没有后继(叶子节点).否则令Next(i)表示i所代表的除了根以外的其它点组成的集合.Root2.Compress NodesnCompress节点j,是原树中以某条路径为根的有根子树的映射.se2.Compress NodesnCompress节点j,是原树中以某条路径为根的有根子树的映射n特别地,如果路径长度为1.那么该Compress节点没有后继.否则令Next(j)表示j代表的路径上的非端点集合.se3.R&C Treesn对于一个非叶子Rake/Compress节点i,Next(i)非空.n对于每个Next(i)中的元素j.我们采用如下方法划分节点i:1 Rake节点的划分n令r表示i的根.rj1 Rake节点的划分n令r表示i的根.n将路径jr作为新的Compress节点.rj1 Rake节点的划分n令r表示i的根.n将路径jr作为新的Compress节点.n将j和j的子孙作为新的Rake节点.rj1 Rake节点的划分n令r表示i的根.n将路径jr作为新的Compress节点.n将j和j的子孙作为新的Rake节点.ni中路径jr的左手方向和右手方向各为一个新的Rake节点.rj2 Compress节点的划分n令s,e分别表示i中路径的头和尾.sej2 Compress节点的划分n令s,e分别表示i中路径的头和尾.nsj,je分别成为新的Compress节点sej2 Compress节点的划分n令s,e分别表示i中路径的头和尾.nsj,je分别成为新的Compress节点nj的其他子孙被划分成两部分,分别作为新的Rake节点.sejR&C Treesn约定,一个有根树对应R&C Trees就是整个树的Rake Node.n这样的分解方式本身就保证度的限制一个节点最多被剖分出4个子节点)n我们以右图为例，展示一棵原树如何被映射到像树。Level 1选取I作为第1层剖分点，原Rake节点被剖分成4个部分，4个部分分别剖分。Level 2选取B,C,D,L作为第2层剖分点。Level 3选取F,E,G,H,K,J作为第3层剖分点。Final这样，我们就构建出了整个树的R&C Trees。Randomizedn决定R&C Trees深度的关键因素在于选择Next(i)中元素的方法.n一个比较好的方法是,随机选择!n如果等概率的选择Next(i)中的元素,可以证明这样的深度下界是(ln2N).问题并没有完全解决.n必须采取更为合理的随机策略.Randomizedn对于Rake节点,仍然采取等概率的方法选择.n对于Compress节点i,j表示Next(i)中的元素,令Weight(j)表示j剩余的子孙个数(包括j).n现在以Weight(j)作为加权,令S表示所有Weight(k)的和,则j有Weight(j)/S的概率被选择.Randomized:Master Theoremn可以证明,通过这样合理的改造,一个N的点的树可以被映射到一个期望深度为O(lnN)的R&C Trees.n基于这种思想,RP-Trees被提出.事实上,这就是Treap7通过R&C思想在任意树的拓展.n通过这种思想,我们可以得到均摊,甚至是严格深度O(lnN)的算法.小结n相比较于Path Decomposing和Divide&Conquer,Rake&Compress具有思想简单,常数小,实现复杂度低等特点.nR&C思想最大的特点是,利用这种思想可以很方便地将各种平衡二叉树技巧拓展到任意树形态上去.n为Dynamic Trees Problem 注入了新的血液.Part III.ApplicationsApplications n网络优化n最大流4 n最小费用流nn动态算法n动态连通性3,8n动态最小生成森林3nn最大流问题n最大流问题是非常经典的图论问题.经典的解决算法有最短增广路,预流推进.n通过改造最短增广路算法并应用动态树,可以得到O(NMlnN)的算法.N为点数,M为边数.n经典的最短增广路算法是通过BFS,每次在残余网络中找到一条最短S-T增广路并进行增广.最大流问题n引理:只需要O(NM)次增广.n证明:考察当前最短路的必要边集An所有ST的最短路全部由A中元素组成n|A|Mn每次都找到最短增广路,增广后A中元素必有一条边被删除(残余量为0).最大流问题n在ST最短路长度被提高之前不可能有边从A外加入到A内.nO(M)次增广后ST增广路长度必被提高.因此最多执行O(NM)次增广.n综上所述,引理得证.最大流问题n原始算法的时间复杂度为O(NM2).n令D(x)表示x到T的最短增广路下界.n对于所有残余量0的边uv,满足D(u)D(v)+1.n如果D(u)=D(v)+1,则称uv为有效边.n引理:全部由有效边组成的到T的路径一定是最短路径!最大流问题n每次在残余网络中沿着让D(x)递减的有效边前进.并不断修正(抬高)D(x).可以证明,该优化方法将寻找最短增广路的时间降为均摊O(N),所以需要的总时间降为O(N2M).n那么,能不能再次改进呢?最大流问题n一个想法是,维护有效边子集的生成森林.n如果S,T不在同一棵树内.令R为S所在树内D值最小的点(最低点),如果存在有效边RR,则执行Link(R,R).n如果此时R没有连出的有效边,显然D(R)是可以改进的.即这个下界是松的,此时我们抬高D(R).并更新有效边集.n当S和T都在同一棵树内时(即最低点为T),对树中的S-T路径进行增广.n每次所增广的路都是最短增广路.Residual networkD(x)x黑色边为有效边Residual networkD(x)x红色边为生成森林边Residual networkD(x)x找到一条由当前最低点A连出的有效边AE执行Link(A,E)Residual networkD(x)x找到最短增广路SAET增广量为8Residual networkD(x)x增广后AE不再有效Cut(A,E)Residual networkD(x)x尝试找到从最低点A连出的边，失败抬高D(A)并更新有效边Residual networkD(x)x找到一条由当前最低点S连出的有效边SC执行Link(S,C)Residual networkD(x)x找到一条由当前最低点C连出的有效边CE执行Link(C,E)Residual networkD(x)x找到最短增广路SCET增广量为3Residual networkD(x)x最大流问题n由于流程和最短增广路一样,故正确性显然.n现在考察其时间复杂度.最大流问题n前面已经说过,最多执行O(NM)次增广,所以在增广上的总时间为O(NMlnN).n因为0D(x)N,所以最多执行O(NM)次Cut操作.因此花费在Cut上的总时间为O(NMlnN).n又因为在一个点的D值被抬高之前,每个点均摊被Link过1次.综上所述,总的时间复杂度为O(NMlnN)时间复杂度算法名称原始算法原始算法引进D(x)动态树找最短路O(NM2)O(N2M)O(NMlnN)进行增广O(N2M)O(N2M)O(NMlnN)调整D(x)O(NM)O(NMlnN)总计O(NM2)O(N2M)O(NMlnN)时间复杂度算法名称原始算法引进引进D(x)动态树找最短路O(NM2)O(N2M)O(NMlnN)进行增广O(N2M)O(N2M)O(NMlnN)调整D(x)O(NM)O(NMlnN)总计O(NM2)O(N2M)O(NMlnN)时间复杂度算法名称原始算法引进D(x)动态树动态树找最短路O(NM2)O(N2M)O(NMlnN)进行增广O(N2M)O(N2M)O(NMlnN)调整D(x)O(NM)O(NMlnN)总计O(NM2)O(N2M)O(NMlnN)小结n使用动态树问题来优化算法的关键是解决瓶颈问题.n在本例中,我们将复杂度摊平到了每一种操作中.使得总的时间复杂度获得了质的飞跃.n其实,不仅仅是最大流算法,网络优化的很多算法都可以使用动态树来优化,并得到很好的理论复杂度.总结n研究动态树问题是非常有价值的,研究它可以加深对图论的理解.其解决方法使用的技巧亦非常具有启发性.更何况其问题本身又是那么美妙!n对于一个复杂的经典问题,我们不应仅仅满足于“知道”或者“学会”,更需要融会贯通.这对我们无论是日常学习或者竞赛活动都是非常有好处的.Referencesn1 Daniel D.Sleator,Robert Endre Tarjan,A data structure for dynamic trees,Journal of Computer and System Sciences,v.26 n.3,p.362-391,June 1983.n2 Daniel Dominic Sleator,Robert Endre Tarjan,Self-adjusting binary search trees,Journal of the ACM(JACM),v.32 n.3,p.652-686,July 1985n3 S.Alstrup,J.Holm,K.de Lichtenberg,and M.Thorup.Minimizing diameters of dynamic trees.In Proceedings of the 24th International Colloquium on Automata,Languages and Programming,pages 270-280,1997.n4 Ahujia,R.K.Dynamic Trees Implementations.Network Flows:Theory,Algorithms,and Applications,p.265-273.n5 R.Tarjan and R.Werneck.Self-adjusting top trees.In Proceedings of the Sixteenth Annual ACM-SIAM Symposium on Discrete Algorithms(SODA),2005.n6 Umut A.Acar,Guy E.Blelloch,Jorge L.Vittes,Separating Structure from Data in Dynamic Trees.n7 C.R.Aragon,R.G.Seidel,Randomized Search Trees,Proc.30th Ann.IEEE Symposium on Foundations of Computer Science,pp.540-545,October 1989.n8 周源,Dynamic Connectivity研究报告,CTSC2005作业-研究报告.Thank you!Questions are welcome!