Part.5【马拉车&扩展KMP】

Manacher（马拉车）是一种求最长回文串的线性算法，复杂度O(n)。网上对其介绍的资料已经挺多了的，请善用搜索引擎。

而扩展KMP说白了就是是求模式串和主串的每一个后缀的最长公共前缀【KMP更像是一个自动机】

题目：

POJ 1159: Palindrome

求原字符串最少增加几个字符后可变成回文串，相当于求最长回文子序列的长度。

解法：直接求串S和反转串Sr的最长公共子序列。

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 10009
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int n, dp[2][5001];
char s[5001];

int main()
{
    scanf("%d%s", &n, s);
    rep(i, 1, n) rep(j, 1, n)
    {
        dp[i%2][j] = max(dp[(i-1)%2][j], dp[i%2][j-1]);
        if (s[i-1]==s[n-j] && dp[(i-1)%2][j-1]+1>dp[i%2][j]) dp[i%2][j] = dp[(i-1)%2][j-1]+1;
    }
    printf("%d\n", n-dp[n%2][n]);
    return 0;
}

HDU 3068: 最长回文

模板题。

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 234567
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int n, l, k[maxn], o, p, ans;
char c[maxn], s[maxn];

int main()
{
    while (scanf("%s", c) == 1)
    {
        l = strlen(c);
        rep(i, 0, l-1) s[i*2+1] = ‘$‘, s[i*2+2] = c[i];
        l = l*2+1; s[l]=s[l+1] = ‘$‘; s[0] = ‘#‘;
        o=p=ans=0;
        rep(i, 0, l)
        {
            k[i] = i<p ? min(k[2*o-i], p-i) : 1;
            while (s[i+k[i]] == s[i-k[i]]) k[i]++;
            if (i+k[i]>p) p = i+k[i], o = i;
            ans = max(ans, k[i]-1);
        }
        printf("%d\n", ans);
    }
    return 0;
}

POJ 3974: Palindrome

同模板题。

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 2000009
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int t, n, l, k[maxn], o, p, ans;
char c[maxn], s[maxn];

int main()
{
    while (++t)
    {
        scanf("%s", c);
        if (c[0]==‘E‘) break; l = strlen(c);
        rep(i, 0, l-1) s[i*2+1]=‘$‘, s[i*2+2]=c[i];
        l=l*2+1; s[0]=‘#‘; s[l]=s[l+1]=‘$‘;
        ans=o=p=0;
        rep(i, 1, l)
        {
            k[i] = i<p ? min(p-i, k[o*2-i]) : 1;
            while (s[i+k[i]] == s[i-k[i]]) k[i]++;
            if (i+k[i]>p) p=i+k[i], o=i;
            ans = max(ans, k[i]-1);
        }
        printf("Case %d: %d\n", t, ans);
    }
    return 0;
}

BZOJ 2342: 双倍回文

求最长双倍回文。

我实在是太弱了，于是写了个暴力+最优性剪枝，成功地用O(n^2)的算法AC！【其实正解是要O(nlogn)的算法的】

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 1000009
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int t, n, l, k[maxn], o, p, ans, b[maxn];
char c[maxn], s[maxn];

int main()
{
    scanf("%d%s", &l, c);
    rep(i, 0, l-1) s[i*2+1]=‘$‘, s[i*2+2]=c[i];
    s[0]=‘#‘; s[l*2+1]=‘$‘;
    ans=o=p=0;
    rep(i, 1, l*2+1)
    {
        k[i] = i<p ? min(p-i, k[o*2-i]) : 1;
        while (s[i+k[i]] == s[i-k[i]]) k[i]++;
        if (i+k[i]>p) p=i+k[i], o=i;
    }
    rep(i, 1, l) b[i+1]=k[i*2-1]/2;
    rep(i, 1, l)
        down(j, b[i]/2, ans+1) if (b[i-j]>=j && b[i+j]>=j) ans = j;
    printf("%d\n", ans*4);
    return 0;
}

BZOJ 2565: 最长双回文串

求最长双回文串。

对于每个字符，利用单调性质求出往左往右所能找到的最长的回文串。

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 200009
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int t, n, l, k[maxn], o, p, ans, lt[maxn], rt[maxn];
char c[maxn], s[maxn];

int main()
{
    scanf("%s", c); l=strlen(c);
    rep(i, 0, l-1) s[i*2+1]=‘$‘, s[i*2+2]=c[i];
    s[0]=‘#‘; s[2*l+1]=‘$‘;
    o = p = 0;
    rep(i, 1, l*2+1)
    {
        k[i] = i<p ? min(p-i, k[2*o-i]) : 1;
        while (s[i+k[i]] == s[i-k[i]]) k[i]++;
        if (i+k[i]>p) p = i+k[i], o = i;
    }
    o = 1; rep(i, 2, l*2+1)
    {
        while (o+k[o]<=i) o++;
        lt[i] = i-o;
    }
    o = l*2+1; down(i, l*2, 1)
    {
        while (i<=o-k[o]) o--;
        rt[i] = o-i;
    }
    rep(i, 1, l*2+1) if (s[i]==‘$‘ && lt[i] && rt[i] && ans < lt[i]+rt[i]) ans = lt[i]+rt[i];
    printf("%d\n", ans);
    return 0;
}

