Kasaraju演算法--強連通圖遍歷及其python實現

来源:https://www.cnblogs.com/thomson-fred/archive/2018/12/20/10151972.html
-Advertisement-
Play Games

在理解有向圖和強連通分量前必須理解與其對應的兩個概念,連通圖(無向圖)和連通分量。 連通圖的定義是:如果一個圖中的任何一個節點可以到達其他節點,那麼它就是連通的。 例如以下圖形: 這是最簡單的一個連通圖,即使它並不閉合。由於節點間的路徑是沒有方向的,符合從任意一個節點出發,都可以到達其他剩餘的節點這 ...


在理解有向圖和強連通分量前必須理解與其對應的兩個概念,連通圖(無向圖)和連通分量。

連通圖的定義是:如果一個圖中的任何一個節點可以到達其他節點,那麼它就是連通的。

例如以下圖形:

這是最簡單的一個連通圖,即使它並不閉合。由於節點間的路徑是沒有方向的,符合從任意一個節點出發,都可以到達其他剩餘的節點這一條件,那麼它就是連通圖了。

 

連通分量

 

顯然這也是一個圖,只不過是由三個子圖組成而已,但這並非一個連通圖。這三個子圖叫做這個圖的連通分量,連通分量的內部歸根還是一個連通圖。

 

有向圖:

在連通圖的基礎上增加了方向,兩個節點之間的路徑只能有單一的方向,即要麼從節點A連向節點B,要麼從節點B連向節點A。有向圖與連通圖(更準確來說是無向圖)最大的區別在於節點之間的路徑是否有方向。

有向圖也分兩種,一種是有環路的有向圖。另外一種是無環路的有向圖,即通常所說的有向無環圖DAG(Directed Acyclic Graph)。嚴格來說,第一種有環路的圖,如果任意一個節點都可以與其他節點形成環路,那麼它也是一個連通圖。

 

例如下麵的就為一個有向圖同時也是連通圖:

 

強連通分量

強連通分量SCCs(strongly connected components)是一種有向的連通圖。

如果一個圖的連通分量是它裡面所有節點到能夠彼此到達的最大子圖,那麼強連通分量SCCs就是一個有向圖中所有節點能夠彼此到達的子圖。

 

顯然由345組成的子圖是無法到達由012組成的子圖的。那麼012和345分別組成兩個強連通分量。

 

在實際的現實問題中,我們考慮問題可能就不會簡單地研究無向圖。例如地圖上的最短路徑規劃,ARP路由演算法等等,考慮的都是有向圖的問題。

如果有這樣一個需求,我們希望用最少的次數遍歷所有節點,怎麼處理呢?

 

時間效應問題,強連通分量間的時間問題。

如果有向圖的各個強連通分量中的元素個數相仿,那麼,它們內部分別進行遍歷的時間量級別是相等的,但實際情況是,這種情況很少發生。一般從一個強連通分量到另一個強連通分量。

 

正如上面的需求:如何用最少的次數遍歷整個有向圖的所有節點。假設我們將0、1、2組成子圖1,將3、4、5組成子圖,子圖1有一條指向子圖2的路徑。這時候,我們從子圖1的任意一點開始遍歷。假設我們從1開始遍歷,那麼遍歷的順序將會是1—2,那麼來到2的時候問題來了,是先走0的路徑還是走子圖1和子圖2之間的路徑去遍歷節點3呢?

如果我們先遍歷節點0,那麼我們遍歷完節點0之後,發現節點1已經遍歷過,就會返回節點2,再沿著子圖1和子圖2之間的路徑去遍歷子圖2。這看起來是挺合理的。

但問題是,如果是先遍歷節點3(也就是說先遍歷子圖2)呢?

假設沿著子圖1和子圖2的路徑去遍歷子圖2,那麼子圖2遍歷完後,子圖1還剩下節點0沒有被遍歷,這時候就會出現很為難的事情,因為之前遍歷的情況無法判斷哪些節點是沒有遍歷的,只能是原路返回,依次去從新遍歷,“發現”哪些節點是還沒去遍歷的。似乎上圖比較簡單,這種方法不會耗費太多的時間。但如果是節點2連接著(並指向)許多個強連通子圖的有向圖,這種“返回式”的遍歷將會是很費勁的一件事。

 

為瞭解決這個問題,Kosaraju演算法提出了它的解決方案。Kosaraju演算法的核心操作是將所有節點間的路徑都翻轉過來。下麵分析一下為什麼這種演算法有它的優勢。

還是拿上面的圖來講述。想象一下上面的有向圖中的所有節點間的路徑都翻轉過來了。讀者可以自己用一張紙簡單畫一下。就像下麵的圖:

 

