[Codeforces]862F - Mahmoud and Ehab and the final stage

题目大意：n个字符串，支持修改一个位置上的字符串和查询一个区间的子区间中长度乘LCP的最大值，输入字符数和询问数不超过10^5。

做法：求出相邻的LCP长度，区间LCP等于区间最小值，查询分几种情况考虑，一种只有一个串，线段树维护长度最大值即可；若有若干个串，设一个阈值k，若答案的LCP<=k，对于小等k的每一个i，若一个位置的相邻LCP大等i，设为1，否则设为0，即求区间最长连续1，每种i开一棵线段树维护即可；若LCP>k，我们把相邻LCP长度超过k的位置存进set，查询的时候拿出来，直接建笛卡尔树计算答案，由于字符总数有限，这样的位置不会超过L/k个。总时间复杂度为O(nklogn+nL/k)，适当调整k，时间复杂度为O(n(nlogn)^0.5)。

代码：

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<queue>
#include<map>
#include<set>
using namespace std;
inline int read()
{
    int x=,f=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
#define MN 100000
#define N 131072
#define K 50
struct node{int l,r,u,mx;}T[K+][N*+];
node operator+(const node&a,const node&b)
{
    return (node){a.l+a.u*b.l,b.r+b.u*a.r,a.u*b.u,max(max(a.mx,b.mx),a.r+b.l)};
}
string s[MN+];
int t[N*+],v[MN+],a[MN+],an,Z[MN+],zn,L[MN+],R[MN+],ans,S[MN+];
set<int> st;
void change(int k,int x){for(t[k+=N]=x;k>>=;)t[k]=max(t[k<<],t[k<<|]);}
int query(int l,int r)
{
    int res=;
    for(l+=N-,r+=N+;l^r^;l>>=,r>>=)
    {
        if(~l&)res=max(res,t[l+]);
        if( r&)res=max(res,t[r-]);
    }
    return res;
}
void change(node*t,int k,node x){for(t[k+=N]=x;k>>=;)t[k]=t[k<<]+t[k<<|];}
node Query(node*t,int l,int r)
{
    node L,R;int ul=,ur=;
    for(l+=N-,r+=N+;l^r^;l>>=,r>>=)
    {
        if(~l&)L=ul?L+t[l+]:(ul=,t[l+]);
        if( r&)R=ur?t[r-]+R:(ur=,t[r-]);
    }
    return ul?ur?L+R:L:R;
}
void renew(int x)
{
    cin>>s[x];change(x,s[x].size());
    if(v[x]>K)st.erase(x);if(v[x+]>K)st.erase(x+);
    for(v[x]=;v[x]<s[x].size()&&v[x]<s[x-].size()&&s[x][v[x]]==s[x-][v[x]];)++v[x];
    for(v[x+]=;v[x+]<s[x].size()&&v[x+]<s[x+].size()&&s[x][v[x+]]==s[x+][v[x+]];)++v[x+];
    if(v[x]>K)st.insert(x);if(v[x+]>K)st.insert(x+);
    node a=(node){,,,},b=(node){,,,};
    for(int i=;i<=K;++i)change(T[i],x,v[x]<i?a:b),change(T[i],x+,v[x+]<i?a:b);
}
void dfs(int x)
{
    S[x]=;
    if(L[x])dfs(L[x]),S[x]+=S[L[x]];
    if(R[x])dfs(R[x]),S[x]+=S[R[x]];
    ans=max(ans,v[x]*(S[x]+));
}
int p(const set<int>::iterator&i){return i!=st.end()?*i:MN+;}
int main()
{
    int n,m,i,l,r,z;
    n=read();m=read();
    for(i=;i<=n;++i)renew(i);
    while(m--)
        if(read()==)
        {
            l=read();r=read();
            ans=query(l,r);
            if(l<r)for(i=;i<=K;++i)z=Query(T[i],l+,r).mx,ans=max(ans,i*(z?z+:));
            set<int>::iterator it=st.upper_bound(l);
            for(z=p(it);z<=r;)
            {
                for(a[an=]=z;(z=p(++it))<=r&&z==a[an]+;)a[++an]=z;
                for(i=,zn=;i<=an;++i)
                {
                    L[a[i]]=R[a[i]]=;
                    while(zn&&v[a[i]]<v[a[Z[zn]]])L[a[i]]=a[Z[zn--]];
                    R[a[Z[zn]]]=a[i];Z[++zn]=i;
                }
                dfs(a[Z[]]);
            }
            printf("%d\n",ans);
        }
        else renew(read());
}