LOJ6041 「雅礼集训 2017 Day7」事情的相似度

题目传送门

分析:
前缀的公共后缀,那直接SAM安排上
两个前缀的公共后缀为两个串对应Parent树上节点的LCA
于是\([l,r]\)之间的公共后缀最长对应这些点两两在Parent树上LCA的最大的len

Set启发式合并做法:
建出了Parent树,我们考虑每一个点\(u\),合并它的子树集合时,两个集合中各取一个点的LCA一定是\(u\)
我们可以找到若干对点ID相近且LCA为\(u\),他们会对答案造成贡献,而每次合并点集,产生的对数为较小的集合的大小
使用启发式合并,空间复杂度是\(O(nlogn)\),时间复杂度\(O(nlog^{2}n)\)
点对预处理出来之后,可以离线+树状数组查询答案
代码可以看yyb这篇博客,很好理解

LCT+树状数组做法:
我们将询问按\(r\)离线,每次按顺序加入一个点\(i\)时,我们将这个点到根的路径全部染成颜色\(i\),而某个点\(u\)原来的颜色\(j\)被覆盖,那么\(u\)就是\(i,j\)的LCA
每个点的颜色都是最近被覆盖的,得到的点对与上面Set做法得到的点对是一样的
到根路径染色实际上是LCT的access过程,于是直接用LCT维护,过程中即可更新树状数组求答案
时间复杂度\(O(nlog^{2}n)\)

代码是第二种做法的:

#include
#include
#include
#include
#include
#include
#include

#define maxn 400005
#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,m;
int Id[maxn];
struct node{
	int fa,nxt[2],len;
}t[maxn];
int lst,tot;
int ch[maxn][2],val[maxn],tag[maxn],fa[maxn];
vectorq[maxn];
int p[maxn];
int stk[maxn],tp;
int mx[maxn],ans[maxn];
char s[maxn];

inline void insert(int c)
{
	int p=lst,np=lst=++tot;
	t[np].len=t[p].len+1;
	while(p&&!t[p].nxt[c])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=++tot;
			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(p&&t[p].nxt[c]==q)t[p].nxt[c]=nq,p=t[p].fa;
		}
	}
}

inline void update(int x,int num)
{for(int i=n-x+1;i<=n;i+=i&(-i))mx[i]=max(mx[i],num);}
inline int query(int x)
{int num=0;for(int i=n-x+1;i;i-=i&(-i))num=max(num,mx[i]);return num;}
inline bool isroot(int x)
{return ch[fa[x]][0]!=x&&ch[fa[x]][1]!=x;}
inline void change(int x,int num){val[x]=tag[x]=num;}
inline void pushdown(int x)
{
	if(!tag[x])return;
	if(ch[x][0])change(ch[x][0],tag[x]);
	if(ch[x][1])change(ch[x][1],tag[x]);
	tag[x]=0;
}
inline void rotate(int x)
{
	int y=fa[x],z=fa[y];
	int l=(ch[y][1]==x),r=l^1;
	if(!isroot(y))
	{
		if(ch[z][0]==y)ch[z][0]=x;else ch[z][1]=x;
	}
	fa[x]=z,fa[y]=x,ch[y][l]=ch[x][r],fa[ch[x][r]]=y,ch[x][r]=y;
}
inline void splay(int x)
{
	stk[++tp]=x;for(int i=x;!isroot(i);i=fa[i])stk[++tp]=fa[i];
	while(tp)pushdown(stk[tp--]);
	while(!isroot(x))
	{
		int y=fa[x],z=fa[y];
		if(!isroot(y))
		{
			if((ch[z][0]==y)^(ch[y][0]==x))rotate(x);
			else rotate(y);
		}
		rotate(x);
	}
}
inline void access(int x,int num)
{
	int tmp=0;
	for(;x;ch[x][1]=tmp,tmp=x,x=fa[x])
		splay(x),update(val[x],t[x].len);
	change(tmp,num);
}

int main()
{
	n=getint(),m=getint();
	scanf("%s",s+1);
	for(int i=1;i<=m;i++)p[i]=getint(),q[getint()].push_back(i);
	lst=tot=1;
	for(int i=1;i<=n;i++)insert(s[i]-'0'),Id[i]=lst;
	for(int i=1;i<=tot;i++)fa[i]=t[i].fa;
	for(int i=1;i<=n;i++)
	{
		access(Id[i],i);
		int sz=q[i].size();
		for(int j=0;j

你可能感兴趣的:(LOJ6041 「雅礼集训 2017 Day7」事情的相似度)