HDU 3613: Best Reward

这题貌似EXKMP和Manacher都可以做。

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 500009
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int tt, l, n[maxn], ex[maxn], v[maxn], t[26];
char c[maxn], s[maxn];

int main()
{
    scanf("%d", &tt);
    while (tt--)
    {
        rep(i, 0, 25) scanf("%d", &t[i]);
        scanf("%s", s); int l = strlen(s);
        rep(i, 0, l-1) c[i] = s[l-i-1];
        int k, p, a; rep(i, 1, l-1) v[i] = 0;

        rep(i, 0, l-1) ex[i] = n[i] = 0;
        n[1] = 0; while (c[n[1]] == c[n[1]+1]) n[1]++;
        n[0] = l; k = 1; p = n[1];
        rep(i, 2, l-1)
        {
            n[i] = i<p ? min(n[i-k], l-i) : 0;
            while (n[i]+i<l && c[n[i]] == c[n[i]+i]) n[i]++;
            if (i+n[i]>p) p = i+n[i], k = i;
        }
        k = p = 0;
        rep(i, 0, l-1)
        {
            ex[i] = i<p ? min(n[i-k], p-i) : 0;
            while (ex[i]+i<l && c[ex[i]] == s[ex[i]+i]) ex[i]++;
            if (i+ex[i]>p) p = i+ex[i], k = i;
        }
        a = 0; rep(i, 1, l-1) { a+=t[s[l-i]-‘a‘]; if (ex[l-i]==i) v[i]+=a; }

        rep(i, 0, l-1) ex[i] = n[i] = 0;
        n[1] = 0; while (s[n[1]] == s[n[1]+1]) n[1]++;
        n[0] = l; k = 1; p = n[1];
        rep(i, 2, l-1)
        {
            n[i] = i<p ? min(n[i-k], p-i) : 0;
            while (n[i]+i<l && s[n[i]] == s[n[i]+i]) n[i]++;
            if (i+n[i]>p) p = i+n[i], k = i;
        }
        k = p = 0;
        rep(i, 0, l-1)
        {
            ex[i] = i<p ? min(n[i-k], p-i) : 0;
            while (ex[i]+i<l && s[ex[i]] == c[ex[i]+i]) ex[i]++;
            if (i+ex[i]>p) p = i+ex[i], k = i;
        }
        a = 0; rep(i, 1, l-1) { a+=t[c[l-i]-‘a‘]; if(ex[l-i]==i) v[l-i]+=a; }

        int ans = 0; rep(i, 1, l-1) ans = max(ans, v[i]); printf("%d\n", ans);
    }
    return 0;
}

HDU 4333: Revolving Digits

扩展KMP的实例。。。

【代码莫名写WA了】

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 300009
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int tt, n[maxn];
char s[maxn];

int main()
{
    scanf("%d", &tt); int t = 0;
    while (t++<tt)
    {
        scanf("%s", s); int l = strlen(s);
        rep(i, 0, l-1) s[l+i] = s[i];
        int k, p;
        n[0] = l*2; n[1] = 0; while (s[n[1]] == s[n[1]+1]) n[1]++; p = n[k=1]+1;
        rep(i, 2, l-1)
        {
            n[i] = i<p ? min(p-i, n[i-k]) : 0;
            while (s[n[i]] == s[n[i]+i]) n[i]++;
            if (i+n[i]>p) p = i+n[i], k = i;
        }
        int q=0, w=0, e=0;
        rep(i, 0, l-1) if (n[i]>=l) w++; else if (s[n[i]] < s[i+n[i]]) e++; else q++;
        printf("Case %d: %d %d %d\n", t, q, w, e);
    }
    return 0;
}

URAL 1297: Palindrome

第三道模板题。

#include <cstdlib>
#include <cstdio>
#include <fstream>
#include <iostream>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <queue>

#define rep(i, l, r) for(int i = l; i <= r; i++)
#define down(i, l, r) for(int i = l; i >= r; i--)
#define ll long long
#define maxn 2009
#define MAX 1<<30
#define Q 1000000007

using namespace std;

int t, n, l, k[maxn], o, p, ans, b[maxn];
char c[maxn], s[maxn];

int main()
{
    scanf("%s", c); l=strlen(c);
    rep(i, 0, l-1) s[i*2+1]=‘$‘, s[i*2+2]=c[i];
    s[0]=‘#‘; s[2*l+1]=‘$‘;
    ans = o = p = 0;
    rep(i, 1, l*2+1)
    {
        k[i] = i<p ? min(p-i, k[2*o-i]) : 1;
        while (s[i+k[i]] == s[i-k[i]]) k[i]++;
        if (i+k[i]>p) p = i+k[i], o = i;
        if (k[ans]<k[i]) ans = i;
    }
    for(int i = ans-k[ans]+2; i <= ans+k[ans]-1; i+=2) printf("%c", s[i]);
    return 0;
}

