2019.8.28后缀数组入门...

8月23的CCPC网赛出了神奇的后缀数组。吓死本宝宝了。
啃了五天板子，还是没太看懂。
那么先用起来，再慢慢体会吧。

两个串的最长公共连续子串

POJ-2774 Long Long Message

HDOJ-1403 Longest Common Substring

题意

给你两个字符串，要求输出最长公共子串长度。

思路

参见大佬的博客
字符串的任何一个子串都是这个字符串的某个后缀的前缀。
求A和B的最长公共子串等价于求A的后缀和B的后缀的最长公共前缀的最大值。
将AB串拼接起来，中间用一个没有出现过的字符隔开。
然后遍历height数组，通过sa[i-1]与sa[i]与插入的隔板字符相对关系判断最长前缀是否来自不同的串。

AC的代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<vector>
#include<queue>
#include<stack>
#include<map>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int maxn = 2e6+10;
struct SuffixArray
{
    int s[maxn];      // 原始字符数组（最后一个字符应必须是0，而前面的字符必须非0）
    int sa[maxn];     // 后缀数组.下标1:n-1有效
    int rak[maxn];    // 名次数组. rak[0]一定是n-1，即最后一个字符.下标1:n-1有效
    int height[maxn]; // height数组.下标2:n-1有效
    int t[maxn], t2[maxn], c[maxn]; // 辅助数组
    int n; // 字符个数,n为ch长度+1
    void clear()
    {
        n=0;
        memset(sa,0,sizeof(sa));
    }
    // m为最大字符值加1。调用之前需设置好s和n
    void build_sa(int m)
    {
        int i, *x = t, *y = t2;
        for(i = 0; i < m; i++) c[i]=0;
        for(i = 0; i < n; i++) c[x[i]=s[i]]++;
        for(i = 1; i < m; i++) c[i]+=c[i-1];
        for(i = n-1; i >= 0; i--) sa[--c[x[i]]]=i;
        for(int k = 1; k <= n; k <<= 1)
        {
            int p = 0;
            for(i = n-k; i < n; i++) y[p++] = i;
            for(i = 0; i < n; i++) if(sa[i] >= k) y[p++] = sa[i]-k;

            for(i = 0; i < m; i++) c[i] = 0;
            for(i = 0; i < n; i++) c[x[y[i]]]++;
            for(i = 1; i < m; i++) c[i] += c[i-1];
            for(i = n-1; i >= 0; i--) sa[--c[x[y[i]]]] = y[i];

            swap(x, y);
            p = 1;
            x[sa[0]] = 0;
            for(i = 1; i < n; i++)
                x[sa[i]] = y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k] ? p-1 : p++;
            if(p >= n) break;
				m = p;
        }
    }
    void build_height()
    {
        int i, j, k = 0;
        for(i = 0; i < n; i++)
			rak[sa[i]] = i;
        for(i = 0; i < n; i++)
        {
            if(k) k--;
            int j = sa[rak[i]-1];
            while(s[i+k] == s[j+k]) k++;
            height[rak[i]] = k;
        }
    }
    void debug()
	{
		for(int i = 1; i < n ; i ++)
		{
			printf("%d ", sa[i]);
			printf("%d\n", height[i]);
		}puts("");
	}
};

SuffixArray sa;
char ch[maxn];

int main()
{
	int len,bar,ans=0;
	while(~scanf("%s",ch))
	{
		bar=strlen(ch);
		ch[bar]='{';
		scanf("%s",ch+bar+1);
		len=strlen(ch);

		sa.clear();
		for(int i=0;i<len;i++)
			sa.s[i]=ch[i]-'a'+1;
		sa.s[len]=0;
		sa.n=len+1;
		sa.build_sa(30);
		sa.build_height();

		ans=0;
		for(int i=1;i<=len;i++)
			if(sa.height[i]>ans&&(ll)(sa.sa[i]-bar)*(sa.sa[i-1]-bar)<0)
				ans=sa.height[i];
//		for(int i=1;i<=len;i++)
//		{
//			cout<<"排名="<<sa.rak[sa.sa[i]]<<",hei="<<sa.height[i]<<"\t:"<<ch+sa.sa[i]<<endl;
//		}
//		sa.debug();
		cout<<ans<<endl;
	}
	return 0;
}

不可重叠最长重复子串

POJ-1743 Musical Theme