這一次,我們還是以0、1、2組成子圖1,以3、4、5組成子圖2。所不同的是,這次遍歷的起始點從子圖1開始。

 

多強連通分量的有向圖

再來看一下這個多子圖的強連通圖,如果像上圖所示,從子圖1開始,就會像上文提到的那樣,遍歷到節點2,會出現多個去向的問題。而在還沒有遍歷完子圖1的前提下,從節點2過渡到子圖2/子圖3,再回溯的時候會引來較大的麻煩。通過Kosaraju演算法之後,從2節點出發的路徑都會變成指向2。此時,遍歷的起點還是從子圖1開始,由於子圖1沒有出路,就不會出現上面所說的問題。再遍歷完子圖1後,繼續遍歷子圖2、子圖3。而子圖2、子圖3的遍歷都是在強連通分量內部實現的。

 

 

演算法實現

鄰接集表示的有向圖

N={
    "a":{"b"},   #a
    "b":{"c"},   #b
    "c":{"a","d","g"},   #c
    "d":{"e"},   #d
    "e":{"f"},   #e
    "f":{"d"},   #f
    "g":{"h"},   #g
    "h":{"i"},   #h
    "i":{"g"}    #i
}

 

翻轉圖實現代碼:

def re_tr(G):
    GT = {}
    for u in G:
        for v in G[u]:
            # print(GT)
            if GT.get(v):
                GT[v].add(u)
            else:
                GT[v] = set()
                GT[v].add(u)

    return GT

 

深度遍歷演算法實現代碼:

#遞歸實現深度優先排序
def rec_dfs(G,s,S=None):
    if S is None:
        #S = set()    #集合存儲已經遍歷過的節點
        S = list()    #用列表可以更方便查看遍歷的次序,而用集合可以方便用difference求差集
    # S.add(s)
    S.append(s)
    print(S)
    for u in G[s]:
        if u in S:continue
        rec_dfs(G,u,S)

return S

 

在強連通圖內遍歷

#遍歷有向圖的強連通分量
def walk(G,start,S=set()):     #傳入的參數S,即上面的seen很關鍵,這避免了通過連通圖之間的路徑進行遍歷
    P,Q = dict(),set()      #list存放遍歷順序,set存放已經遍歷過的節點     
    P[start] = None
    Q.add(start)
    while Q:
        u = Q.pop()                      #選擇下一個遍歷節點(隨機性)
        for v in G[u].difference(P,S):         #返回差集
            Q.add(v)
            P[v] = u
    print(P)    
    return P

 

獲得強連通分量

#獲得各個強連通圖
def scc(G):
    GT = re_tr(G)
    sccs,seen = [],set()
    for u in rec_dfs(G,"a"):    #以a為起點
        if u in seen:continue
        C = walk(GT,u,seen)
        seen.update(C)
        sccs.append(C)
    return sccs

 

單元測試

print(scc(N))

結果:
{'a': None, 'c': 'a', 'b': 'c'}
{'d': None, 'f': 'd', 'e': 'f'}
{'g': None, 'i': 'g', 'h': 'i'}
[{'a': None, 'c': 'a', 'b': 'c'}, {'d': None, 'f': 'd', 'e': 'f'}, {'g': None, 'i': 'g', 'h': 'i'}]

 

這是本人學習過程所寫的第一篇關於圖的演算法文章,供大家一起學習討論。其中難免會有錯誤。如有錯誤之處,請各位指出,萬分感謝!

 


您的分享是我們最大的動力!

