算法模板——哈希单模板字符串匹配

实现功能：第一行输入模板串；第二行输入N；接下来N行每行一个字符串，将每个字符串中出现的模板串的起始位置找出

原理：字符串双值哈希啦啦啦，和KMP其实差不太多，但是字符串双值哈希绝对是个字符串题乱搞神器！！！

const pa=314159;pb=951413;
var
   i,j,k,l,m,n:longint;
   ap,bp:array[0..100000] of int64;
   ma,mb,va,vb:int64;
   s1,s2,s3:ansistring;
begin
     ap[0]:=1;bp[0]:=1;
     for i:=1 to 100000 do
         begin
              ap[i]:=(ap[i-1]*pa) mod pb;
              bp[i]:=(bp[i-1]*pb) mod pa;
         end;
     readln(s1);
     ma:=0;mb:=0;
     for i:=1 to length(s1) do
         begin
              ma:=(ma+ord(s1[i])*ap[i]) mod pb;
              mb:=(mb+ord(s1[i])*bp[i]) mod pa;
         end;
     readln(n);
     for j:=1 to n do
         begin
              readln(s2);
              if length(s2)<length(s1) then continue;
              va:=0;vb:=0;
              for i:=0 to length(s1)-1 do
                  begin
                       va:=(va+ord(s2[i])*ap[i]) mod pb;
                       vb:=(vb+ord(s2[i])*bp[i]) mod pa;
                  end;
              for i:=length(s1) to length(s2) do
                  begin
                       va:=(va-ord(s2[i-length(s1)])*ap[i-length(s1)]+ord(s2[i])*ap[i]) mod pb;
                       va:=(va+(abs(va) div pb+1)*pb) mod pb;
                       vb:=(vb-ord(s2[i-length(s1)])*bp[i-length(s1)]+ord(s2[i])*bp[i]) mod pa;
                       vb:=(vb+(abs(vb) div pa+1)*pa) mod pa;
                       if (((ma*ap[i-length(s1)]) mod pb)=va) and (((mb*bp[i-length(s1)]) mod pa)=vb) then write(i-length(s1)+1,' ');
                  end;
              writeln;
         end;
end.