#4214. 谢特

题意

内存限制:512 MiB
时间限制:1000 ms

由于你成功地在 1 s \text{1 s} 1 s 内算出了上一题的答案,英雄们很高兴并邀请你加入了他们的游戏。然而进入游戏之后你才发现,英雄们打的游戏和你想象的并不一样……

英雄们打的游戏是这样的:首先系统会产生(注意不一定是随机产生)一个字符串,然后每个英雄就会开始根据自己分到的任务计算这个字符串的某些特征,谁先算出自己的答案谁就是胜者。

由于打游戏的英雄比较多,因此英雄们分到的任务也就可能很奇怪。比如你分到的这个任务就是这样:

定义这个字符串以第 i i i 个字符开头的后缀为后缀 i i i (编号从 1 1 1 开始),每个后缀 i i i 都有一个权值 w i w_i wi ,同时定义两个后缀 i , j i,j i,j ( i ≠ j i \ne j i̸=j) 的贡献为它们的最长公共前缀长度加上它们权值的异或和,也就是 L C P ( i , j ) + ( w i xor w j ) \mathrm{LCP}(i,j)+(w_i \mathbin{\text{xor}} w_j) LCP(i,j)+(wixorwj) 。而你的任务就是,求出这个字符串的所有后缀两两之间贡献的最大值。

n ≤ 1 0 5 n\le 10^5 n105 w i &lt; n w_i&lt; n wi<n

题解

考虑反串,建立 S A M \mathrm{SAM} SAM ,可以考虑到 p a r e n t parent parent 树上每个节点所能接受的最长的长度即为这个节点的子树内的所有原串中的后缀的 L C P \mathrm{LCP} LCP

所以每个点建立 t r i e trie trie 树,在 t r i e trie trie 树上寻找最优解,然后 t r i e trie trie 树合并即可

#include <bits/stdc++.h>
#define I inline
using namespace std;
const int N=2e5+5;char s[N];
int n,lst=1,sz=1,w[N],ans,hd[N],V[N];
int nx[N],t,tt,T[N],ch[N*20][2];
struct SAM{
	int link,len;
	map<int,int>nx;
}a[N];
I void build(int x){
	int p=lst,np=++sz;
	a[np].len=a[p].len+1;
	while(p && !a[p].nx.count(x))
		a[p].nx[x]=np,p=a[p].link;
	if (!p) a[np].link=1;
	else{
		int q=a[p].nx[x];
		if (a[q].len==a[p].len+1)
			a[np].link=q;
		else{
			int nq=++sz;a[nq]=a[q];
			a[nq].len=a[p].len+1;
			a[np].link=a[q].link=nq;
			while(p && a[p].nx[x]==q)
				a[p].nx[x]=nq,p=a[p].link;
		}
	}
	lst=np;
}
I void add(int u,int v){V[++t]=v;nx[t]=hd[u];hd[u]=t;}
I void ins(int &x,int v,int d){
	if (!x) x=++tt;
	if (~d) ins(ch[x][(v>>d)&1],v,d-1);
}
I int query(int x,int y,int d){
	if (!~d) return 0;
	if (ch[x][0] && ch[x][1]){
		if (ch[y][0] && ch[y][1])
			return max(query(ch[x][0],ch[y][1],d-1),query(ch[x][1],ch[y][0],d-1))|(1<<d);
		else{
			bool v=ch[y][1];
			return query(ch[x][v^1],ch[y][v],d-1)|(1<<d);
		}
	}
	else{
		bool v=ch[x][1];
		return ch[y][v^1]?(query(ch[x][v],ch[y][v^1],d-1)|(1<<d)):query(ch[x][v],ch[y][v],d-1);
	}
}
I int merge(int x,int y){
	if (!x || !y) return x|y;int z=++tt;
	ch[z][0]=merge(ch[x][0],ch[y][0]);
	ch[z][1]=merge(ch[x][1],ch[y][1]);
	return z;
}
I void dfs(int x){
	for (int v,i=hd[x];i;i=nx[i]){
		dfs(v=V[i]);
		if (T[x] && T[v])
			ans=max(ans,a[x].len+query(T[x],T[v],18));
		T[x]=merge(T[x],T[v]);
	}
}
int main(){
	scanf("%d%s",&n,s+1);
	for (int i=1;i<=n;i++) scanf("%d",&w[i]);
	for (int i=n;i;i--) build(s[i]),ins(T[lst],w[i],18);
	for (int i=2;i<=sz;i++) add(a[i].link,i);dfs(1);
	return printf("%d\n",ans),0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值