正题
题目链接:https://www.luogu.com.cn/problem/P2178
题目大意
长度为nnn的字符串,每个位置有一个权值。
如果对于q,pq,pq,p有sq,q+r−1=sp,p+r−1s_{q,q+r-1}=s_{p,p+r-1}sq,q+r−1=sp,p+r−1那么就称q,pq,pq,p为rrr相似。
对于每个rrr求有多少对q,pq,pq,p是rrr相似并且求一对rrr相似的p,qp,qp,q似的权值的乘积最大。
解题思路
如果两个后缀的LCP>=rLCP>=rLCP>=r那么这两个后缀就是rrr相似的。
然后考虑排序后的heightheightheight,如果我们按照heightheightheight从大到小合并那么就说明这次合并时可以使得左右两个集合heightiheight_iheighti相似。
开一个并查集然后按照heightheightheight从大到小排序然后一个一个合并即可
codecodecode
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=3e5+10;
ll n,m,p[N],siz[N],val[N],wal[N];
ll c[N],x[N],sa[N],y[N],rk[N],h[N];
ll ans[N],ans1[N],ans2[N],fa[N];
char s[N];
void Qsort(){for(ll i=1;i<=m;i++)c[i]=0;for(ll i=1;i<=n;i++)c[x[i]]++;for(ll i=1;i<=m;i++)c[i]+=c[i-1];for(ll i=n;i>=1;i--)sa[c[x[y[i]]]--]=y[i],y[i]=0;return;
}
void Get_SA(){m=26;for(ll i=1;i<=n;i++)x[i]=s[i]-'a'+1,y[i]=i;Qsort();for(ll w=1;w<=n;w<<=1){ll p=0;for(ll i=n-w+1;i<=n;i++)y[++p]=i;for(ll i=1;i<=n;i++)if(sa[i]>w)y[++p]=sa[i]-w;Qsort();swap(x,y);x[sa[1]]=p=1;for(ll i=2;i<=n;i++)x[sa[i]]=(y[sa[i]]==y[sa[i-1]]&&y[sa[i]+w]==y[sa[i-1]+w])?p:++p;if(p==n)break;m=p;}return;
}
void Get_Height(){ll k=0;for(ll i=1;i<=n;i++)rk[sa[i]]=i;for(ll i=1;i<=n;i++){if(rk[i]==1)continue;if(k)k--;ll j=sa[rk[i]-1];while(i+k<=n&&j+k<=n&&s[i+k]==s[j+k])k++;h[rk[i]]=k;}return;
}
bool cmp(ll x,ll y)
{return h[x]>h[y];}
ll find(ll x)
{return fa[x]==x?(x):(fa[x]=find(fa[x]));}
void unionn(ll x,ll y,ll h){ans1[h]+=siz[x]*siz[y];ans[x]=max(ans[x],ans[y]);ans[x]=max(ans[x],max(val[x]*val[y],wal[x]*wal[y]));ans[x]=max(ans[x],max(val[x]*wal[y],wal[x]*val[y]));val[x]=min(val[x],val[y]);wal[x]=max(wal[x],wal[y]);fa[y]=x;siz[x]+=siz[y];ans2[h]=max(ans2[h],ans[x]);return;
}
int main()
{scanf("%lld",&n);scanf("%s",s+1);for(ll i=1;i<=n;i++)scanf("%lld",&val[i]),p[i]=i,siz[i]=1,wal[i]=val[i],fa[i]=i;Get_SA();Get_Height();sort(p+2,p+1+n,cmp);memset(ans2,0xcf,sizeof(ans2));memset(ans,0xcf,sizeof(ans));for(ll i=2;i<=n;i++){ll x=find(sa[p[i]-1]),y=find(sa[p[i]]);if(x==y)continue;unionn(x,y,h[p[i]]);}for(ll i=n-1;i>=0;i--)ans1[i]+=ans1[i+1],ans2[i]=max(ans2[i],ans2[i+1]);for(ll i=0;i<n;i++)printf("%lld %lld\n",ans1[i],(!ans1[i])?0:ans2[i]);
}