扩展kmp(hdu 6153-加上注释的代码）

最新推荐文章于 2019-11-01 18:27:57 发布

我要去看得最远的地方

最新推荐文章于 2019-11-01 18:27:57 发布

阅读量146

点赞数

分类专栏： KMP

本文链接：https://blog.csdn.net/qq_42434171/article/details/83317396

版权

KMP 专栏收录该内容

6 篇文章 0 订阅

订阅专栏

://wenku.baidu.com/view/8e9ebefb0242a8956bece4b3.html//这个是扩展kmp的讲解讲得不错

例题（hdu 6153）

第一类

if((i-1)+nxt[i-a]<p)//解析此代码

//a......p
// i....
//i-a个已匹配(不包括第i个) 然后nxt[i-a]从t(i-a....lent) 与t[0....lent] 的最长公共长度假如为w

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long ll;
#define mod 1000000007 
ll extend[1000009];
ll nxt[1000009];
char t[1000009],s[1000009];
int lent,lens;
//nxt[i]表示t[i...m]与t的最长公共前缀长度 
//extend[i]=s[i....n]与t的最长公共前缀长度； 
void getnext()
{//此时t既是子串又是母串 
	nxt[0]=lent;//t的[0...lent]与t的最长公共在前缀长度 
	ll a,p;
	a=1;
	while(a<lent&&t[a]==t[a-1])a++;//匹配找nxt[1]; 
	nxt[1]=a-1;
	a=1;
	for(ll i=2;i<lent;i++)
	{// 第一种 
		p=a+nxt[a]-1;            //a......p
					 //   i....
					 //i-a个已匹配(不包括第i个) 然后nxt[i-a]从t(i-a....lent) 与t[0....lent] 的最长公共长度假如为w 
		if((i-1)+nxt[i-a]<p)//nxt[i-a]
		nxt[i]=nxt[i-a];
		else
		{//第二种 
			ll j=(p-i+1)>0?(p-i+1):0;
			while(i+j<lent&&t[i+j]==t[j]) j++;
			nxt[i]=j;
			a=i;//更新 
		}
	}
//核心一个起始位置a,一个最远匹配位置p,然后next和extend数组 
}
void exkmp()
{
	ll a,p;
	a=p=0;
	ll len=min(lens,lent);
	while(p<len&&t[p]==s[p])p++;
	extend[0]=p;
	for(ll i=1;i<lens;i++)
	{
		p=a+extend[a]-1;
		if((i-1)+nxt[i-a]<p)
		{
			extend[i]=nxt[i-a];
		}
		else
		{
			ll j=(p-i+1)>0?(p-i+1):0;
			while(j<lent&&i+j<lens&&s[i+j]==t[j])
			{
				j++;
			}
			extend[i]=j;
			a=i;
		}
	}
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
	scanf("%s%s",s,t);
	lent=strlen(t);
	lens=strlen(s);
	reverse(s,s+lens);
	reverse(t,t+lent);
	memset(nxt,0,sizeof(nxt));
	memset(extend,0,sizeof(extend));
	getnext();
	exkmp();
/*	ll sum=0;
	for(int i=0;i<lens;i++)
	{
		if(extend[i])
		sum=(sum%mod+(extend[i]%mod*(extend[i]+1)/2)%mod)%mod;
		//printf("%lld",extend[i]);
	}
	printf("%lld\n",sum);*/
}
	return 0;
	
}