• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            C++分析研究  
            C++
            日歷
            <2011年9月>
            28293031123
            45678910
            11121314151617
            18192021222324
            2526272829301
            2345678
            統計
            • 隨筆 - 92
            • 文章 - 4
            • 評論 - 4
            • 引用 - 0

            導航

            常用鏈接

            留言簿

            隨筆檔案

            文章檔案

            搜索

            •  

            最新評論

            閱讀排行榜

            評論排行榜

             
              題意:給出一個字符串,求出最長回文字串。
             
               思路:一開始我直接上了后綴數組DC3的解法,然后MLE了。看了DISCUSS發現還有一種計算回文字串更加優越的算法,就是manacher算法。就去學習了一下,
             
               這個算法要解決的就是一個字符串中最長的回文子串有多長。這個算法可以在O(n)的時間復雜度內既線性時間復雜度的情況下,求出以每個字符為中心的最長回文有多長,
             
               這個算法有一個很巧妙的地方,它把奇數的回文串和偶數的回文串統一起來考慮了。這一點一直是在做回文串問題中時比較煩的地方。這個算法還有一個很好的地方就是充分利用了字符匹配的特殊性,避免了大量不必要的重復匹配。
             
               算法大致過程是這樣。先在每兩個相鄰字符中間插入一個分隔符,當然這個分隔符要在原串中沒有出現過。一般可以用'#'分隔。這樣就非常巧妙的將奇數長度回文串與偶數長度回文串統一起來考慮了(見下面的一個例子,回文串長度全為奇數了),然后用一個輔助數組P記錄以每個字符為中心的最長回文串的信息。P[id]記錄的是以字符str[id]為中心的最長回文串,當以str[id]為第一個字符,這個最長回文串向右延伸了P[id]個字符。
               原串: w aa bwsw f d
               新串: # w# a # a # b# w # s # w # f # d #
               輔助數組P: 1 2 1 2 3 2 1 2 1 2 1 4 1 2 1 2 1 2 1
               這里有一個很好的性質,P[id]-1就是該回文子串在原串中的長度(包括'#')。如果這里不是特別清楚,可以自己拿出紙來畫一畫,自己體會體會。當然這里可能每個人寫法不盡相同,不過我想大致思路應該是一樣的吧。
               好,我們繼續。現在的關鍵問題就在于怎么在O(n)時間復雜度內求出P數組了。只要把這個P數組求出來,最長回文子串就可以直接掃一遍得出來了sat答案
               由于這個算法是線性從前往后掃的。那么當我們準備求P[i]的時候,i以前的P[j]我們是已經得到了的。我們用mx記在i之前的回文串中,延伸至最右端的位置。同時用id這個變量記下取得這個最優mx時的id值。(注:為了防止字符比較的時候越界,我在這個加了'#'的字符串之前還加了另一個特殊字符'$',故我的新串下標是從1開始的)托福答案
               /*****************************************************************************************************************************************************************/
               CODE:
               #include <set>
               #include <map>
               #include <stack>
               #include <cmath>
               #include <queue>
               #include <cstdio>
               #include <string>
               #include <vector>
               #include <iomanip>
               #include <cstring>
               #include <iostream>
               #include <algorithm>
               #define Max 2505
               #define FI first
               #define SE second
               #define ll long long
               #define PI acos(-1.0)
               #define inf 0x3fffffff
               #define LL(x) ( x 《 1 )
               #define bug puts("here")
               #define PII pair<int,int>
               #define RR(x) ( x 《 1 | 1 )
               #define mp(a,b) make_pair(a,b)
               #define mem(a,b) memset(a,b,sizeof(a))
               #define REP(i,s,t) for( int i = ( s ) ; i <= ( t ) ; ++ i )
               #define N 2000055
               using namespace std;
               char s[N] ;
               char str[N] ;
               int rad[N] ;
               int manacher () {
               int len = strlen(s) ;
               int max = 0;
               str[0] = '$';
               str[1] = '#';
               int i = 0 ;
             
               for (; i < len; i++) {
               str[i * 2 + 2] = s[i];
               str[i * 2 + 3] = '#';
               }
               str[2 * len + 2] = 0;
               for (int i = 1; i < 2 * len + 2 ; i++) {
               rad[i] = 0;
               }
               int id = 0;
               for (i = 1; i < 2 * len + 2; i++) {
               if (max > i)
               rad[i] = min(rad[2 * id - i], rad[id] + id - i) ;
               else
               rad[i] = 1 ;
               while (str[i + rad[i]] == str[i - rad[i]])
               rad[i] ++ ;
               if (rad[i] + i > max) {
               max = rad[i] + i;
               id = i;
               }
               }
               int mx = 0;
               for (i = 1; i < 2 * len + 2 ; i++) {
               if (mx < rad[i] - 1)
               mx = rad[i] - 1;
               }
               return mx;
               }
               int main() {
               int ca = 0 ;
               while(scanf("%s",s) != EOF) {
               if(strcmp(s , "END") == 0)break ;
               printf("Case %d: " , ++ ca) ;
               printf("%d\n",manacher()) ;
               }
               return 0 ;
               }
            posted on 2013-11-17 12:19 HAOSOLA 閱讀(462) 評論(0)  編輯 收藏 引用
             
            Copyright © HAOSOLA Powered by: 博客園 模板提供:滬江博客
            PK10開獎 PK10開獎
            三上悠亚久久精品| 久久精品国产72国产精福利| 久久露脸国产精品| 久久综合九色综合网站| 成人久久免费网站| 97久久精品午夜一区二区| 久久国产免费| 亚洲午夜久久久久妓女影院| 久久久久久狠狠丁香| 日韩影院久久| 色综合久久天天综合| 久久精品国产99国产精品亚洲| 波多野结衣中文字幕久久 | 色综合久久无码五十路人妻| 久久九九兔免费精品6| 亚洲国产成人久久精品动漫| 日韩欧美亚洲综合久久影院Ds| 久久精品国产亚洲77777| 亚洲国产小视频精品久久久三级| www性久久久com| 久久久无码精品亚洲日韩京东传媒| av午夜福利一片免费看久久| 亚洲国产精品无码久久久不卡 | 国产精品内射久久久久欢欢| 亚洲va国产va天堂va久久| 久久久精品国产亚洲成人满18免费网站 | 久久国产精品成人影院| 久久婷婷五月综合国产尤物app | 亚洲精品成人网久久久久久| 国产无套内射久久久国产| 国产一级持黄大片99久久| 漂亮人妻被黑人久久精品| 久久国产色av免费看| 伊色综合久久之综合久久| 久久伊人中文无码| 欧美日韩精品久久久免费观看| 99久久www免费人成精品| 7国产欧美日韩综合天堂中文久久久久| 国产成人久久精品一区二区三区| 久久久久亚洲AV片无码下载蜜桃| 久久久久人妻精品一区|