hdu 5918 Sequence I （CCPC长春，KMP）

最新推荐文章于 2019-09-24 16:09:12 发布

原创最新推荐文章于 2019-09-24 16:09:12 发布 · 442 阅读

0 ·

CC 4.0 BY-SA版权

字符串——KMP 专栏收录该内容

8 篇文章

订阅专栏

本文解析了一道ACM竞赛题目，通过将序列按步长分组，并使用KMP算法进行模式匹配来解决给定两个序列及步长求匹配子串数量的问题。

题目：http://acm.split.hdu.edu.cn/showproblem.php?pid=5918
题意：

给你a,b两个序列和一个p ，求有多少个 q恰好满足 b1,b2,b3….bm 等于a[q],a[q+p],a[q+2p]……a[q+(m-1)p];

分析：

将a序列，每隔p位置分成一组，这样最多有p组，个数和是n
将每组和b序列跑kmp计算答案
这题暴力也能A

代码：

#include<bits/stdc++.h>
using namespace std;
const int N=1e6+5;
int s[N],t[N],slen,tlen,p;
int nex[N];
void getnex() {
    int j=0,k;
    k=nex[0]=-1;
    while(j<tlen) {
        if(k==-1||t[j]==t[k])nex[++j]=++k;
        else k=nex[k];
    }
}
vector<int>g[N];
int main() {
    int T;
    //freopen("f.txt","r",stdin);
    scanf("%d",&T);
    for(int cas=1; cas<=T; cas++) {
        scanf("%d%d%d",&slen,&tlen,&p);
        for(int i=0; i<slen; i++)scanf("%d",&s[i]);
        for(int i=0; i<tlen; i++)scanf("%d",&t[i]);
        for(int i=0; i<p; i++)g[i].clear();
        for(int i=0; i<slen; i++)g[i%p].push_back(s[i]);
        getnex();
        int ans=0;
        for(int k=0; k<p; k++) {
            int j=0;
            for(int i=0; i<g[k].size(); i++) {
                while(j&&g[k][i]!=t[j])j=nex[j];
                if(g[k][i]==t[j])j++;
                if(j==tlen) {
                    ans++;
                    j=nex[j];
                }
            }
        }
        printf("Case #%d: %d\n",cas,ans);
    }
    return 0;
}