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

最新推荐文章于 2022-09-19 11:18:35 发布

我要吃熊猫

最新推荐文章于 2022-09-19 11:18:35 发布

阅读量364

点赞数

分类专栏：后缀自动机

本文链接：https://blog.csdn.net/chai_jing/article/details/72187513

版权

后缀自动机专栏收录该内容

20 篇文章 0 订阅

订阅专栏

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

Time Limit: 10 Sec Memory Limit: 512 MB
Submit: 1381 Solved: 811
[ Submit][ Status][ Discuss]

Description

幽香是全幻想乡里最受人欢迎的萌妹子，这天，是幽香的2600岁生日，无数幽香的粉丝到了幽香家门前的太阳花田上来为幽香庆祝生日。

粉丝们非常热情，自发组织表演了一系列节目给幽香看。幽香当然也非常高兴啦。

这时幽香发现了一件非常有趣的事情，太阳花田有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

HINT

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

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

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

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

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

Source

【分析】

摔。为什么大家写的广义后缀自动机都不一样啊QAQ

【代码】

//bzoj 诸神眷顾的幻想乡 
#include<iostream>
#include<cstring>
#include<cstdio>
#define ll long long
#define M(a) memset(a,0,sizeof a)
#define fo(i,j,k) for(i=j;i<=k;i++)
using namespace std;
const int mxn=4000005;
ll ans;
int n,m,len,tot,cnt,root;
int head[mxn],du[mxn],s[mxn],step[mxn],pre[mxn],son[mxn][10];
struct edge {int to,next;} f[mxn]; 
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9') {if(ch=='-') f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9') x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
	return x*f;
}
inline void add(int u,int v)
{
	f[++cnt].to=v,f[cnt].next=head[u],head[u]=cnt;
}
inline int sam(int p,int c)
{
	int np=(++tot);
	step[np]=step[p]+1;
	while(p && !son[p][c])
	  son[p][c]=np,p=pre[p];
	if(!p) {pre[np]=root;return np;}
	int q=son[p][c];
	if(step[q]==step[p]+1) pre[np]=q;
	else
	{
		int nq;
		step[nq=(++tot)]=step[p]+1;
		memcpy(son[nq],son[q],sizeof son[q]);
		pre[nq]=pre[q];
		pre[np]=pre[q]=nq;
		while(p && son[p][c]==q)
		  son[p][c]=nq,p=pre[p];
	}
	return np;
}
inline void dfs(int u,int fa,int p)
{
	p=sam(p,s[u]);
	for(int i=head[u];i;i=f[i].next)
	{
		int v=f[i].to;
		if(v!=fa) dfs(v,u,p);
	}
}
int main()
{
	int i,j,u,v;
	n=read(),m=read();
	fo(i,1,n) s[i]=read();
	fo(i,2,n)
	{
		u=read(),v=read();
		add(u,v),add(v,u);
		du[u]++,du[v]++;
	}
	tot=root=1;
	fo(i,1,n) if(du[i]==1) dfs(i,0,1);
	fo(i,2,tot) ans+=step[i]-step[pre[i]];
	printf("%lld\n",ans);
	return 0;
}