主頁 > 後端開發 > 最長回文子串問題-Manacher演算法【建議收藏】

最長回文子串問題-Manacher演算法【建議收藏】

2021-09-27 09:18:52 後端開發

前面我們講過一個關于字串的演算法:KMP演算法,今天我們來講另外一個字串演算法:Manacher演算法,這個演算法是用于解決一個問題叫:最長回文子串

前期文章:KMP演算法

牛客網OJ鏈接

image-20210926152653053

說的簡單一點,給定一個字串,回傳的值是這個字串的最長回文子串的長度,顧名思義,即是回文串,也是子串,

文章目錄

  • 一、BF演算法
  • 二、Manacher演算法

一、BF演算法

那上圖的示例2為例:abc1234321ab,

最簡單的思路就是從左到右遍歷每一個字符,每來到一個字符位置,我們可以向左右兩邊進行擴展,分別比較左右兩邊的字符,如果相等,就繼續向兩邊擴展;如果不相等,就停止,計算以當前字符,向兩邊擴展出的長度,就是以當前字符為中心的回文子串,比如:

img

imgimg

就像上圖這樣,從左往右依次遍歷即可,

上面這種思路確實能夠解題,但是還有一個很重要的點,那就是假設給定的字串是偶數個字符,那么這種方式就會錯過一些回文子串的匹配,因為此時對于偶數個字符來說,對稱點是在中間兩個字符之間的,如下圖:

img

所以以每個字符為中心點,向兩邊擴展,還是存在著一定的問題,如何解決呢?

那就是將原字串進行處理,加工為一個含有特殊字符的字串,比如原字串為:123321,;加工后的字串為:#1#2#3#3#2#1#;

也就是說,在每個字符的中間,加入其它字符,這樣就能使一個偶數個字符的字串,轉換為奇數個字符的字串,這樣就可以遍歷,向兩邊擴展了,

問題:我們所加入的字符,必須是原字符中沒有的字符嗎?

這個問題留作大家思考,

public static int getLengthOfSubString(String str) {
    if (str == null) {
        return 0;
    }

    char[] generateStr = generateString(str);
    int length = generateStr.length;
    int max = 0; //答案
    
    for (int i = 0; i < length; i++) {
        int tmp = 1; //每個字符都能以自己本身的字符作為回文子串,所以初始值是1
        int radius = 1; //回文半徑,也就是以i位置為中心,半徑radius的范圍內

        while (i - radius >= 0 && i + radius < length) { //左右兩邊都在陣列的范圍內,回圈繼續
            if (generateStr[i - radius] == generateStr[i + radius]) {
                tmp += 2; //左右兩個字符相等的情況
                radius++; //回文半徑加1
            } else {
                break;
            }
        }

        max = Math.max(max, tmp); //判斷當前的tmp是否是最長的回文子串
    }
    return max / 2; //因為我們比較的處理后的字串,計算出的回文串要除以2.才是最終的答案
}

public static char[] generateString(String str) {
    char[] res = new char[str.length() * 2 + 1]; //原2倍長度,再加1
    int index = 0;
    for (int i = 0; i < res.length; i++) {
        //奇數位置放#,偶數位置放原字符
        res[i] = (i % 2) == 1? str.charAt(index++) : '#';
    }
    return res;
}

以上代碼就是BF演算法,暴力解,每個字符都需要遍歷一次,而每次字符都需要向外擴展,最壞情況下,就是向外擴展一直到整個字串結束,比如:1111111; 這種情況就是每個字符向外擴展,都會擴展很長,甚至是擴展至字串結束,所以這個BF演算法的時間復雜度是O(N2) ,

二、Manacher演算法

Manacher演算法也是在BF演算法的基礎之上,做了優化,所以大家看Manacher演算法之前,先理解BF暴力解的流程,

Manacher演算法引入了三個概念:

  • 當前回文子串的中心點 :C
  • 當前已經遍歷到最長回文子串的最右邊界下標:R
  • 回文半徑陣列;(用于存盤已經擴展完成的回文子串的半徑)

