【poj 3461】Oulipo hash求子串hash值

最新推荐文章于 2023-06-25 08:42:03 发布

pbihao

最新推荐文章于 2023-06-25 08:42:03 发布

阅读量271

点赞数

分类专栏： poj 字符串

本文链接：https://blog.csdn.net/pbihao/article/details/53122418

版权

字符串同时被 2 个专栏收录

11 篇文章 0 订阅

订阅专栏

poj

10 篇文章 0 订阅

订阅专栏

我靠，这道题卡时限太厉害了吧，我用大质数取模居然能T，最好无赖改成了unsigned long long自然溢出才过，不过这样写hash很容易被卡啊。用hash过一次

字串hash公式 hash[l,r]=(hash[r]-hash[l-1]*p^(r-l+1)+Mod)%Mod

#include<cstdio>
#include<cstring>
#include<iostream>
#define LL unsigned long long
using namespace std;
int T;
const int p=31,Mod=1000050;
LL pp[1000020],h[1000010];
char s[1000020],c[1000020];

int main(){
	scanf("%d",&T);
	pp[1]=p;for(int i=2;i<=1000000;i++)pp[i]=pp[i-1]*p;
	while(T--){
		LL hashs=0;
		int ans=0;
		scanf("%s%s",s,c);
		LL lens=strlen(s),lenc=strlen(c);
		for(int i=0;i<lens;i++)hashs=hashs*p+s[i];
		for(int i=0;i<lenc;i++){
			if(i==0)h[i]=c[i];
			else h[i]=(h[i-1]*p+c[i]);
		}
		
		
		for(int i=0;i<=lenc-lens;i++){
			if(hashs==(h[i+lens-1]-(i==0 ? 0 : h[i-1])*pp[lens]))
				ans++;
		}
		printf("%d\n",ans);
	}
	return 0;
}

KMP:

#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
char c[10020],s[1000020];
int sl,cl,nex[20030],m;

void get_nex(){
	nex[1]=nex[0]=0;
	for(int i=1;i<cl;i++){
		int j=nex[i];
		while(j&&c[j]!=c[i])j=nex[j];
		nex[i+1]=c[i]==c[j]?j+1:0;
	}
}
int solve(){
	int ans=0,j=0;
	for(int i=0;i<sl;i++){
		while(j&&c[j]!=s[i])j=nex[j];
		if(c[j]==s[i])j++;
		if(j==cl){
			ans++;
			j=nex[j];
		}
	}
	return ans;
}

int main(){
	scanf("%d",&m);
	while(m--){
		scanf("%s%s",c,s);
		sl=strlen(s),cl=strlen(c);
		get_nex();
		printf("%d\n",solve());
	}
	return 0;
}