bzoj 3620: 似乎在梦中见过的样子 kmp暴力

       这道题目感觉上好像是可以用SA过掉的。。但是看一下N15000然后Statue里面一个个都6000+ms,翻了一下题解发现是kmp暴力。

       打了一半发现不会然后回去看题目,发现位置相同但是结构不懂的子串只算一次!!!这大概就是为什么不能用SA的原因吧。。不过这样kmp就好写多了。

      首先枚举头结点i,同kmp构造S[i,n]的失配函数,然后我们发现对于一个右端点j能更新答案当且仅当存在一个节点x,是j的失配指针上的某一个祖先(就是可以从j通过失配指针走到x),满足k<(i+j)/2且x-i+1>=k,也就是可以以[i,x]为A,[x+1,j-(x-i)]为B(+1-1的什么没有细考错了请指出)。于是可以用一个辅助数组g[j]表示满足x-i+1>=k的x中最小的x,当然如果fail[j]-i+1<k那么g[j]=j(fail[]为失配指针)。

AC代码如下:

#include<iostream>
#include<cstdio>
#include<cstring>
#define N 20005
using namespace std;

int n,m,f[N],g[N]; char ch[N];
int main(){
	scanf("%s",ch+1); n=strlen(ch+1);
	scanf("%d",&m); int i,j,k,ans=0;
	for (i=1; i<=n-(m<<1); i++){
		f[i]=k=i-1; g[i]=i;
		for (j=i+1; j<=n; j++){
			while (k!=i-1 && ch[j]!=ch[k+1]) k=f[k];
			if (ch[j]==ch[k+1]) k++; f[j]=k;
			g[j]=(k<i+m-1)?j:g[k];
			if (g[j]<((i+j)>>1)) ans++;
		}
	}
	printf("%d\n",ans);
	return 0;
}

by lych

2016.3.18

    原文作者:KMP算法
    原文地址: https://blog.csdn.net/lych_cys/article/details/50926124
    本文转自网络文章,转载此文章仅为分享知识,如有侵权,请联系博主进行删除。
点赞