#1449 : 后缀自动机三·重复旋律6 后缀自动机 某个长度的字串,在原串中出现的次数的最大值

先用SAM建出的后缀树的性质,求出每个endpos集合的大小。

然后用线段树或线性更新求出每个长度出现次数的最大值即可。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int M=1e6+7;

int head[M*2],cnt=1;//点集个数为2n 
void init(int n){cnt=1;for(int i=0;i<=n;i++)head[i]=0;}
struct EDGE{int to,nxt,w;}ee[M*2*2];
void add(int x,int y,int w){ee[++cnt].nxt=head[x],ee[cnt].w=w,ee[cnt].to=y,head[x]=cnt;}

struct SAM {
	int maxlen[M*2], trans[M*2][26], link[M*2], tg[M*2], size, last;
	SAM() {
	 	size = last = 1; 
		memset(tg,0,sizeof(tg));//注意这里要初始化tg 
	}
	inline void extend(int id) {
		int cur = (++ size), p;
		tg[cur]=1;//标记当前点,表示当前集合包含原串的前缀 
		maxlen[cur] = maxlen[last] + 1;
		for (p = last; p && !trans[p][id]; p = link[p]) trans[p][id] = cur;
		if (!p) link[cur] = 1;
		else {
			int q = trans[p][id];
			if (maxlen[q] == maxlen[p] + 1) link[cur] = q;
			else {
				int clone = (++ size);
				maxlen[clone] = maxlen[p] + 1;
				for(int i=0;i<26;i++)trans[clone][i]=trans[q][i];
				link[clone] = link[q];
				for (; p && trans[p][id] == q; p = link[p]) trans[p][id] = clone;
				link[cur] = link[q] = clone;
			}
		} 
		last = cur;
	}
	/*inline void gao(){//本质不同字符字串个数 
        ll ans=0;
        for(int i=2;i<=last;++i)ans+=maxlen[i]-maxlen[link[i]];
        printf("%lld\n",ans);
    }*/ 
    int nm[M*2],mx[M];
    void dfs(int x,int fa){
    	nm[x]=tg[x];
    	for(int i=head[x];i;i=ee[i].nxt){
    		int y=ee[i].to;
    		if(y==fa)continue;
    		dfs(y,x);
    		nm[x]+=nm[y];
		}
	}
	void gao(int n){
		for(int i=2;i<=size;i++)add(i,link[i],1),add(link[i],i,1);
		dfs(1,0);
		for(int i=2;i<=size;i++)
			mx[maxlen[i]]=max(mx[maxlen[i]],nm[i]);
			//cout<<i<<" "<<maxlen[i]<<" "<<nm[i]<<endl;;
		for(int i=n-1;i>=1;i--)mx[i]=max(mx[i],mx[i+1]);
		for(int i=1;i<=n;i++)printf("%d\n",mx[i]);
	} 
}sam;
char s[M];
int main(){
	scanf("%s",s+1);
	int n=strlen(s+1);
	for(int i=1;i<=n;i++)sam.extend(s[i]-'a');
    sam.gao(n);
    return 0;
}

 

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值