【BZOJ1396】识别子串&【BZOJ2865】字符串识别（后缀自动机）

【BZOJ1396】识别子串&【BZOJ2865】字符串识别（后缀自动机）

题面

自从有了DBZOJ
终于有地方交权限题了

题解

很明显，只出现了一次的串
在\(SAM\)的\(right/endpos\)集合大小一定为\(1\)
换句话说，在\(parent\)树上是叶子节点
找到所有这样的节点，
假设它的\(len=r\)，它父亲的\(len=p\)，它的结束位置为显然就是\(r\)
令\(l=r-p\)
以\(r\)结尾，
并且只出现了一次的串的左端点
为\(1..l\)，那么，他们的答案可以更新为\(r+1-i\)
剩下的位置\(l+1..r\)，他们无法作为左端点，只能包含在这些串中
于是找到一个最短的包含他们的串\(S[l..r]\)
所以，这段区间的答案可以更新为\(r-l+1\)

显然不好一起维护，
于是开两棵线段树，一个维护\(r+1-i\),先不考虑\(i\)，最后减去就好
另一个直接维护\(r-l+1\)

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define lson (now<<1)
#define rson (now<<1|1)
#define MAX 111111
inline int read()
{
    RG int x=0,t=1;RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
int n;
bool vis[MAX<<1];
struct SAM
{
	struct Node
	{
		int son[26];
		int ff,len;
	}t[MAX<<1];
	int last,tot;
	void init(){last=tot=1;}
	void extend(int c)
	{
		int p=last,np=++tot;last=np;
		t[np].len=t[p].len+1;
		while(p&&!t[p].son[c])t[p].son[c]=np,p=t[p].ff;
		if(!p)t[np].ff=1;
		else
		{
			int q=t[p].son[c];
			if(t[q].len==t[p].len+1)t[np].ff=q;
			else
			{
				int nq=++tot;
				t[nq]=t[q];t[nq].len=t[p].len+1;
				t[np].ff=t[q].ff=nq;
				while(p&&t[p].son[c]==q)t[p].son[c]=nq,p=t[p].ff;
			}
		}
	}
}SAM;
char ch[MAX];
struct SegMentTree
{
	struct Node{int v;}t[MAX<<2];
	void Build(int now,int l,int r)
	{
		t[now].v=1e9;if(l==r)return;
		int mid=(l+r)>>1;
		Build(lson,l,mid);Build(rson,mid+1,r);
	}
	void puttag(int now,int w){t[now].v=min(t[now].v,w);}
	void pushdown(int now)
	{
		if(t[now].v==1e9)return;
		puttag(lson,t[now].v);puttag(rson,t[now].v);
		t[now].v=1e9;
	}
	void Modify(int now,int l,int r,int L,int R,int w)
	{
		if(L>R)return;
		if(L<=l&&r<=R){puttag(now,w);return;}
		pushdown(now);int mid=(l+r)>>1;
		if(L<=mid)Modify(lson,l,mid,L,R,w);
		if(R>mid)Modify(rson,mid+1,r,L,R,w);
	}
	int Query(int now,int l,int r,int p)
	{
		if(l==r)return t[now].v;
		pushdown(now);int mid=(l+r)>>1;
		if(p<=mid)return Query(lson,l,mid,p);
		else return Query(rson,mid+1,r,p);
	}
}A,B;
int main()
{
	SAM.init();
	scanf("%s",ch+1);n=strlen(ch+1);
	for(int i=1;i<=n;++i)SAM.extend(ch[i]-97);
	for(int i=1;i<=SAM.tot;++i)vis[SAM.t[i].ff]=true;
	A.Build(1,1,n);B.Build(1,1,n);
	for(int i=1;i<=SAM.tot;++i)
		if(!vis[i])
		{
			int l=SAM.t[i].len-SAM.t[SAM.t[i].ff].len,r=SAM.t[i].len;
			A.Modify(1,1,n,l,r,r-l+1);
			B.Modify(1,1,n,1,l-1,r+1);
		}
	for(int i=1;i<=n;++i)
		printf("%d\n",min(A.Query(1,1,n,i),B.Query(1,1,n,i)-i));
	return 0;
}