-Advertisement-
Play Games
更多相關文章
  • 運行下麵的代碼: 出錯: 這是因為換成了Python3.6.5,Python3.6.5已經刪除了has_key()方法,改成了下麵的寫法: 就可以正常運行了。 ...
  • 1- 什麼是棧? 一個棧是一個項的有序集合。添加項和移除項都在同一端,這一端被稱為‘棧頂’。另一端被稱為‘棧底’。 棧使用的是後進先出原則即‘LIFO’原則,也就是說最新添加的項在移除時是第一個被移除的。在日常生活中有很多例子比如說在餐廳中有一堆餐盤,我們拿走的是最頂上的一個,排在我們後面的人將拿走 ...
  • 1、描述符的定義 描述符是與特定屬性互相綁定的一種協議,通過方法被觸發修改屬性,這些方法包括__get__(),__set__(),__delete__().將這些方法定義在類中,即可實現描述符 2、屬性與__dict__ Python中類有屬於自己的字典屬性,經過類的實例化的對象也同樣有自己的字典 ...
  • 1.反射 主要是用到了4個函數( 用的最多的就是getattr()和 hasattr() ): getattr() 從xxx對象中獲取到xxx屬性值 hasattr() 判斷xxx對象中是否有xxx屬性值delattr() 從xxx對象中刪除xxx屬性setattr() 設置xxx對象中的xxx屬性 ...
  • https://docs.scipy.org/doc/scipy/reference/interpolate.html#module-scipy.interpolate https://stackoverflow.com/questions/31464345/fitting-a-closed-cur ...
  • python爬蟲+數據可視化項目(一) 爬取目標:中國天氣網(起始url:http://www.weather.com.cn/textFC/hb.shtml#) 爬取內容:全國實時溫度最低的十個城市氣溫排行榜 使用工具:requests庫實現發送請求、獲取響應。 beautifulsoup實現數據解 ...
  • 例題import lxml.html test_data = """ <div> <ul> <li class="item-0"><a href="link1.html" id="places_neighbours__row">9,596,960first item</a></li> <li cla ...
  • sql註入後可以通過該資料庫獲取所有表的欄位信息 1. COLLATIONS表 提供有關每個字元集的排序規則的信息。 COLLATIONS表包含以下列: COLLATION_NAME 排序規則名稱。 CHARACTER_SET_NAME 與排序規則關聯的字元集的名稱。 ID 排序規則ID。 IS_D... ...
一周排行
    -Advertisement-
    Play Games
  • 移動開發(一):使用.NET MAUI開發第一個安卓APP 對於工作多年的C#程式員來說,近來想嘗試開發一款安卓APP,考慮了很久最終選擇使用.NET MAUI這個微軟官方的框架來嘗試體驗開發安卓APP,畢竟是使用Visual Studio開發工具,使用起來也比較的順手,結合微軟官方的教程進行了安卓 ...
  • 前言 QuestPDF 是一個開源 .NET 庫,用於生成 PDF 文檔。使用了C# Fluent API方式可簡化開發、減少錯誤並提高工作效率。利用它可以輕鬆生成 PDF 報告、發票、導出文件等。 項目介紹 QuestPDF 是一個革命性的開源 .NET 庫,它徹底改變了我們生成 PDF 文檔的方 ...
  • 項目地址 項目後端地址: https://github.com/ZyPLJ/ZYTteeHole 項目前端頁面地址: ZyPLJ/TreeHoleVue (github.com) https://github.com/ZyPLJ/TreeHoleVue 目前項目測試訪問地址: http://tree ...
  • 話不多說,直接開乾 一.下載 1.官方鏈接下載: https://www.microsoft.com/zh-cn/sql-server/sql-server-downloads 2.在下載目錄中找到下麵這個小的安裝包 SQL2022-SSEI-Dev.exe,運行開始下載SQL server; 二. ...
  • 前言 隨著物聯網(IoT)技術的迅猛發展,MQTT(消息隊列遙測傳輸)協議憑藉其輕量級和高效性,已成為眾多物聯網應用的首選通信標準。 MQTTnet 作為一個高性能的 .NET 開源庫,為 .NET 平臺上的 MQTT 客戶端與伺服器開發提供了強大的支持。 本文將全面介紹 MQTTnet 的核心功能 ...
  • Serilog支持多種接收器用於日誌存儲,增強器用於添加屬性,LogContext管理動態屬性,支持多種輸出格式包括純文本、JSON及ExpressionTemplate。還提供了自定義格式化選項,適用於不同需求。 ...
  • 目錄簡介獲取 HTML 文檔解析 HTML 文檔測試參考文章 簡介 動態內容網站使用 JavaScript 腳本動態檢索和渲染數據,爬取信息時需要模擬瀏覽器行為,否則獲取到的源碼基本是空的。 本文使用的爬取步驟如下: 使用 Selenium 獲取渲染後的 HTML 文檔 使用 HtmlAgility ...
  • 1.前言 什麼是熱更新 游戲或者軟體更新時,無需重新下載客戶端進行安裝,而是在應用程式啟動的情況下,在內部進行資源或者代碼更新 Unity目前常用熱更新解決方案 HybridCLR,Xlua,ILRuntime等 Unity目前常用資源管理解決方案 AssetBundles,Addressable, ...
  • 本文章主要是在C# ASP.NET Core Web API框架實現向手機發送驗證碼簡訊功能。這裡我選擇是一個互億無線簡訊驗證碼平臺,其實像阿裡雲,騰訊雲上面也可以。 首先我們先去 互億無線 https://www.ihuyi.com/api/sms.html 去註冊一個賬號 註冊完成賬號後,它會送 ...
  • 通過以下方式可以高效,並保證數據同步的可靠性 1.API設計 使用RESTful設計,確保API端點明確,並使用適當的HTTP方法(如POST用於創建,PUT用於更新)。 設計清晰的請求和響應模型,以確保客戶端能夠理解預期格式。 2.數據驗證 在伺服器端進行嚴格的數據驗證,確保接收到的數據符合預期格 ...