【后缀自动机】SPOJ 1812-LCSII

题意：

　　给出最多10个长度不超过100000的字符串，求他们的LCS的长度。时限是鬼畜的0.25s 。

 

后缀自动机练习...虽然有人这么说但我并不觉得hash能过。

 

本题可以说是【论SAM中按step排序更新pre的重要性】：

　　总的来说做法和1811-LCS有点类似，不同的是因为有多个字符串，因此每个字符串都需要在SAM上跑一次。

　　记录下每个节点最长能容纳长度为多少的字符，最后取个MAX就行了。

 

　　用nans[i]表示匹配当前字符串时，i号点能容纳的子串长度，如果i的pre上的当前答案比i还要小的话就用nans[i]来更新nans[i.pre] 。

　　用ans表示所有nans的最小值，在计算每个字符串时用nans来更新ans，最后在ans中找一个MAX就可以了。

 

一开始忘记了黑体的那句话，结果拍了很久又没拍出来，提交的时候在第10个点WA了简直难过...

 

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <iostream>
#define maxn 105000*2
using namespace std;

struct node { int s[26],pre; } Suf[maxn];
char S[maxn],T[maxn];
int i,j,n,m,k,len1,len2,len,last = 1,root = 1,cnt = 1;
int step[maxn],nans[maxn],ans[maxn],f[maxn];

void add(int nxt)
{
    int now = ++cnt, p = last;
    step[now] = step[p] + 1;
    for ( ; p != 0 && Suf[p].s[nxt] == 0; p = Suf[p].pre) Suf[p].s[nxt] = now;
    if (p == 0) Suf[now].pre = root;
    else
    {
        int q = Suf[p].s[nxt];
        if (step[q] == step[p] + 1) Suf[now].pre = q;
        else
        {
            int nq = ++cnt;
            step[nq] = step[p] + 1;
            Suf[nq] = Suf[q];
            Suf[q].pre = Suf[now].pre = nq;
            for ( ; p != 0 && Suf[p].s[nxt] == q; p = Suf[p].pre) Suf[p].s[nxt] = nq;
        }
    }
    last = now;
}

bool cmp(int a,int b) { return step[a] < step[b]; }

int main()
{
    scanf("%s\n",S+1);
    len1 = strlen(S+1);
    for (i = 1; i <= len1; i++) add(S[i]-'a');
    
    for (i = 1; i <= cnt; i++) f[i] = i;
    stable_sort(f+1,f+1+cnt,cmp);

    memcpy(ans,step,sizeof(ans));
    
    while (scanf("%s\n",T+1) != EOF)
    {
        memset(nans,0,sizeof(nans));
        
        len2 = strlen(T+1);
        int now = root,len = 0;
        for (i = 1; i <= len2; i++)
        {
            int nxt = T[i]-'a';
            if (Suf[now].s[nxt] != 0) 
            {
                len++;
                now = Suf[now].s[nxt];
                nans[now] = max(nans[now],len);
            }
            else
            {
                while (Suf[now].s[nxt] == 0 && now != 0) now = Suf[now].pre;
                if (now == 0) now = root,len = 0;
                else len = step[now] + 1,now = Suf[now].s[nxt],nans[now] = max(nans[now],len);
            }
        }
        for (i = cnt; i >= 1; i--)
        {
            ans[f[i]] = min(ans[f[i]],nans[f[i]]);
            nans[Suf[f[i]].pre] = max(nans[Suf[f[i]].pre],nans[f[i]]);
        }
    }
    
    int ANS = 0;
    for (i = 1; i <= cnt; i++)
        ANS = max(ANS,ans[i]);
    
    printf("%d\n",ANS);
    return 0;
}