BNUOJ 34490 Justice String，字符串哈希，二分

BNUOJ 34490 Justice String

这是今年北京邀请赛的一个题目。
问串1是否有与串2“匹配”的子串。这里“匹配”是可以有最多两个字符不同的匹配。

“用两秒钟构思”，这不是后缀数组水题么。
然后我就水啊水，用的DC3,时限2秒，n=200000，nlogn的算法。
尼玛，什么破OJ啊，这样都被卡。

好吧。穷则思变。
记得后缀自动机陈立杰的PPT里说的那个二分哈希求解LCS的方法么。
O(n)预处理，O(1)可以得到所有哈希值。哈希值相同的两串视为相同（实际上有很小几率犯错，忽略）
二分得到最大匹配长度，判定长度是否大于串2长，跳1，重复三次。搞定。
同样是nlogn，常数小的算法还是很有学习的必要的。

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
#define Maxn 101000

unsigned int base[Maxn];
unsigned int h1[Maxn],h2[Maxn];
const unsigned int MAGIC=39;

char s1[Maxn],s2[Maxn];

inline void init_hash(int l,char *s,unsigned int *h){//s 0..l-1 -> h 1..l
    int i;
    h[0]=0;
    for(i=1;i<=l;++i) h[i]=h[i-1]*MAGIC+(s[i-1]-'a');
    base[0]=1;
    for(i=1;i<=l;++i) base[i]=base[i-1]*MAGIC;
}

inline unsigned int string_hash(unsigned int *h,int l,int r){//(l,r]
    return h[r]-h[l]*base[r-l];
}

inline bool comp(unsigned int *h1,int l1,int r1,unsigned int *h2,int l2,int r2){
    unsigned int u1=string_hash(h1,l1,r1+1);
    unsigned int u2=string_hash(h2,l2,r2+1);
    if (u1==u2) return true;
    else return false;
}

int findmax(unsigned int *h1,int l1,unsigned int *h2,int l2,int len){
    int l=1,r=len-l2,m,ans=0;
    while(l<=r){
        m=(l+r)/2;
        if (comp(h1,l1,l1+m-1,h2,l2,l2+m-1)){
            if (ans<m) ans=m;
            l=m+1;
        }
        else r=m-1;
    }
    return ans;
}



int main(){
    //freopen("jin.txt","r",stdin);
    int tcas,cas,l1,l2,i,tim,flag,p1,p2,ans,tmp1;

    scanf("%d",&tcas);
    for(cas=1;cas<=tcas;++cas){
        scanf("%s%s",s1,s2);
        l1=strlen(s1);
        l2=strlen(s2);
        if (l1<l2) {
            printf("Case #%d: %d\n",cas,-1);
            continue;
        }
        init_hash(l1,s1,h1);
        init_hash(l2,s2,h2);

        ans=-1;
        for(i=0;i<l1;++i){
            if (i+l2>l1) break;
            tim=0;
            flag=0;
            p1=0;
            p2=0;
            while(tim<=2){
                if (p2>=l2) {flag=1;break;}
                tmp1=findmax(h1,p1+i,h2,p2,l2);
                p1=p1+tmp1;
                p2=p2+tmp1;
                if (p2>=l2) {flag=1;break;}
                p1+=1;
                p2+=1;
                tim++;
            }
            if (flag==1) {ans=i;break;}
        }
        printf("Case #%d: %d\n",cas,ans);
    }
    return 0;
}