loj6198谢特 后缀数组+并查集+Trie

先把问题放在后缀数组上考虑

已知两个数组a b，求min(a[i],...,a[j])+(b[i]^b[j])的最大值

套路题

初始每个点都是一个小连通块

把a按从大到小的顺序加入，计算当前加入边作为min的贡献：

每次加入会把两个连通块联通，答案就是两边连通块各出一个数能得到的异或和最大值

我：这不是线性基吗

miaom：mdzz，只能有两个数

我：蛤，好难啊，怎么做啊

miaom：Trie啊

我：哦

没了

 #include <bits/stdc++.h>
 #define N 500001
 #define MAX 18
 using namespace std;
 int NODE,n;char ch;
 int a[N],b[N],c[N*MAX][],size[N],rt[N],lef[N];
 pair<int,int> so[N];
 int X[N],Y[N],SA[N],Rk[N],Ht[N],v[N];
 void GetSA(int a[],int n,int m=)
 {
     int *x=X,*y=Y,i,p,j;
     memset(v,,sizeof v);
     for(i=;i<=n;i++)v[x[i]=a[i]]++;
     for(i=;i<=m;i++)v[i]+=v[i-];
     for(i=n;i>=;i--)SA[v[x[i]]--]=i;
     for(i=;i<=n;i<<=,m=p)
     {
         p=;
         memset(v,,sizeof v);
         for(j=n-i+;j<=n;j++)y[++p]=j;
         for(j=;j<=n;j++)if(SA[j]>i)y[++p]=SA[j]-i;
         for(j=;j<=n;j++)v[x[y[j]]]++;
         for(j=;j<=m;j++)v[j]+=v[j-];
         for(j=n;j>=;j--)SA[v[x[y[j]]]--]=y[j];
         swap(x,y);
         p=;x[SA[]]=;
         for(j=;j<=n;j++)
             if(y[SA[j-]]==y[SA[j]]&&y[SA[j-]+i]==y[SA[j]+i]) x[SA[j]]=p;
                 else x[SA[j]]=++p;
         if(p>=n)break;
     }
 }
 void GetHt(int a[],int n)
 {
     int k=;
     for(int i=;i<=n;i++)Rk[SA[i]]=i;
     for(int i=;i<=n;i++)
     {
         if(k)k--;
         int j=SA[Rk[i]-];
         while(i+k<=n&&j+k<=n&&a[i+k]==a[j+k])k++;
         Ht[Rk[i]]=k;
     }
 }
 int build(int x)
 {
     int rt=++NODE,now=rt;
     for(int i=MAX;i>=;i--)
         c[now][(x>>i)&]=++NODE,now=NODE;
     return rt;
 }
 int que(int x,int y,int z=)
 {
     int best=z;
     if(c[x][])
          if(c[y][]) best=que(c[x][],c[y][],z*+);
          else best=que(c[x][],c[y][],z*);
     if(c[x][])
         if(c[y][]) best=max(best,que(c[x][],c[y][],z*+));
         else best=max(best,que(c[x][],c[y][],z*));
     return best;
 }
 void merge(int x,int y)
 {
     if(c[x][])
          if(c[y][]) merge(c[x][],c[y][]);
          else c[y][]=c[x][];
     if(c[x][])
         if(c[y][]) merge(c[x][],c[y][]);
         else c[y][]=c[x][];
 }
 int getfa(int x)
 {
     if(lef[x]==x) return x;
     else return lef[x]=getfa(lef[x]);
 }
 int main()
 {
     scanf("%d",&n);
     for(ch=getchar();!isalpha(ch);ch=getchar());
     for(int i=;i<=n;i++,ch=getchar())
         a[i]=ch;
     for(int i=;i<=n;i++)
         scanf("%d",&b[i]);
     GetSA(a,n);
     GetHt(a,n);
     for(int i=;i<=n;i++)
         size[i]=,lef[i]=i,rt[i]=build(b[SA[i]]);
     for(int i=;i<=n;i++)
         so[i-]=make_pair(n-Ht[i],i);
     sort(so+,so+n);
     int ret=;
     for(int i=;i<n;i++)
     {
         int x=so[i].second,y=x-,z=n-so[i].first;
         x=getfa(x);y=getfa(y);
         if(size[x]<size[y]) swap(x,y);
         int bes=que(rt[y],rt[x])+z;
         if(bes>ret)
             ret=bes;
         merge(rt[y],rt[x]);
         size[x]+=size[y];
         lef[y]=x;
     }
     printf("%d\n",ret);
     return ;
 }