后缀数组 LCP--模板题

题意：

给你S串和T串，用T串的所有前缀去匹配S串（匹配值是最长公共子串）。

问你总值相加是多少。

思路：

先把两个S，T串倒过来，再拼接 S#T 合成一串，跑一下后缀数组

在排序好的rank里计算每个T后缀的最长匹配长度。（前后两个for即可）

最后dp对后缀取max，累计答案。（因为后缀从pos开始的ans1肯定被后缀从pos-1开始的ans2包含，所以如果ans2<ans1，那要对ans2取max（ans1）

#define IOS ios_base::sync_with_stdio(0); cin.tie(0);
#include <cstdio>//sprintf islower isupper
#include <cstdlib>//malloc  exit strcat itoa system("cls")
#include <iostream>//pair
#include <fstream>//freopen("C:\Users\13606\Desktop\草稿.txt","r",stdin);
#include <bitset>
//#include <map>
//#include<unordered_map>   https://www.nitacm.com/problem_show.php?pid=585
#include <vector>
#include <stack>
#include <set>
#include <string.h>//strstr substr
#include <string>
#include <time.h>//srand(((unsigned)time(NULL))); Seed n=rand()%10 - 0~9;
#include <cmath>
#include <deque>
#include <queue>//priority_queue<int, vector<int>, greater<int> > q;//less
#include <vector>//emplace_back
//#include <math.h>
//#include <windows.h>//reverse(a,a+len);// ~ ! ~ ! floor
#include <algorithm>//sort + unique : sz=unique(b+1,b+n+1)-(b+1);+nth_element(first, nth, last, compare)
using namespace std;//next_permutation(a+1,a+1+n);//prev_permutation
#define fo(a,b,c) for(register int a=b;a<=c;++a)
#define fr(a,b,c) for(register int a=b;a>=c;--a)
#define mem(a,b) memset(a,b,sizeof(a))
#define pr printf
#define sc scanf
#define ls rt<<1
#define rs rt<<1|1
typedef long long ll;
#define rint register int
void swapp(int &a,int &b);
double fabss(double a);
int maxx(int a,int b);
int minn(int a,int b);
int Del_bit_1(int n);
int lowbit(int n);
int abss(int a);
//const long long INF=(1LL<<60);
const double E=2.718281828;
const double PI=acos(-1.0);
const int inf=(1<<30);
const double ESP=1e-9;
const int mod=(int)1e9+7;
const int N=(int)2e6+10;
void PR(int _[],int n)
{
    for(int i=1;i<=n;++i)
        pr("%d ",_[i]);
    pr("
");
}
char s[N],s1[N],s2[N];
int a[N];
int sa[N],rk[N],s_a[N],t[N];//桶的大小要为N,因为放的是rank;
int height[N],h[N];//h是位置i的长度,heidgt是ranki的长度;
void Init(int _[],int n)
{
    for(rint i=0;i<=n;++i)
        _[i]=0;
}
bool cmp(int i,int j,int k)
{
    return s_a[i]==s_a[j]&&s_a[i+k]==s_a[j+k];
}
void Sort(int len)
{
    int m=30;//字符集大小;
    for(rint i=1;i<=len;++i) ++t[a[i]],rk[i]=a[i];
    for(rint i=1;i<=m;++i) t[i]+=t[i-1];
    for(rint i=len;i>=1;--i) sa[t[rk[i]]--]=i;
    for(rint k=1;k<=len;k<<=1)
    {
        int cnt=0;
        //按第二个rank排;
        for(rint i=len-k+1;i<=len;++i) s_a[++cnt]=i;
        for(rint i=1;i<=len;++i)if(sa[i]>k) s_a[++cnt]=sa[i]-k;
        //按第一个rank排;
        Init(t,m);
        for(rint i=1;i<=len;++i) ++t[rk[s_a[i]]];
        for(rint i=1;i<=m;++i) t[i]+=t[i-1];
        for(rint i=len;i>=1;--i) sa[t[rk[s_a[i]]]--]=s_a[i];

        swap(rk,s_a);rk[sa[1]]=cnt=1;
        for(rint i=2;i<=len;++i)
            rk[sa[i]]=cmp(sa[i],sa[i-1],k)?cnt:++cnt;
        if(cnt==len)break;
        m=cnt;
    }
    //求height数组;
    for(rint i=1;i<=len;++i)
    {
        h[i]=max(0,h[i-1]-1);
        if(rk[i]==1)continue;
        while(a[i+h[i]]==a[sa[rk[i]-1]+h[i]]) ++h[i];
    }
    for(rint i=1;i<=len;++i) height[i]=h[sa[i]];
}

int ans[N],res[N];

int main()
{
//    freopen("D:\Chrome Download\testdata (2).in","r",stdin);
    int len1,len2;
    sc("%d%d",&len1,&len2);
    sc("%s%s",s1+1,s2+1);
    int len=0;
    for(int i=len1;i>=1;--i)
        a[++len]=s1[i]-'a'+1;
    a[++len]=27;
    for(int i=len2;i>=1;--i)
        a[++len]=s2[i]-'a'+1;
    Sort(len);
//    PR(height,len);
//    PR(sa,len);
//    PR(rk,len);
    int temp=0;
    for(int i=1;i<=len;++i)
    {
    //    temp=min(temp,height[i+1]);
        if(sa[i]<=len1)
            temp=height[i+1];
        else
            ans[i]=temp,temp=min(temp,height[i+1]);
    }
    temp=0;
    for(int i=len;i>=1;--i)
    {
    //    temp=min(temp,height[i+1]);
        if(sa[i]<=len1)
            temp=height[i];
        else
            ans[i]=max(ans[i],temp),temp=min(temp,height[i]);
    }
//    PR(ans,len);
    ll Ans=0;
    for(int i=len;i>=len1+2;--i)
    {
        res[i]=max(res[i+1],ans[rk[i]]);
        Ans+=res[i];
    }
    pr("%lld
",Ans);
    return 0;
}

/**************************************************************************************/

int maxx(int a,int b)
{
    return a>b?a:b;
}

void swapp(int &a,int &b)
{
    a^=b^=a^=b;
}

int lowbit(int n)
{
    return n&(-n);
}

int Del_bit_1(int n)
{
    return n&(n-1);
}

int abss(int a)
{
    return a>0?a:-a;
}

double fabss(double a)
{
    return a>0?a:-a;
}

int minn(int a,int b)
{
    return a<b?a:b;
}