LOJ6198 谢特

题目传送门

分析:
题目名称好评!!!(错乱
考虑两个串的LCP,一一对应整个串反过来的LCS,两个串的LCS是其在SAM的Parent树上的LCA
我们把思路指向SAM,在Parent树上求答案
如何求两两异或最大值,Trie树遍历就可以了
每次将儿子的(w),合并到父亲上,合并之前询问两点联通块最大值,LCP是父亲节点的len
合并的时候启发式一下,复杂度(O(nlog^{2}n))

#include<cstdio>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>

#define maxn 200005
#define INF 0x3f3f3f3f

using namespace std;

inline int getint()
{
	int num=0,flag=1;char c;
	while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
	while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
	return num*flag;
}

int n;
struct sam{
	int nxt[26];
	int fa,len;
}t[maxn];
char s[maxn];
int lst=1,cnt=1;
int ch[maxn<<5][2],w[maxn];
int rt[maxn],tot;
vector<int>V[maxn];
int pos[maxn],b[maxn],ans;

inline void insert(int c)
{
	int p=lst,np=lst=++cnt;
	t[np].len=t[p].len+1;
	while(!t[p].nxt[c]&&p)t[p].nxt[c]=np,p=t[p].fa;
	if(!p)t[np].fa=1;
	else
	{
		int q=t[p].nxt[c];
		if(t[q].len==t[p].len+1)t[np].fa=q;
		else
		{
			int nq=++cnt;
			memcpy(t[nq].nxt,t[q].nxt,sizeof t[q].nxt);
			t[nq].len=t[p].len+1;
			t[nq].fa=t[q].fa;
			t[q].fa=t[np].fa=nq;
			while(t[p].nxt[c]==q&&p)t[p].nxt[c]=nq,p=t[p].fa;
		}
	}
}

inline void add(int &now,int k,int num)
{
	if(!now)now=++tot;
	if(!~k)return;
	if(num&(1<<k))add(ch[now][1],k-1,num);
	else add(ch[now][0],k-1,num);
}
inline int query(int now,int k,int num)
{
	if(!now||!~k)return 0;
	int d=(num>>k)&1;
	if(ch[now][!d])return (1<<k)+query(ch[now][!d],k-1,num);
	return query(ch[now][d],k-1,num);
}
inline int merge(int x,int y)
{
	if(!x||!y)return x|y;
	ch[x][0]=merge(ch[x][0],ch[y][0]);
	ch[x][1]=merge(ch[x][1],ch[y][1]);
	return x;
}

int main()
{
	n=getint();
	scanf("%s",s+1);
	for(int i=1;i<=n;i++)w[i]=getint();
	for(int i=n;i;i--)
	{
		insert(s[i]-97);
		V[lst].push_back(w[i]),add(rt[lst],16,w[i]);
	}
	for(int i=1;i<=cnt;i++)b[t[i].len]++;![](https://img2020.cnblogs.com/blog/1891964/202006/1891964-20200609165630220-842229293.png)


	for(int i=1;i<=t[lst].len;i++)b[i]+=b[i-1];
	for(int i=1;i<=cnt;i++)pos[b[t[i].len]--]=i;
	for(int i=cnt;i>1;i--)
	{
		int u=pos[i];
		if(V[t[u].fa].size()<V[u].size())swap(rt[t[u].fa],rt[u]),swap(V[t[u].fa],V[u]);
		int mx=0;
		for(int j=0;j<V[u].size();j++)
		{
			V[t[u].fa].push_back(V[u][j]);
			mx=max(mx,query(rt[t[u].fa],16,V[u][j]));
		}
		rt[t[u].fa]=merge(rt[t[u].fa],rt[u]);
		ans=max(ans,t[t[u].fa].len+mx);
	}
	printf("%d
",ans);
}

原文地址:https://www.cnblogs.com/Darknesses/p/13073729.html