BZOJ3926&&lg3346 ZJOI诸神眷顾的幻想乡(广义后缀自动机)

BZOJ3926&&lg3346 ZJOI诸神眷顾的幻想乡(广义后缀自动机)

题面

自己找去

HINT

我们可以把题目拆解成几个部分，首先我们手玩一个结论，从所有的叶子节点出发，遍历整颗树，这里面一定包括了所有的文本串，既然如此，就上广义后缀自动机维护，然后本质不同子串个数，也就是\(\sum_{u=2}^{tot}{node[u].len-node[fa].len}\)，然后这题就OK了

#include<bits/stdc++.h>

using namespace std;

inline int read(){

    int w=0,f=1;

    char ch=getchar();

    while(ch<'0'||ch>'9'){

        if(ch=='-')f=-1;

        ch=getchar();

    }

    while(ch>='0'&&ch<='9'){

        w=(w<<3)+(w<<1)+ch-48;

        ch=getchar();

    }

    return w*f;

}

int n,m,val[2000010];

struct SUFFIXAUTOMATON{

    struct Node{

        int len,fa;

        int ch[20];

    }node[2000010];

    int lst,tot,root;

    inline void init(){

        lst=tot=root=1;return;

    }

    inline void extend(int now){

        int p=lst;tot++;lst=tot;int np=tot;

        node[np].len=node[p].len+1;

        while(p&&!node[p].ch[now]){

            node[p].ch[now]=np;

            p=node[p].fa;

        }

        if(!p) node[np].fa=1;

        else{

            int q=node[p].ch[now];

            if(node[q].len==node[p].len+1){

                node[np].fa=q;

            }

            else{

                int nq=++tot;node[nq]=node[q];

                node[nq].len=node[p].len+1;

                node[q].fa=nq;node[np].fa=nq;

                while(p&&node[p].ch[now]==q){

                    node[p].ch[now]=nq;

                    p=node[p].fa;

                }

            }

        }

    }

}SAM;

int cnt,head[2000010],du[2000010];

struct Edge{

    int from,to,next;

}edge[4000010];

inline void addedge(int u,int v){

    cnt++;

    edge[cnt].from=u;

    edge[cnt].to=v;

    edge[cnt].next=head[u];

    head[u]=cnt;du[u]++;

}

inline void dfs(int u,int fa){

    SAM.extend(val[u]);int cur=SAM.lst;

    for(int i=head[u];i;i=edge[i].next){

        int v=edge[i].to;if(v==fa) continue;

        SAM.lst=cur;dfs(v,u);

    }

    return;

}

int main(){

    n=read();m=read();SAM.init();

    for(int i=1;i<=n;i++){

        val[i]=read();val[i]++;

    }

    for(int i=1;i<n;i++){

        int x=read();int y=read();

        addedge(x,y);addedge(y,x);

    }

    for(int i=1;i<=n;i++){

        if(du[i]!=1) continue;

        SAM.lst=1;dfs(i,0);

    }

    long long ans=0;

    for(int i=1;i<=SAM.tot;i++){

        int f=SAM.node[i].fa;

        ans+=SAM.node[i].len-SAM.node[f].len;

    }

    cout<<ans<<endl;

    return 0;

}