题链:
http://www.joyoi.cn/problem/tyvj-2301(非权限OI患者,苟且在joyoi。。。)
题解:
后缀自动机,线段树
先对原串建立后缀自动机,不难发现,
会影响答案是那些right集合大小恰好为1的状态。
考虑这些状态是如何影响答案的。
对于一个right集合大小为1的状态s,
令其允许的最大长度为maxs[s],其允许的最小长度为maxs[parent[s]]+1,其right集合里唯一的元素是minr(这里minr表示该状态对应的串在该位置结束)
我们可以得到对应的l=minr-maxs[s]+1,r=minr-maxs[parent[s]],即该状态对应的子串集就是S[l~r,minr]。
显然对于S[r~minr]这些字符,该状态可以给他们贡献一个maxs[parent[s]]+1的答案(对应着只出现过一次的子串S[r,minr])。
而对于S[l~r]这些字符,该状态可以给他们贡献的答案并不是一个相同的值,但是有一个共同点,
就是贡献给第i个字符的答案对应的子串是S[i,minr],即他们有着相同的结尾位置。
考虑到上面的两种贡献都是区间贡献,所以用线段树分别维护两种贡献即可。
(线段树的实现可以考虑先给线段树区间打上永久化标记,最后再通过一遍dfs得出每个叶子节点的最优答案。)
代码:
#include<bits/stdc++.h>
#define MAXN 100005
#define INF 0x3f3f3f3f
using namespace std;
struct SGT{
int rt,size;
int ls[MAXN*2],rs[MAXN*2],lazy[2][MAXN*2];
void Build(int &u,int l,int r){
u=++size;
lazy[0][u]=lazy[1][u]=INF;
if(l==r) return;
int mid=(l+r)/2;
Build(ls[u],l,mid);
Build(rs[u],mid+1,r);
}
void Modify(int u,int l,int r,int al,int ar,int val,int k){
if(al<=l&&r<=ar) return (void)(lazy[k][u]=min(lazy[k][u],val));
int mid=(l+r)/2;
if(al<=mid) Modify(ls[u],l,mid,al,ar,val,k);
if(mid<ar) Modify(rs[u],mid+1,r,al,ar,val,k);
}
void Answer(int u,int l,int r,int minlazy0,int minlazy1){
minlazy0=min(minlazy0,lazy[0][u]);
minlazy1=min(minlazy1,lazy[1][u]);
if(l==r){
printf("%d\n",min(minlazy0,minlazy1-l+1));
return;
}
int mid=(l+r)/2;
Answer(ls[u],l,mid,minlazy0,minlazy1);
Answer(rs[u],mid+1,r,minlazy0,minlazy1);
}
}DT;
struct SAM{
int size;
int maxs[MAXN*3],trans[MAXN*3][26],parent[MAXN*3],minr[MAXN*3],right[MAXN*3];
int Newnode(int a,int b){
++size; maxs[size]=a; minr[size]=INF;
memcpy(trans[size],trans[b],sizeof(trans[b]));
return size;
}
int Extend(int last,int x){
static int p,np,q,nq;
p=last; np=Newnode(maxs[p]+1,0);
for(;p&&!trans[p][x];p=parent[p]) trans[p][x]=np;
if(!p) parent[np]=1;
else{
q=trans[p][x];
if(maxs[p]+1!=maxs[q]){
nq=Newnode(maxs[p]+1,q);
parent[nq]=parent[q];
parent[q]=parent[np]=nq;
for(;p&&trans[p][x]==q;p=parent[p]) trans[p][x]=nq;
}
else parent[np]=q;
}
return np;
}
void Build(char *S){
static int p=1,last,len,tmp[MAXN],order[MAXN*3];
memset(trans[0],0,sizeof(trans[0]));
size=0; last=Newnode(0,0); len=strlen(S);
for(int i=0;i<len;i++) last=Extend(last,S[i]-'a');
for(int i=0;i<len;i++) p=trans[p][S[i]-'a'],minr[p]=i,right[p]=1;
for(int i=1;i<=size;i++) tmp[maxs[i]]++;
for(int i=1;i<=len;i++) tmp[i]+=tmp[i-1];
for(int i=1;i<=size;i++) order[tmp[maxs[i]]--]=i;
for(int i=size;i;i--){
p=order[i];
minr[parent[p]]=min(minr[parent[p]],minr[p]);
right[parent[p]]+=right[p];
}
}
void Solve(int len){
for(int i=1,l,r;i<=size;i++) if(right[i]==1){
l=minr[i]-maxs[i]+1; r=minr[i]-maxs[parent[i]];
DT.Modify(DT.rt,0,len-1,r,minr[i],maxs[parent[i]]+1,0);
DT.Modify(DT.rt,0,len-1,l,r,minr[i],1);
}
}
}SUF;
int main(){
static char S[MAXN];
scanf("%s",S);
int len=strlen(S);
DT.Build(DT.rt,0,len-1);
SUF.Build(S);
SUF.Solve(len);
DT.Answer(DT.rt,0,len-1,INF,INF);
return 0;
}