总之，Manacher、KMP和ExKMP所利用的就是之前已经计算过的信息，以减少无意义的匹配。

Link

时间： 2024-03-12 15:00:50

Part.5【马拉车&扩展KMP】的相关文章

KMP、扩展KMP、MANACHER

一.KMP 作用:用于在一个文本串S内查找一个模式串P出现的位置 string str = "bacbababadababacambabacaddababacasdsd"; string ptr = "ababaca"; 如上图,可得在第10与26处包含ptr数组: 暴力做法:暴力for,碰到不一样的直接返回,从后一个开始继续for,最差能到O(n * m) KMP 做法: 主要的思路是跳,比如你一开始从上面例子里的bacbababadababacambabacad

hdu2328 Corporate Identity 扩展KMP

Beside other services, ACM helps companies to clearly state their "corporate identity", which includes company logo but also other signs, like trademarks. One of such companies is Internet Building Masters (IBM), which has recently asked ACM for

扩展KMP算法

扩展KMP,用于求s的后缀的最长前缀.用extand数组表示第i个后缀的最长前缀的字符个数. 注意几点:1.next数组是对T的 2.extand数组是对S的 3.应用:回文,重复串等代码如下: 1 #include<iostream> 2 #include<string> 3 #include<cstdio> 4 using namespace std; 5 const int MM=100005; //长度最大值 6 int next[MM],extand[M

扩展kmp——原创

扩展kmp LRH 所谓扩展kmp指的是与kmp相似的求辅助数组的原理,但是本身与kmp关系不大. 1.exkmp的用途:给定一个主串s和一个子串t,求出s中每一个后缀和子串t的最长公共前缀. 2.算法推导: 给定一个主串:S=aaaaaaaaaabaaa T=aaaaaaaaaaa (下标都是从零开始!!!) 第一步需要有两个辅助数组:extand[i]和next[i]: extand[i]:表示主串S以i开始的后缀与子串T

扩展KMP模板

扩展KMP: 给出模板串A和子串B,长度分别为lenA和lenB,要求在线性时间内,对于每个A[i](0 <= i < lenA),求出A[i..lenA-1]与B的最长公共前缀长度,记为ex[i](或者说,ex[i]为满足A[i..i + z - 1]==B[0 .. z - 1]的最大的z值). 扩展KMP可以用来解决很多字符串问题,如求一个字符串的最长回文子串和最长重复子串.[算法] 设next[i]为满足B[i..i + z - 1] == B[0..z - 1]的最

扩展KMP - HDU 4333 Revolving Digits

Revolving Digits Problem's Link: http://acm.hdu.edu.cn/showproblem.php?pid=4333 Mean: 给你一个字符串,你可以将该字符串的任意长度后缀截取下来然后接到最前面,让你统计所有新串中有多少种字典序小于.等于.大于原串. analyse: KMP的经典题. 首先我们将原串扩展成两倍,算一遍扩展KMP(自匹配),时间复杂度O(n). 这样一来,我们就得到了eKMP[i],eKMP[i]代表s[i...len-1]与s的最长

HDU 4300 （扩展KMP或KMP）

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4300 题意:说实话,看半天题目,愣是没看懂,百度题意才看懂,大概意思就是:第一个串串为匹配串,第i个代表字母顺序中的第i个,比如第一个是q,那么,q就代表a. 第二个串是密文+明文,密文全部都有,但明文可能不全,输出完整的密文+明文. 题解:有两种做法,第一种,用扩展KMP: 1.先把s2砍半,前面一半必定为密文,后面一半可能是密文+明文. 2.把前面的一半转化成明文. 3.用后面的和前面的找最长公

浅谈Manacher算法与扩展KMP之间的联系

首先,在谈到Manacher算法之前,我们先来看一个小问题:给定一个字符串S,求该字符串的最长回文子串的长度.对于该问题的求解,网上解法颇多,时间复杂度也不尽相同,这里列述几种常见的解法. 解法一通过枚举S的子串,然后判断该子串是否为回文,由于S的子串个数大约为,加上每次判断需要的时间,所以总的时间复杂度为,空间复杂度为. bool check(string &S, int left, int right) { while (left < right && S[left]

KMP与扩展KMP

原文转自:http://www.cppblog.com/MatoNo1/archive/2011/04/17/144390.aspx KMP:给出两个字符串A(称为模板串)和B(称为子串),长度分别为lenA和lenB,要求在线性时间内,对于每个A[i] (0<=i<lenA),求出A[i]往前和B的前缀匹配的最大匹配长度,记为ex[i](或者说,ex[i]为满足A[i- z+1..i]==B[0..z-1]的最大的z值).KMP的主要目的是求B是不是A的子串,以及若是,B在A中所有出现的位置