通過上面三個變數,我們就能解決這一難題了,話不多說,講述推導程序,整體分為兩個大步驟,C和R的初始值都是-1,也就是陣列最左邊的外面,

  1. 當i位置(當前遍歷的字符)不在R(最右邊界)內時

    img

    此時這種情況,我們只能向左右兩邊進行擴展,這個沒辦法,重要的是第2種情況,

  2. 當i位置(當前遍歷的字符)在R(最右邊界)內時

    image-20210926165520939

    以7為中心,向兩邊擴展出來的回文子串,就是橙色括號圈起來的范圍,此時的i就是在R邊界的里面,

    當我們以中心點為對稱點,做出i的對稱點,如下圖:

    img

    做出來的對稱點,我們就能得到這個點的下標,然后去回文半徑陣列里查這個下標對應的回文半徑,就能得到關于這個對稱點的回文子串,例如上圖中黑色虛線框中的值

    對于黑色虛線框的值,我們又可以分為三種小情況:

    • 黑色虛線框與以C中心點擴展的回文子串壓線

      壓線的情況,就是上圖中這種情況,黑色虛線框的左邊界與橙色線重合, 根據對稱性,因為黑色虛線框的值是回文子串,那么右邊以i為中心,也能擴展出回文子串,如下圖所示:

      img

      所以我們可以直接通過對稱點i得到已經完成匹配的回文子串,然后我們可以直接從i位置的已經計算好的回文子串外開始擴展,比如:左邊值7和右邊值1做比較,如果相等,當前回文半徑加1,然后繼續比較下一對字符,

    • 黑色虛線框的左邊界,超過了以C中心點擴展的回文子串的左邊界(超出):如下圖:

      img

      對稱點i,以它為中心對應的回文子串正如左邊的黑色虛線框所示:2,3,4,3,2,此時虛線框已經超出了橙色線的范圍,又因為橙色線范圍內是一個回文子串,所以我們可以推匯出當前i位置,至少有回文子串,就是(R-i)為半徑的范圍,即上圖右邊黑色虛線框內,

      此時我們只需要在此基礎之上,比較R右邊的值5 和 黑色虛線框左邊的2,看是否相等,若相等,則再次比較下一對字符,依次類推,

    • 黑色虛線框整體,都是在以C中心點擴展的回文子串的左半部分(即沒壓線,也沒超出):如下圖:

      img

      此時以i位置為中心,向左右兩邊擴展,就可以從黑色虛線框兩邊開始比較字符了,

    上面三種情況,都是由對稱點i得到關于該點的回文子串;再對稱到右邊i位置,以此為基礎,繼續向外擴展比較字符,那可能有同學就會疑惑,為什么就能從左邊對稱點i,就能推匯出右邊i位置的回文子串呢? 證明如下:

    img

上述所有,就是Manacher的推導程序,就是通過對稱,拿到C點左邊的對稱點,就能從回文半徑陣列中拿到該位置的回文子串,因此就能對應到C點右邊的回文子串,在此基礎之上進行字符比較,節省了一些已經比較過的字符的時間,

public static int manacherStr(String str) {
    if (str == null || str.length() < 1) {
        return 0;
    }

    char[] s = generateString(str);
    int length = s.length;
    int C = -1; //回文子串的中心點
    int R = -1; //最長回文子串的右邊界
    int[] pArr = new int[length]; //回文半徑陣列
    int max = 0; //答案
    for (int i = 0; i < length; i++) {
        //判斷i是否在R的范圍內,如果不在,選擇相應的對稱點,
        pArr[i] = i < R? Math.min(R - i, pArr[2 * C - i]) : 1; 

        //以下回圈,就是上面上面分析的3種小情況,也可以自己用if else陳述句
        while (i - pArr[i] >= 0 && i + pArr[i] < length) {
            if (s[i - pArr[i]] == s[i + pArr[i]]) { //左右兩邊的字符
                pArr[i]++; //回文半徑加1
            } else {
                break;
            }
        }

        //更新 新的回文子串的右邊界和 C中心點
        if (i + pArr[i] > R) {
            R = i + pArr[i];
            C = i;
        }
        max = Math.max(max, pArr[i]); //判斷是否是最長回文半徑
    }

    return max - 1; //最終的答案,與max的值,相差1
}

public static char[] generateString(String str) {
    char[] res = new char[str.length() * 2 + 1]; //原2倍長度,再加1
    int index = 0;
    for (int i = 0; i < res.length; i++) {
        //奇數位置放#,偶數位置放原字符
        res[i] = (i % 2) == 1? str.charAt(index++) : '#';
    }
    return res;
}

上述所有,就是Manacher演算法的全部,Manacher就是在BF演算法基礎之上,新加了回文半徑陣列,對于這個陣列來,可以解決很多關于字串的問題,所以很好的掌握這個演算法,對以后刷題有很大的幫助,

好啦,本期更新就到此結束啦!!!我們下期見!!!

img

轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/303303.html

標籤:java

上一篇:?每日演算法&面試題?一起肥學9??

下一篇:??蘇州程式大白一文教你學會微信小程式開發??《??記得收藏??》

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more