HDU 1711 Number Sequence

大意是给两个序列，求第二个序列在第一个序列中的最小位置，无解输出-1。

朴素算法O(m*n)，m<=10000，n<=1000000.肯定超时

所以就要用到KMP算法，俗称看毛片算法(我是不喜欢这么叫的。。)

KMP算法是预处理出目标串的失配边，然后遇到不匹配的就顺着失配边走，直到找到匹配。

看起来很难理解，那末举个例子：

1<=i<=11

b[i]=(1 2 3 1 2 5 6 1 2 3 4)

那末它的失配数组是 

1<=i<=12

f[i]=(1 1 1 1 2 3 1 1 2 3 4 1)

意思是，当i这个位置不匹配时，就转到f[i]继续进行比较。

为什么多一个呢？因为匹配完之后要回到开头。

构造和比较很像，也是顺着失配边走

代码：

<span style="font-size:18px;">#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int f[1000020],b[1000020],a[1000010];
int j,n,m,t;
void kmp(int m)
{
f[1]=1;f[2]=1;
int j=0;
int i=0;
for(i=2;i<=m;++i)
{
j=f[i];
while(j>1&&b[i]!=b[j])j=f[j];//跟着失配边走，直到回到开头或找到匹配
if(b[i]==b[j])f[i+1]=j+1;
else f[i+1]=1;//如果相等就从下一个开始比较
}//构造失配边
j=1;
for(i=1;i<=n;++i)
{
while(j>1&&b[j]!=a[i])j=f[j];//跟着失配边走，直到回到开头或找到匹配
if(b[j]==a[i])++j;//如果相等就从下一个开始比较
if(j==m+1)
{
printf("%d\n",i-m+1);
return;
}
}
printf("%d\n",-1);
return;
}
int main()
{
//freopen("1711.in","r",stdin);
scanf("%d",&t);
for(int i=1;i<=t;++i)
{
memset(f,0,sizeof(f));
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
scanf("%d%d",&n,&m);
for(int j=1;j<=n;++j)
scanf("%d",&a[j]);
for(int j=1;j<=m;++j)
scanf("%d",&b[j]);
//if (kmp(m)==-1) printf("%d",-1);
kmp(m);
}
return 0;
}</span>