BNUOJ 34990 Justice String (基于hash的LCP)

最新推荐文章于 2024-09-20 21:00:12 发布

转载最新推荐文章于 2024-09-20 21:00:12 发布 · 226 阅读

---hash 专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一种字符串匹配算法，通过枚举起点求子串与目标字符串的最长公共前缀（LCP），并使用哈希技巧来预处理字符串以降低计算复杂度。通过二分或枚举方法找到满足条件的子串起始位置。

转自http://blog.youkuaiyun.com/w20810/article/details/48438047

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

分析：枚举起点求s和B的LCP。看s最多修改两次后是否满足条件。

hash的用法挺多的，原理都是预处理降低复杂度。
然后就可以二分或者枚举了。。

    #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;  
    }