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