2017.3.6[hihocoder#1415]后缀数组三·重复旋律3

http://hihocoder.com/problemset/problem/1415

#include<cmath>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 200010
#define RG register
#define inf 0x3f3f3f3f
#define Inf 99999999999999999LL
using namespace std;
typedef long long LL;
char c[2*N],c1[N],c2[N];
int ans,len,len1,len2,a[N],b[N],sa[N],ssa[N],cnta[N],cntb[N],heit[N],order[N];
inline int Abs(RG const int &a){return a>0?a:-a;}
inline int Max(RG const int &a,RG const int &b){return a>b?a:b;}
inline int Min(RG const int &a,RG const int &b){return a>b?b:a;}
inline int gi(){
    RG int x=0;RG bool flag=0;RG char c=getchar();
    while((c<'0'||c>'9')&&c!='-') c=getchar();
    if(c=='-') c=getchar(),flag=1;
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return flag?-x:x;
}
inline void getsa(){
    for (RG int i=1;i<=len;++i) ++cnta[c[i]];
    for (RG int i=1;i<=200;++i) cnta[i]+=cnta[i-1];
    for (RG int i=len;i;--i)    sa[cnta[c[i]]--]=i;
    order[sa[1]]=1;
    if(len>1){
    for (RG int i=2;i<=len;++i){
        order[sa[i]]=order[sa[i-1]];
        if(c[sa[i]]!=c[sa[i-1]]) ++order[sa[i]];
    }
    for (RG int now=1;order[sa[len]]<len;now<<=1){
        memset(cnta,0,sizeof(cnta));
        memset(cntb,0,sizeof(cntb));
        for (RG int i=1;i<=len;++i){
        ++cnta[a[i]=order[i]];
        ++cntb[b[i]=(i+now<=len)?order[i+now]:0];
        }
        for (RG int i=1;i<=len;++i) cntb[i]+=cntb[i-1];
        for (RG int i=len;i;--i)    ssa[cntb[b[i]]--]=i;
        for (RG int i=1;i<=len;++i) cnta[i]+=cnta[i-1];
        for (RG int i=len;i;--i)    sa[cnta[a[ssa[i]]]--]=ssa[i];
        order[sa[1]]=1;
        for (RG int i=2;i<=len;++i){
        order[sa[i]]=order[sa[i-1]];
        if(a[sa[i]]!=a[sa[i-1]]||b[sa[i]]!=b[sa[i-1]]) ++order[sa[i]];
        }
    }
    }
    for (RG int i=1,j=0;i<=len;++i){
    if(j) --j;
    while(c[i+j]==c[sa[order[i]-1]+j]) ++j;
    heit[order[i]]=j;
    }
}
inline void work(){
    scanf("%s%s",c1+1,c2+1);
    len1=strlen(c1+1);len2=strlen(c2+1);
    for (RG int i=1;i<=len1;++i) c[i]=c1[i];
    c[len1+1]='z'+1;
    for (RG int i=1;i<=len2;++i) c[len1+1+i]=c2[i];
    /*scanf("%s",c+1);*/
    //cout << c+1 << endl;
    len=strlen(c+1);
    getsa();
    /*for (RG int i=1;i<=len;++i) cout<<order[i]<<' ';
    cout<<endl;
    for (RG int i=1;i<=len;++i) cout<<sa[i]<<' ';
    cout<<endl;
    for (RG int i=1;i<=len;++i) cout<<heit[i]<<' ';
    cout<<endl;*/
    for (RG int i=2;i<=len;++i)
    if(heit[i]>ans){
        RG int pren=sa[i-1],bacn=sa[i];
        if((pren<=len1&&bacn>=len1+2)||(bacn<=len1&&pren>=len1+2))
        ans=heit[i];
    }
    printf("%d
",ans);
    /*for(int i=1;i<=len;i++){
    for(int j=sa[i];j<=len;j++){
        cout << c[j] << " ";
    }
    cout << " " << heit[i] << endl;
    }*/
}
int main(){
    work();
    return 0;
}