【spoj1811 & spoj1812 - LCS1 & LCS2】sam

spoj1811  给两个长度小于100000的字符串 A 和 B，求出他们的最长公共连续子串。

先将串 A 构造为 SAM ，然后用 B 按如下规则去跑自动机。
用一个变量 lcs 记录当前的最长公共子串，初始化为0。
设当前状态结点为 p，要匹配的字符为 c，若 go[c] 中有边，说明能够转移状态，则转移并 lcs++；
若不能转移则将状态移动到 p 的 par ，如果仍然不能转移则重复该过程直到 p 回到根节点，并将 lcs 置为 0；
如果在上一个过程中进入了能够转移的状态，则设 lcs 为当前状态的 val。
为什么失配后要移向 par 呢？因为在状态 p 上失配说明该状态的 [min,max] 所表示字符串都不是 B 中的子串，但是比它们短的后缀仍有可能是 B 的子串，而 par 指针恰好指向了该状态的后缀。

spoj1812  要求多个串的最长公共子串,n<=100000,10个串。

本来想用广义sam，然后对每个点位压10位表示它能到哪些串，最后dfs一遍。
20*10^5=10^6个点，时间压得很紧直接超时。

参照上题做法，就把一个串A放进去，nlcs[x]维护当前串与A在自动机的x节点匹配的最大长度，lcs[x]维护前i个串。做完一个串后for一遍，x没有走到过的点lcs[x]=0，其余的lcs[x]=minn(lcs[x],nlcs[x]);

spoj1811

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;

const int N=*;
char a[N],b[N];
int al,bl,tot,last,son[N][],step[N],pre[N];

int maxx(int x,int y){return x>y ? x:y;}

int add_node(int x)
{
    step[++tot]=x;
    return tot;
}

void extend(int ch)
{
    int p=last,np=add_node(step[p]+);
    while(p && !son[p][ch]) son[p][ch]=np,p=pre[p];
    if(!p) pre[np]=;
    else
    {
        int q=son[p][ch];
        if(step[q]==step[p]+) pre[np]=q;
        else
        {
            int nq=add_node(step[p]+);
            memcpy(son[nq],son[q],sizeof(son[q]));
            pre[nq]=pre[q];
            pre[q]=pre[np]=nq;
            while(son[p][ch]==q) son[p][ch]=nq,p=pre[p];
        }
    }
    last=np;
}

int main()
{
    freopen("a.in","r",stdin);
    // freopen("me.out","w",stdout);
    memset(son,,sizeof(son));
    memset(pre,,sizeof(pre));
    memset(step,,sizeof(step));
    tot=;add_node();last=;
    scanf("%s%s",a+,b+);
    al=strlen(a+),bl=strlen(b+);
    for(int i=;i<=al;i++) extend(a[i]-'a'+);
    int ch,x=,ans=,len=;
    for(int i=;i<=bl;i++)
    {
        ch=b[i]-'a'+;
        while(x && !son[x][ch]) x=pre[x],len=step[x];
        if(x==) x=;
        if(son[x][ch]) x=son[x][ch],len++;

        ans=maxx(ans,len);
    }
    printf("%d\n",ans);
    return ;
}

spoj1812

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;

const int N=*;
char s[][N];
int n,tot,last,son[N][],step[N],pre[N],lcs[N],nlcs[N],l[];
bool vis[N];
struct node{
    int x,y,next;
}a[*N];

int maxx(int x,int y){return x>y ? x:y;}
int minn(int x,int y){return x<y ? x:y;}

int add_node(int x)
{
    step[++tot]=x;lcs[tot]=x;
    return tot;
}

void extend(int ch)
{
    int p=last,np=add_node(step[p]+);
    while(p && !son[p][ch]) son[p][ch]=np,p=pre[p];
    if(!p) pre[np]=;
    else
    {
        int q=son[p][ch];
        if(step[q]==step[p]+) pre[np]=q;
        else
        {
            int nq=add_node(step[p]+);
            memcpy(son[nq],son[q],sizeof(son[q]));
            pre[nq]=pre[q];
            pre[q]=pre[np]=nq;
            while(son[p][ch]==q) son[p][ch]=nq,p=pre[p];
        }
    }
    last=np;
}

void solve(int id)
{
    memset(nlcs,,sizeof(nlcs));
    memset(vis,,sizeof(vis));
    int x=,len=,ch;
    for(int i=;i<=l[id];i++)
    {
        ch=s[id][i]-'a'+;
        while(x && !son[x][ch])
        {
            x=pre[x],len=step[x],vis[x]=;
            nlcs[x]=maxx(nlcs[x],len);
        }
        if(x==) x=,vis[]=;
        if(son[x][ch]) x=son[x][ch],vis[x]=,len++;
        nlcs[x]=maxx(nlcs[x],len);
    }
    for(int i=;i<=tot;i++)
        if(!vis[i]) lcs[i]=;
        else lcs[i]=minn(lcs[i],nlcs[i]);
}

int main()
{
    freopen("a.in","r",stdin);
    // freopen("me.out","w",stdout);
    memset(son,,sizeof(son));
    memset(pre,,sizeof(pre));
    memset(step,,sizeof(step));
    tot=;add_node();last=;n=;
    while(scanf("%s",s[++n]+)!=EOF)
    {
        l[n]=strlen(s[n]+);
    }n--;
    for(int i=;i<=l[];i++) extend(s[][i]-'a'+);
    for(int i=;i<=n;i++) solve(i);
    int ans=;
    for(int i=;i<=tot;i++) ans=maxx(ans,lcs[i]);
    printf("%d\n",ans);
    return ;
}