HDU 6153 KMP

Posted #WoNderlAnd#

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了HDU 6153 KMP相关的知识,希望对你有一定的参考价值。

HDU6153:http://acm.hdu.edu.cn/showproblem.php?pid=6153

Orz 果然字符串很神奇……(神奇到人哭出来啊TAT)果然算法理解的不够透彻的话,是无法灵活运用的QAQ

其实只要简单的改一点点KMP的地方……

KMP的next[i]存储的是模式串下标0至i处相同的前缀和后缀的最长长度,每次与目标串匹配的时候,会略过前后缀相同的部分以减少回溯,提高效率。

由于每次匹配都是从前向后匹配的,因此对于模式串每次的匹配子串都是[0……i],是前缀,但是根据题目的要求,我们要找的是s2的子串[i……n],是后缀。

其实只要把模式串reverse一下,就可以把后缀转化为前缀,然后在进行处理就方便了。

对kmp理解的比较透彻的就不要听我瞎扯了……写的比较……emmmmm……不过还是很欢迎给窝这个连门槛都摸不到的小白提提建议的(^?^*)

可以用一个ll数组记录长度为i的前缀(也就是转化前的后缀)出现的次数(这里数组用crt表示)。(用ll……不要用int……会TLE)

因为next数组本身记录的就是相同子串的长度,因此可以在每个循环内部对匹配过的字符的下标进行计数

先码个代码:

while (i < len1)
    {
        if(j==-1||s1[i]==s2[j])
          {  i++; j++;   }
        else j = nextt[j];
        num[j]++;
        if (j == len2) j = nextt[j];
    }

 

下面拿题目的例子来解释一下出现次数的计数:

             b     a    b    a    b    a    b    a

next      -1    0    0    1    0    0    0    0

            a                                                        ——模式串目标串匹配过程

crt        0     0    0     0       匹配失败,j=next[j]=-1;这里crt数组不更新

crt        1     0    0     0

                   a

crt       1      1    0     0

                   a    b

crt       1      1     1                                         

                  a    b     a                  ——j==s2长度,但由于我们要在整个目标串中匹配,这里不能跳出,要继续匹配

crt       1     1    1     1                                      因此又加了  if(j==len2)  j=next[j];

                               a    b

crt       1     1     2    1

                               a    b    a

crt       1     1     2    2

              …………

                                           a   b    a

crt      1     1    3    3

可以发现,把crt计数放在每次j处理过之后,下标刚好就为目前子串的长度,所以就可以安心的记录啦;

前面有提到,kmp的算法为了减少回溯,省略了相同部分的匹配,因此我们每次记录的只是当前相同前缀长度的个数,从上面的例子不难发现,

除了第一次‘a‘以外的其他的‘a‘的匹配都被省略了,其他都是从’ab‘开始匹配的,每次匹配失败再次匹配时,省略了next[i]之前部分的计数,因此还需要进一步对计数数组进行处理。

for (int i = len2; i >0; i--)                
       crt[nextt[i]] += crt[i];

既然少了加回来就好了,匹配失败以后再次匹配时,都是从next[i]开始匹配的,因此,crt[next[i]]未计数的部分=crt[i],循环把len~1的部分加回来就OK了

(记得是从len倒着往前加)

下面看完整的代码。

#include<iostream>
#include<algorithm>
using namespace std;
const int MAX = 1e6;
const int mod = 1e9 + 7;
int nextt[MAX + 10];
char s1[MAX + 10];
char s2[MAX + 10];
long long num[MAX + 10];
int t,len1,len2;
void getnext()
{
    int p=0, k = -1;
    nextt[0] = -1;
    while (p < len2)
    {
        if (k == -1 || s2[p] == s2[k])
        {    
                    p++;k++;
                     nextt[p] =k;
                }
        else k = nextt[k];
    }
}
void kmp()
{
    getnext();
    int i = 0, j = 0;
    while (i < len1)
    {
        if(j==-1||s1[i]==s2[j])
        {i++; j++;}
        else j = nextt[j];
        num[j]++;
        if (j == len2) j = nextt[j];
    }
}
int main()
{
    cin >> t;
    while (t--)
    {
        memset(num,0, sizeof(num));
        cin >> s1;
        cin >> s2;
        len1 = strlen(s1);
        len2 = strlen(s2);
        reverse(s1, s1 + len1);
        reverse(s2, s2 + len2);
        kmp();
        long long ans = 0;
        for (int i = len2; i > 0; i--)
        {                
            num[nextt[i]] += num[i];
            ans = (ans + i*num[i]) % mod;
        }
            printf("%d\n", ans);
    }
    return 0;
}

有什么不对的地方或者解释的不好的地方还请大家指出来( ̄▽ ̄)"

以上是关于HDU 6153 KMP的主要内容,如果未能解决你的问题,请参考以下文章

HDU 6153 A Secret (KMP)

hdu6153KMP

HDU 6153 A Secret KMP

HDU 6153 拓展KMP (2017CCPC)

HDU 6153 A Secret(扩展kmp)

HDU 6153 A Secret(扩展KMP模板题)