【距离GDOI：137天】 扩展KMP...字符串QAQ

上次和黄神两人一合计，干脆我学字符串他学图论，然后两人相互教...但以蒟蒻最近这状态来看，估计会到时候也教不了QAQ

扩展KMP大概是下面这个课件讲的这样

http://wenku.baidu.com/view/8e9ebefb0242a8956bece4b3.html

然后需要注意的是初始化问题：

对于next数组，根据定义next[0]=len  ，而因为k<>0，否则next[i]=next[i-k]不成立...所以next[1]也要预处理出来，才能继续求出剩余的next

对于extend数组，还是要预处理出ex[0]，理由同上。暴力预处理即可...

然后是P党的ansistring问题，有人说这东西会慢，而且不方便调试，蒟蒻干脆改成了char数组去模拟，还可以自己定义起点从0还是从1，我的代码是从0开始的，因为从1开始总是各种跪...估计是各种加加减减没处理好，干脆改成0了...反正以后滚去C++也是从0开始的...

然后就是课件里面讲的那样了..

主体代码如下...嗯有些变量名可能没声明...也没主程序...

 const maxn=;
 type
  extype=array[..maxn] of longint;
  strtype=array[..maxn] of char;
 var
  ex0,ex1:extype;
  s1,s2:strtype;
  next,sum:array[..maxn] of longint;

 procedure makenext(t:strtype;len:longint);
 var p,i,k:longint;
 begin
  k:=;
  next[]:=len;
  while (k+<len) and (t[k]=t[k+]) do inc(k);
  next[]:=k;
  k:=; i:=;
  while i<len do
   begin
    p:=k+next[k]-;
    next[i]:=min(next[i-k],max(p-i+,));
    while (i+next[i]<len) and (t[next[i]]=t[i+next[i]]) do inc(next[i]);
    if i+next[i]>k+next[k] then k:=i;
    inc(i);
   end;
 end;

 procedure makeext(t,s:strtype;var ex:extype;len:longint);
 var p,i,k:longint;
 begin
  makenext(t,len);
  k:=;
  while (k<len) and (s[k]=t[k]) do inc(k);
  ex[]:=k;
  k:=; i:=;
  while i<len do
   begin
    p:=k+ex[k]-;
    ex[i]:=min(next[i-k],max(p-i+,));
    while (i+ex[i]<len) and (t[ex[i]]=s[i+ex[i]]) do inc(ex[i]);
    if i+ex[i]>k+ex[k] then k:=i;
    inc(i);
   end;
 end;