【bzoj2946】[Poi2000]公共串 后缀自动机

Description

   给出几个由小写字母构成的单词,求它们最长的公共子串的长度。

任务:
l 读入单词
l 计算最长公共子串的长度
l 输出结果

Input

文件的第一行是整数 n,1<=n<=5,表示单词的数量。接下来n行每行一个单词,只由小写字母组成,单词的长度至少为1,最大为2000。

Output

仅一行,一个整数,最长公共子串的长度。

Sample Input

3

abcb

bca

acbc

Sample Output

HINT

Source

.


把第一个串建后缀自动机

然后拿其它串在上面跑,记录每个串匹配到这个点的最大长度,然后在parent树上保存子树最小值(还有和val取min),然后对于所有串取最大值…

说的简单,代码…

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

const int SZ = 100010;
const int INF = 1000000010;

struct sam_node{
    sam_node *ch[30],*par;
    int val,ans[5];
}T[SZ], *root, *last;

int Tcnt = 0;

sam_node* newnode(int x)
{
    sam_node *k = T + (Tcnt ++);
    k -> val = x;
    memset(k -> ch,0,sizeof(k -> ch));
    k -> par = 0;
    memset(k -> ans,0,sizeof(k -> ans));
    return k;
}

void insert(int x)
{
    sam_node *p = last,*np = newnode(last -> val + 1);
    while(p && !p -> ch[x])
        p -> ch[x] = np,p = p -> par;
    if(!p)
        np -> par = root;
    else
    {
        sam_node *q = p -> ch[x];
        if(q -> val == p -> val + 1)
            np -> par = q;
        else
        {
            sam_node *nq = newnode(p -> val + 1);
            memcpy(nq -> ch,q -> ch,sizeof(nq -> ch));
            nq -> par = q -> par;
            q -> par = np -> par = nq;
            while(p && p -> ch[x] == q)
                p -> ch[x] = nq,p = p -> par;
        }
    }
    last = np;
}


void find(char s[],int id)
{
    int l = strlen(s);
    sam_node *p = root;
    for(int i = 0,len = 0;i < l;i ++)
    {
        int c = s[i] - 'a' + 1;
        if(p -> ch[c])
        {
            len ++;
            p = p -> ch[c];
        }
        else
        {
            while(p != root && !p -> ch[c])
                p = p -> par;
            len = p -> val;
            if(p -> ch[c])
                p = p -> ch[c],len ++;
        }
        p -> ans[id] = max(p -> ans[id],min(p -> val,len)); 
    }
    for(int i = Tcnt - 1;i >= 1;i --)
    {
        sam_node *k = T + i;
        k -> par -> ans[id] = min(k -> par -> val,max(k -> par -> ans[id],k -> ans[id]));
    }
}

void init()
{
    root = newnode(0);
    last = root;
}

int n;

int dfs(sam_node *p)
{
    int ans = INF;
    for(int i = 1;i < n;i ++)
        ans = min(ans,p -> ans[i]);
    for(int i = 1;i <= 26;i ++)
        if(p -> ch[i])
            ans = max(ans,dfs(p -> ch[i]));
    return ans;
}

char s[SZ];

int main()
{
    init();
    scanf("%d",&n);
    scanf("%s",s);
    int l = strlen(s);
    for(int i = 0;i < l;i ++)
        insert(s[i] - 'a' + 1);
    for(int i = 1;i < n;i ++)
    {
        scanf("%s",s);
        find(s,i);
    }
    printf("%d",dfs(root));
    return 0;
}

/*

3
abdab
dab
abd

*/
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值