CF700E Cool Slogans

CF700E Cool Slogans

• 给定一个字符串 $S$，要求构造字符串序列 $s_1,s_2,\dots ,s_k$满足任意$s_i$ 都是$S$的子串，且$s_i$在$s_{i-1}$中出现至少两次。
• 求最大$k$
• $|S|\le 2e5$

Solution

• 首先所有的$s_i$的序列都可以缩成$s_i$是$s_{i-1}$后缀的形式，所以考虑直接在后缀树上进行DP。
• 需要证明一个重要的结论：对于SAM上的节点$x$，只有它的$maxlen$是可能作为转移状态的。
  • 只需要证明对于$x$的祖先节点$y$，任意$y$的$[minlen,maxlen]$中的字符串在$x$代表的最长字符串中出现次数都相同。
  • 考虑如果$a,b$都是$y$代表的字符串，$|b|=|a|+1$，$b$比$a$少出现一次，多出来的字符一定在$st{r}_{max}x$的外面，那么由于已经是最长的了，所以这个长度为$|st{r}_{max}x|+1$的串属于$z$，$z$为$x$子树内节点，那么存在一个$endpos=endpos(z)-endpos(x)$使得存在$b$而不存在$a$，与$a,b$相同$endpos$矛盾。
• 那么接下来就直接可持久化线段树合并$endpos$，再根据单调性将转移到当前节点$x$的$y$贪心地移动即可。判断只需要在随便选择一个$endpos(x)$，判断这个区间内的情况就好了。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define maxn 400005
#define maxm 10000005
using namespace std;

int n,i,j,k;
struct node{
	int to[26],len,end,fa;
} T[maxn];
int tot,las,Rt;
void add(int c,int end){
	int np=++tot,p=las; las=np;
	T[np].len=T[p].len+1,T[np].end=end;
	for(;p&&!T[p].to[c];p=T[p].fa) T[p].to[c]=np;
	if (!p) T[np].fa=Rt; else {
		int q=T[p].to[c];
		if (T[q].len==T[p].len+1) T[np].fa=q; else {
			int nq=++tot; T[nq]=T[q],T[nq].end=0;
			T[nq].len=T[p].len+1,T[q].fa=T[np].fa=nq;
			for(;p&&T[p].to[c]==q;p=T[p].fa) T[p].to[c]=nq;
		}
	}
}

int em,e[maxn],nx[maxn],ls[maxn];
void insert(int x,int y){
	em++; e[em]=y; nx[em]=ls[x]; ls[x]=em;
}

int t[maxm],tl[maxm],tr[maxm],cnt,rt[maxn];
void change(int &x,int l,int r,int p){
	if (!x) x=++cnt; t[x]++;
	if (l==r) return;
	int mid=(l+r)>>1;
	if (p<=mid) change(tl[x],l,mid,p);
	else change(tr[x],mid+1,r,p);
}

int find(int x,int l,int r,int L,int R){
	if (!x||l>R||r<L) return 0;
	if (L<=l&&r<=R) return t[x];
	int mid=(l+r)>>1;
	return find(tl[x],l,mid,L,R)+find(tr[x],mid+1,r,L,R);
}

void merge(int &x,int y,int l,int r){
	t[++cnt]=t[x]+t[y],tl[cnt]=tl[x],tr[cnt]=tr[x],x=cnt;
	if (l==r) return;
	int mid=(l+r)>>1;
	if (tl[x]&&tl[y]) merge(tl[x],tl[y],l,mid);
	else if (tl[y]) tl[x]=tl[y];
	if (tr[x]&&tr[y]) merge(tr[x],tr[y],mid+1,r);
	else if (tr[y]) tr[x]=tr[y];
}

void cover(int x){
	if (tl[x]) cover(tl[x]);
	if (tr[x]) cover(tr[x]);
}

void dfs(int x){
	rt[x]=++cnt;
	if (T[x].end) change(rt[x],1,n,T[x].end);
	for(int i=ls[x];i;i=nx[i]) {
		dfs(e[i]);
		merge(rt[x],rt[e[i]],1,n);
		T[x].end=T[e[i]].end;
	}
}

int d[maxn],f[maxn],ans;
int check(int y,int x){
	return find(rt[y],1,n,T[x].end-T[x].len+T[y].len,T[x].end-1)>0;
}

void dp(int x,int y){
	while (y<d[0]&&check(d[y+1],x)) y++;
	f[x]=f[d[y]]+1;
	ans=max(ans,f[x]);
	d[++d[0]]=x;
	for(int i=ls[x];i;i=nx[i]) dp(e[i],y);
	d[0]--;
}

int main(){
	freopen("ceshi.in","r",stdin);
	scanf("%d",&n);
	char ch=getchar(); while (ch<'a'||ch>'z') ch=getchar();
	tot=Rt=las=1;
	for(i=1;i<=n;i++) add(ch-'a',i),ch=getchar();
	for(i=1;i<=tot;i++) if (T[i].fa) 
		insert(T[i].fa,i);
	dfs(Rt);
	for(i=ls[Rt];i;i=nx[i]) {
		dp(e[i],0);
	}
	printf("%d",ans);
}