bzoj3926[Zjoi2015]诸神眷顾的幻想乡 广义后缀自动机

广义后缀自动机裸题,具体可以参考15年国家队论文。
其实就是在trie上建SAM,一边dfs一边建就可以了。

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
const int N=5e5+5;
typedef long long ll;
int cnt,tot,n,m,l[N*10],head[N],col[N],d[N];
int go[N*2],next[N*2],ch[N*10][15],fa[N*10];
inline void add(int x,int y)
{
    go[++cnt]=y;
    next[cnt]=head[x];
    head[x]=cnt;
}
int ins(int p,int x)
{
    int np,q,nq;
    np=++cnt;l[np]=l[p]+1;
    for (;!ch[p][x]&&p;p=fa[p]) ch[p][x]=np;
    if (!p) fa[np]=1;
    else
    {
        q=ch[p][x];
        if (l[q]==l[p]+1) fa[np]=q;
        else
        {
            nq=++cnt;l[nq]=l[p]+1;
            memcpy(ch[nq],ch[q],sizeof(ch[q]));
            fa[nq]=fa[q];
            fa[q]=fa[np]=nq;
            for (;ch[p][x]==q;p=fa[p]) ch[p][x]=nq;
        }
    }
    return np;
}
inline void dfs(int x,int fa,int p)
{
    int t=ins(p,col[x]);
    for(int i=head[x];i;i=next[i])
    {
        int v=go[i];
        if (v!=fa)
        dfs(v,x,t);
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    fo(i,1,n)scanf("%d",&col[i]);
    fo(i,1,n-1)
    {
        int x,y;
        scanf("%d%d",&x,&y);
        add(x,y);
        add(y,x);
        d[x]++,d[y]++; 
    }
    cnt=1;
    fo(i,1,n)if (d[i]==1)dfs(i,0,1);
    ll ans=0;
    fo(i,1,cnt)ans+=l[i]-l[fa[i]];
    printf("%lld\n",ans);
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值