[BZOJ3926][ZJOI2015]诸神眷顾的幻想乡(后缀自动机)

 

日,无数幽香的粉丝到了幽香家门前的太阳花田上来为幽香庆祝生日。 

粉丝们非常热情,自发组织表演了一系列节目给幽香看。幽香当然也非常高兴啦。 
这时幽香发现了一件非常有趣的事情,太阳花田有n块空地。在过去,幽香为了方便,在这n块空地之间修建了n-1条边将它们连通起来。也就是说,这n块空地形成了一个树的结构。 
有n个粉丝们来到了太阳花田上。为了表达对幽香生日的祝贺,他们选择了c中颜色的衣服,每种颜色恰好可以用一个0到c-1之间的整数来表示。并且每个人都站在一个空地上,每个空地上也只有一个人。这样整个太阳花田就花花绿绿了。幽香看到了,感觉也非常开心。 
粉丝们策划的一个节目是这样的,选中两个粉丝A和B(A和B可以相同),然后A所在的空地到B所在的空地的路径上的粉丝依次跳起来(包括端 点),幽香就能看到一个长度为A到B之间路径上的所有粉丝的数目(包括A和B)的颜色序列。一开始大家打算让人一两个粉丝(注意:A,B和B,A是不同 的,他们形成的序列刚好相反,比如红绿蓝和蓝绿红)都来一次,但是有人指出这样可能会出现一些一模一样的颜色序列,会导致审美疲劳。 
于是他们想要问题,在这个树上,一共有多少可能的不同的颜色序列(子串)幽香可以看到呢? 
太阳花田的结构比较特殊,只与一个空地相邻的空地数量不超过20个。 

Input

 第一行两个正整数n,c。表示空地数量和颜色数量。 

第二行有n个0到c-1之间,由空格隔开的整数,依次表示第i块空地上的粉丝的衣服颜色。(这里我们按照节点标号从小到大的顺序依次给出每块空地上粉丝的衣服颜色)。 
接下来n-1行,每行两个正整数u,v,表示有一条连接空地u和空地v的边。 

Output

 一行,输出一个整数,表示答案。 

Sample Input

7 3
0 2 1 2 1 0 0
1 2
3 4
3 5
4 6
5 7
2 5

Sample Output

30

HINT

对于所有数据,1<=n<=100000, 1<=c<=10。 


对于15%的数据,n<=2000。 

另有5%的数据,所有空地都至多与两个空地相邻。 

另有5%的数据,除一块空地与三个空地相邻外,其他空地都分别至多与两个空地相邻。 

另有5%的数据,除某两块空地与三个空地相邻外,其他空地都分别至多与两个空地相邻

Source

[Submit][Status][Discuss]

SAM模板题,只要看出来字符串模型就行。注意要从叶子开始。

UPD:哦这好像是个广义后缀自动机,没什么区别吧。

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<algorithm>
 4 #define rep(i,l,r) for (int i=l; i<=r; i++)
 5 typedef long long ll;
 6 using namespace std;
 7 
 8 const int N=200100,M=4000100;
 9 ll ans;
10 int tot=1,n,c,u,v,cnt,ind[N],a[N],son[M][11],fa[M],mx[M],nxt[N],h[N],to[N];
11 void add(int u,int v){ to[++cnt]=v; nxt[cnt]=h[u]; h[u]=cnt; }
12 
13 int extend(int p,int c){
14     int np=++tot; mx[np]=mx[p]+1;
15     while (!son[p][c] && p) son[p][c]=np,p=fa[p];
16     if (!p) fa[np]=1;
17     else{
18         int q=son[p][c];
19         if (mx[q]==mx[p]+1) fa[np]=q;
20         else{
21             int nq=++tot; mx[nq]=mx[p]+1;
22             memcpy(son[nq],son[q],sizeof(son[q]));
23             fa[nq]=fa[q]; fa[np]=fa[q]=nq;
24             while (p && son[p][c]==q) son[p][c]=nq,p=fa[p];
25         }
26     }
27     return np;
28 }
29 
30 void solve(){ rep(i,1,tot) ans+=mx[i]-mx[fa[i]]; }
31 
32 void dfs(int x,int fa,int p){
33     int t=extend(p,a[x]);
34     for (int i=h[x]; i; i=nxt[i]) if (to[i]!=fa) dfs(to[i],x,t);
35 }
36 
37 int main(){
38     freopen("bzoj3926.in","r",stdin);
39     freopen("bzoj3926.out","w",stdout);
40     scanf("%d%d",&n,&c);
41     rep(i,1,n) scanf("%d",&a[i]);
42     rep(i,1,n-1) scanf("%d%d",&u,&v),add(u,v),add(v,u),ind[u]++,ind[v]++;
43     rep(i,1,n) if (ind[i]==1) dfs(i,0,1);
44     solve(); printf("%lld\n",ans);
45     return 0;
46 }

 

posted @ 2018-03-22 08:52  HocRiser  阅读(206)  评论(0编辑  收藏  举报