CodeForces700E Cool Slogans

感谢dalaoWJZ的讲解。

我们对于每一个串a[i]相当于在他parent的right集合里找一个出现位置在id-len[x]+len[parent]到id[x]-1区间的

用主席树判存在性即可。

至于我问了很久这个区间为什么不能往前判，原因是你每个节点都代表唯一一个子串，所以我们定性查询就可以啦。

By：大奕哥

 #include<bits/stdc++.h>
 using namespace std;
 const int N=4e5+;
 int rt[N],num,n;
 char s[N];
 struct tree{
     int l,r;
 }t[N*];
 void change(int &x,int l,int r,int p)
 {
     if(!x)x=++num;
     if(l==r)return;
     int mid=l+r>>;
     if(p>mid)change(t[x].r,mid+,r,p);
     else change(t[x].l,l,mid,p);
 }
 void merge(int &x,int y)
 {
     if(!x||!y){x=x+y;return;}
     ++num;t[num]=t[x];x=num;
     merge(t[x].l,t[y].l);
     merge(t[x].r,t[y].r);
 }
 bool query(int x,int l,int r,int L,int R)
 {
     if(!x)return ;
     if(l==L&&r==R)return ;
     int mid=l+r>>;
     if(mid<L)return query(t[x].r,mid+,r,L,R);
     else if(mid>R)return query(t[x].l,l,mid,L,R);
     else return query(t[x].l,l,mid,L,mid)|query(t[x].r,mid+,r,mid+,R);
 }
 struct SAM
 {
     int cnt,root,last,c[N][],f[N],dp[N],tt[N],pos[N],r[N],id[N],l[N];
     SAM(){cnt=;last=root=++cnt;}
     void add(int x,int dd)
     {
         int now=last,a=++cnt;last=a;id[a]=dd;
         l[a]=l[now]+;
         for(;now&&!c[now][x];now=f[now])c[now][x]=a;
         if(!now)f[a]=root;
         else{
             int q=c[now][x];
             if(l[q]==l[now]+)f[a]=q;
             else{
                 int b=++cnt;id[b]=dd;
                 l[b]=l[now]+;
                 f[b]=f[q];
                 f[a]=f[q]=b;
                 memcpy(c[b],c[q],sizeof(c[q]));
                 for(;now&&c[now][x]==q;now=f[now])c[now][x]=b;
             }
         }
         return;
     }
     void sort(){
         for(int i=;i<=cnt;++i)tt[l[i]]++;
         for(int i=;i<=n;++i)tt[i]+=tt[i-];
         for(int i=;i<=cnt;++i)r[tt[l[i]]--]=i;
         return;
     }
     void work()
     {
         for(int i=cnt;i>=;--i)
         {
             int x=r[i];
             change(rt[x],,n,id[x]);
             merge(rt[f[x]],rt[x]);
         }int ans=;
         for(int i=;i<=cnt;++i)
         {
             int x=r[i];
             if(f[x]==){pos[x]=x,dp[x]=;continue;}
             if(query(rt[pos[f[x]]],,n,id[x]-l[x]+l[pos[f[x]]],id[x]-))
             dp[x]=dp[pos[f[x]]]+,pos[x]=x;
             else dp[x]=dp[pos[f[x]]],pos[x]=pos[f[x]];
             ans=max(ans,dp[x]);
         }
         printf("%d\n",ans);
         return;
     }
 }A;
 int main()
 {
     scanf("%d%s",&n,s+);
     for(int i=;i<=n;++i)A.add(s[i]-'a',i);
     A.sort();A.work();
     return ;
 }