LA 4513 Stammering Aliens

最新推荐文章于 2021-08-04 15:15:41 发布

squee_spoon

最新推荐文章于 2021-08-04 15:15:41 发布

阅读量510

点赞数

CC 4.0 BY-SA版权

分类专栏：数据结构

本文链接：https://blog.youkuaiyun.com/squee_spoon/article/details/40350099

数据结构专栏收录该内容

73 篇文章

订阅专栏

本文介绍了一种基于哈希值的LCP算法，用于找出字符串中至少出现m次的最长子串。通过使用后缀哈希，可以实现子串比较的O(1)时间复杂度。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：给一个串，找出至少出现m次的最长子串。

思路：基于哈希值的LCP算法。从白书上学来的，真是一个神奇的技能。。利用后缀哈希H[i]=H[i+1]x+s[i]，为每个长度为l的字符串计算一个哈希值，就可以在O(1)的时间内比较两个子串。

#include <iostream>       
#include <stdio.h>       
#include <cmath>       
#include <algorithm>       
#include <iomanip>       
#include <cstdlib>       
#include <string>       
#include <string.h>       
#include <vector>       
#include <queue>       
#include <stack>       
#include <map>     
#include <set>     
#include <ctype.h>       
#include <sstream>   
#define INF 1000000000   
#define ll long long   
#define min3(a,b,c) min(a,min(b,c))
#define max3(a,b,c) max(a,max(b,c))
     
using namespace std;

int m;
char str[40010];
int siz;

unsigned long long hash[40010];
unsigned long long H[40010],xn[40010];

int rank[40010];

bool cmp(const int& a,const int& b){
	return hash[a]<hash[b]||(hash[a]==hash[b]&&a<b);
}

int judge(int ans){
	for(int i=0;i<=siz-ans;i++){
		hash[i]=H[i]-H[i+ans]*xn[ans];
		rank[i]=i;
	}
	sort(rank,rank+siz-ans+1,cmp);
	int c=0;
	int pos=-1;
	for(int i=0;i<=siz-ans;i++){
		if(i==0||hash[rank[i]]!=hash[rank[i-1]])c=0;
		if(++c>=m)pos=max(pos,rank[i]);
	}
	return pos;
}

int main(){
	const int x=123;
	xn[0]=1;
	for(int i=1;i<=40000;i++)xn[i]=xn[i-1]*x;
	
	
	while(cin>>m){
		if(!m)break;
		memset(H,0,sizeof(H));
		
		scanf("%s",str);
		siz=strlen(str);
		
		for(int i=siz-1;i>=0;i--){
			H[i]=H[i+1]*x+str[i]-'a';
		}
		
		int l=1;
		int r=siz-m+1;
		
		int ans,pos=-1;
		while(l<=r){
			int mid=(l+r)/2;
			int tmp=judge(mid);
			if(tmp!=-1){
				l=mid+1;
				ans=mid;
				pos=tmp;
			}else{
				r=mid-1;
			}
		}
		if(pos!=-1){
			cout<<ans<<" "<<pos<<endl;
		}else{
			cout<<"none"<<endl;
		}
	}
	return 0;
}