码迷,mamicode.com
首页 > 其他好文 > 详细

3926: [Zjoi2015]诸神眷顾的幻想乡

时间:2018-03-07 17:33:21      阅读:180      评论:0      收藏:0      [点我收藏+]

标签:span   后缀自动机   none   幻想乡   printf   int   AC   last   getchar   

传送门

一个广义后缀自动机模板。

技术分享图片
//Achen
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<vector>
#include<cstdio>
#include<queue>
#include<cmath>
const int N=4000007;
typedef long long LL;
using namespace std;
int n,c,col[N];

template<typename T>void read(T &x)  {
    char ch=getchar(); x=0; T f=1;
    while(ch!=-&&(ch<0||ch>9)) ch=getchar();
    if(ch==-) f=-1,ch=getchar();
    for(;ch>=0&&ch<=9;ch=getchar()) x=x*10+ch-0; x*=f;
}

int ecnt,fir[N],nxt[N],to[N],in[N];
void add(int u,int v) {
    nxt[++ecnt]=fir[u]; fir[u]=ecnt; to[ecnt]=v; in[v]++;
    nxt[++ecnt]=fir[v]; fir[v]=ecnt; to[ecnt]=u; in[u]++;
}

int tot,rt,p,np,ch[N][11],fa[N],sz[N],l[N];
int insert(int c,int last) {
    p=last; np=++tot; 
    l[np]=l[p]+1;
    for(;p&&!ch[p][c];p=fa[p]) ch[p][c]=np;
    if(!p) fa[np]=rt;
    else {
        int q=ch[p][c];
        if(l[p]+1==l[q]) fa[np]=q;
        else {
            int nq=++tot; l[nq]=l[p]+1;
            memcpy(ch[nq],ch[q],sizeof(ch[q]));
            fa[nq]=fa[q]; fa[q]=fa[np]=nq;
            for(;p&&ch[p][c]==q;p=fa[p]) ch[p][c]=nq;
        }
    }
    sz[np]++; return np;
}

void dfs(int x,int fa,int pa) {
    int last=insert(col[x],pa);
    for(int i=fir[x];i;i=nxt[i]) if(to[i]!=fa) {
        dfs(to[i],x,last);
    }
}

int main() {
    read(n); read(c); rt=++tot;
    for(int i=1;i<=n;i++) read(col[i]);
    for(int i=1;i<n;i++) {
        int u,v;
        read(u); read(v);
        add(u,v);
    }
    for(int i=1;i<=n;i++) if(in[i]==1) 
        dfs(i,0,rt);
    LL ans=0;
    for(int i=2;i<=tot;i++) ans+=(LL)l[i]-l[fa[i]];
    printf("%lld\n",ans); 
    return 0;
}
View Code

 

3926: [Zjoi2015]诸神眷顾的幻想乡

标签:span   后缀自动机   none   幻想乡   printf   int   AC   last   getchar   

原文地址:https://www.cnblogs.com/Achenchen/p/8523595.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!