說到樹鏈剖分,其實故事還挺多的。我記得在高中的OI經(jīng)歷中,我曾無數(shù)次聽到這個名詞,各種省賽、邀請賽貌似都會考這個東西,那時我覺得樹鏈剖分深不可測,是我等蒟蒻不能理解的東西……然后我還記得,某年(好像是NOip2014?)有一題貌似也要用樹鏈剖分。總之是被虐的一B,然后現(xiàn)在覺得這個東西也沒什么難的……
好吧,言歸正傳。所謂樹鏈剖分其實就是樹的輕重鏈剖分,而樹鏈就是指從某一父節(jié)點一直到葉子節(jié)點的一條鏈。下面先介紹幾個概念:
1、重兒子:父節(jié)點的所有兒子中,子樹節(jié)點數(shù)目最多的兒子
2、重邊:連接父節(jié)點和兒子的邊
3、重鏈:從父節(jié)點往下所有重兒子連在一起的鏈
有了這些定義,我們就能顯而易見的發(fā)現(xiàn):任何節(jié)點都只屬于一條樹鏈。然后按照這樣把一棵樹剖分成這樣一條一條的樹鏈,樹鏈剖分就完成了,如圖:
其中這棵樹被剖分成了紅、綠、藍三條鏈,但是要注意,圈出來的兩條邊其實是不屬于樹鏈的。是的完成了這個,樹鏈剖分就已經(jīng)做完了。其實剖分很容易理解,但是問題是這樣做了有什么用呢?這個坑我們后面填~接下來先說說如何實現(xiàn)它。
第一步當(dāng)然是對于每個非葉子節(jié)點都要選取一條重邊以及重兒子,這就是重鏈的延伸方向。如你所想,用dfs很容易實現(xiàn),依次遍歷該節(jié)點的子樹,然后挑選可達深度最大的點作為重兒子,具體代碼如下:
inline void dfs1(int u, int f, int d) //u為當(dāng)前節(jié)點,f為父親,d為當(dāng)前深度{ dep[u] = d; //先記下深度 size[u] = 1; //size[]表示子樹的節(jié)點數(shù)目 son[u] = 0; //son[]表示重兒子 fa[u] = f; efor(i,u) //遍歷子樹 { int ff = g[i].y; if (ff == f) continue; dfs1(ff, u, d + 1); size[u] += size[ff]; if (size[son[u]] < size[ff]) son[u] = ff; //選取節(jié)點數(shù)最多的點作為重兒子 }}標(biāo)記完重兒子后,我們要把每條鏈連起來,真正形成鏈。那么,我們用什么表示一條鏈呢?我們知道,樹狀數(shù)組、線段樹、平衡樹都是能支持在一段區(qū)間里維護性質(zhì),若我們能把一條鏈上的節(jié)點弄在一段連續(xù)的區(qū)間中,那我們就容易處理這些問題了。于是便想到了對每個節(jié)點進行重(chong)標(biāo)號,一條鏈上的點的新編號一定是連續(xù)的,具體的還是利用dfs來實現(xiàn),代碼:
inline void dfs2(int u, int tp) //tp表示當(dāng)前樹鏈的鏈頭,其實這個操作應(yīng)該叫l(wèi)abel{ top[u] = tp; //top[]表示節(jié)點i所屬樹鏈的鏈頭,方便查詢和修改 id[u] = ++num; //id[]表示一個節(jié)點新的id(編號) Rank[id[u]]=u; //Rank[]表示新編號i對應(yīng)原來的編號,即Rank[id[i]]=i if (son[u]) dfs2(son[u], tp); //先對同一條鏈上的重兒子標(biāo)號 efor(i,u) { int ff = g[i].y; if (ff == fa[u] || ff == son[u]) continue; dfs2(ff, ff); //對其他兒子標(biāo)號 }}這樣通過dfs1和dfs2,我們就完成了樹鏈剖分,現(xiàn)在開始填坑……
樹鏈剖分適合解決樹上的一些修改詢問以及求LCA之類的問題,同過把樹進行剖分,我們可以用很多數(shù)據(jù)結(jié)構(gòu)維護每一條鏈,然后再通過求LCA,我們可以解決對于樹上任意兩點的詢問問題。例如不斷的對邊權(quán)修改,求任意兩點之間的最大邊或最大值,這就可以用線段樹對每個樹鏈進行維護。另外,還值得一提的是樹鏈剖分后求LCA的方法,由于我們在dfs2的時候記錄了top[]數(shù)組,所以用樸素法就能很快的求出LCA,代碼:
inline int LCA(int u, int v){ int tp1 = top[u], tp2 = top[v]; //tp為鏈頭 while (tp1 != tp2) { if (dep[tp1] < dep[tp2]) //判斷深度 { swap(tp1, tp2); swap(u, v); } u = fa[tp1]; //每次向上走了tp1所以速度很快 tp1 = top[u]; } return tp1;}令我不滿意的是,我好像還是沒有說清楚具體如何進行修改和查詢,好吧,下次結(jié)合具體的題目再說把,再給自己挖一個坑……
最后說說我最不擅長的復(fù)雜度吧,由于是按照輕重鏈剖分的,所以有這樣的定理:
1、對于輕邊(u,v),size[v]<=size[u]/2
2、從根到某一點的路徑上,不超過logN條輕邊
3、重鏈總數(shù)不超過logN
因此剖分的時間復(fù)雜度為O(N);用數(shù)據(jù)結(jié)構(gòu)去維護每一條鏈時,每次的時間復(fù)雜度為O(logN),當(dāng)然你若不是用高級數(shù)據(jù)結(jié)構(gòu)我也沒辦法咯=_=……,所以總的來說時間復(fù)雜度是O(NlogN)級別的。
最后,我又看了一下百度,我覺的百科說的還是更有道理,樹鏈剖分是用來維護樹的路徑信息的一種算法?數(shù)據(jù)結(jié)構(gòu)?我還是不來定義了……
對了,那個坑我會填的^_^~~~
新聞熱點
疑難解答