其他分享
首页 > 其他分享> > SAM代补题

SAM代补题

作者:互联网

Hacker

对模式串建立 SAM ,将匹配串的字符一个个走下去,没有该字符就向上跳 parent tree 上的父亲继续找,如此得到对于每个前缀 b1,i 的可最长匹配的后缀,加个线段树维护权值前缀和的最小值即可。

#include<bits/stdc++.h>
#define IL inline
#define LL long long
using namespace std;
const int N=2e5+3;
const LL inf=1e18;
struct node{
	int ch[26],fa,len;
};
int n,m,k,ml,pos;
LL val[N],ans;
char s[N];
struct segment{
	LL Min[N<<2];
	#define ls k<<1
	#define rs k<<1|1
	IL void pushup(int k){Min[k]=min(Min[ls],Min[rs]);}
	void build(int k,int l,int r){
		if(l==r){Min[k]=val[l];return;}
		int mid=l+r>>1;
		build(ls,l,mid),build(rs,mid+1,r);
		pushup(k);
	}
	LL query(int k,int l,int r,int ll,int rr){
		if(l>=ll&&r<=rr) return Min[k];
		int mid=l+r>>1;
		if(ll>mid) return query(rs,mid+1,r,ll,rr);
		if(rr<=mid) return query(ls,l,mid,ll,rr);
		return min(query(ls,l,mid,ll,rr),query(rs,mid+1,r,ll,rr));
	}
}T;
struct SAM{
	node a[N];int las=1,cnt=1;
	IL void ins(int pos,int c){
		int p=las,np;a[np=las=++cnt].len=pos;
		while(p&&!a[p].ch[c]) a[p].ch[c]=np,p=a[p].fa;
		if(!p){a[np].fa=1;return;}
		int q=a[p].ch[c];
		if(a[q].len==a[p].len+1){a[np].fa=q;return;}
		int nq=++cnt;a[nq]=a[q],a[nq].len=a[p].len+1;
		a[np].fa=a[q].fa=nq;
		while(p&&a[p].ch[c]==q) a[p].ch[c]=nq,p=a[p].fa;
	}
	void go(int wei,int c){
		while(pos^1&&!a[pos].ch[c]) pos=a[pos].fa,ml=a[pos].len;
		if(a[pos].ch[c]) pos=a[pos].ch[c],++ml;
		ans=max(ans,val[wei]-T.query(1,0,m,wei-ml,wei));
	}
}S;
IL int in(){
	char c;int f=1;
	while((c=getchar())<'0'||c>'9')
	  if(c=='-') f=-1;
	int x=c-'0';
	while((c=getchar())>='0'&&c<='9')
	  x=x*10+c-'0';
	return x*f;
}
int main()
{
	n=in(),m=in(),k=in();
	scanf("%s",s+1);
	for(int i=1;i<=n;++i) S.ins(i,s[i]-'a');
	for(int i=1;i<=m;++i) val[i]=val[i-1]+in();
	T.build(1,0,m);
	while(k--){
		scanf("%s",s+1),pos=1,ml=0,ans=0;
		for(int i=1;i<=m;++i) S.go(i,s[i]-'a');
		printf("%lld\n",ans);
	}
	return 0;
}

标签:const,rr,SAM,int,LL,mid,补题,ll
来源: https://www.cnblogs.com/liang302/p/16644847.html