BNUOJ 34490 Justice String
这是今年北京邀请赛的一个题目。
问串1是否有与串2“匹配”的子串。这里“匹配”是可以有最多两个字符不同的匹配。
“用两秒钟构思”,这不是后缀数组水题么。
然后我就水啊水,用的DC3,时限2秒,n=200000,nlogn的算法。
尼玛,什么破OJ啊,这样都被卡。
好吧。穷则思变。
记得后缀自动机陈立杰的PPT里说的那个二分哈希求解LCS的方法么。
O(n)预处理,O(1)可以得到所有哈希值。哈希值相同的两串视为相同(实际上有很小几率犯错,忽略)
二分得到最大匹配长度,判定长度是否大于串2长,跳1,重复三次。搞定。
同样是nlogn,常数小的算法还是很有学习的必要的。
这是今年北京邀请赛的一个题目。
问串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;
}