BNUOJ 34990 Justice String (基于hash的LCP)

转自http://blog.csdn.net/w20810/article/details/48438047

题意：给定字符串A和B,在A里面找一个子串s，在s里面最多修改2个字符后使得s和B相同。求s的起始位置。

还是那种套路，不过这个地方的f1，f2，f3表示得挺好的。。

#include <iostream>#include <cstdio>#include <cstring>#include <algorithm>using namespace std;typedef unsigned long long ULL;const int maxn = 2e5+6;const int seed = 131;ULL H[maxn],X[maxn];char s1[maxn],s2[maxn];int Len1,Len2,L;void Init(){    H[L]=0;    for(int i=L-1;i>=0;i--)        H[i]=H[i+1]*seed+s1[i]-'a';    X[0]=1;    for(int i=1;i<L;i++)        X[i]=X[i-1]*seed;}inline ULL GetHash(int i,int L) //{    return H[i]-H[i+L]*X[L];}int LCP(int a,int b){    int ret=0,down=1,up=L-b,mid;    while(down<=up)    {        mid=(down+up)>>1;        if(GetHash(a,mid)==GetHash(b,mid))        {            down=mid+1;            if(mid>ret)                 ret=mid;        }        else            up=mid-1;    }    return ret;}int main(){    int ncase,i,j,f1,f2,f3,ans;    scanf("%d",&ncase);    for(int T=1;T<=ncase;T++)    {        scanf("%s%s",s1,s2);        Len1=strlen(s1);        Len2=strlen(s2);        L=Len1+Len2;        strcat(s1,s2);        Init();        for(i=0;i<=Len1-Len2;i++)        {            f1=LCP(i,Len1);            if(f1==Len2)                break;            f1++;            f2=LCP(i+f1,Len1+f1);            if(f1+f2==Len2)                break;            f2++;            f3=LCP(i+f1+f2,Len1+f1+f2);            if(f1+f2+f3==Len2)                break;        }        printf("Case #%d: %d\n",T,((i<=Len1-Len2)?i:-1));    }    return 0;}