并不对劲的loj3049:p5284:[十二省联考]字符串问题

题目大意

给出字符串(S(|S|leq2 imes10^5))，
(na(naleq2 imes 10^5))个区间([l_i,r_i])表示(S_{l_i},S_{l_i+1},...,S_{r_i})组成的这个(S)的子串是第(i)个A类串
(nb(nbleq2 imes 10^5))个区间([l_i,r_i])表示(S_{l_i},S_{l_i+1},...,S_{r_i})组成的这个(S)的子串是第(i)个B类串
(m(mleq2 imes10^5))个限制条件((x,y))表示第(x)个A类串后面可以接以第(y)个B类串为前缀的A类串
现在想要生成字符串(T)，满足(T)由若干个A类串拼接而成，而且相邻的A类串之间满足限制条件
求(T)的最长长度，或判断(T)可以无限长

题解

发现将每个A类串看成一个点，点权为它的长度，将每个A类串向它后面能接的A类串连边
问题就变成了找环或求点权之和最大的链
暴力连边肯定不行，考虑后缀树优化建图
连出的边的形态大概是：(A_i->B_i->以B_i为前缀的A)
那么考虑将后缀树上的每个点(x)拆成两个点(x_A,x_B)
要连的边就是：(x_B->x_A,x_{A_i}->x_{B_i})和从(x_B)连向它后缀树的儿子的(x_B)的边
还有个问题：已知区间([l,r])，在(S)的后缀树中找到(S[l:r])对应的点
后缀树中每个叶子结点是(S)的后缀，父亲是儿子的前缀
那么对于(forall iin[1,|S|])可以在建后缀树时记(S[i:|S|])的位置
找(S[l:r])对应的位置就变成了先找(S[l:|S|])的位置，再在这个点的祖先中找到(dis)值不小于(r-l+1)的最高点
这一步可以预处理每个点第(2^k)个祖先+倍增
还有一点要注意的是，后缀树中的每个点可能代表了多个串，所以还要再把每个点内部拆一次
并没有想象中那么难写

代码

#include<algorithm>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<ctime>
#include<iomanip>
#include<iostream>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<vector>
#define rep(i,x,y) for(register int i=(x);i<=(y);++i)
#define dwn(i,x,y) for(register int i=(x);i>=(y);--i)
#define view(u,k) for(int k=fir[u];~k;k=nxt[k])
#define maxn 200010
#define maxm (maxn<<1)
#define LL long long
#define pb push_back
using namespace std;
int read()
{
	int x=0,f=1;char ch=getchar();
	while(!isdigit(ch)&&ch!='-')ch=getchar();
	if(ch=='-')f=-1,ch=getchar();
	while(isdigit(ch))x=(x<<1)+(x<<3)+ch-'0',ch=getchar();
	return x*f;
}
void write(LL x)
{
	if(x==0){putchar('0'),putchar('
');return;}
	int f=0;char ch[20];
	if(x<0)putchar('-'),x=-x;
	while(x)ch[++f]=x%10+'0',x/=10;
	while(f)putchar(ch[f--]);
	putchar('
');
	return;
}
char s[maxn];
int dfn[maxm<<2],in[maxm<<2],yes,n,ch[maxm][26],dis[maxm],fa[maxm],posa[maxn],posb[maxn],pos[maxn],rt,cntnd,lst,c[maxn],ord[maxm];
int m,T,na,nb,fir[maxm<<2],nxt[maxm<<3],v[maxm<<3],cnte,anc[maxm][20],low[maxm<<2],ins[maxm<<2],tim,stk[maxm<<2],tp,tot;
int Q[maxm<<2],hd,tl,ltha[maxn],lthb[maxn];
LL maxd[maxm<<2],ans;
int w[maxm<<2];
vector <int >val[maxm],id[maxm];
#define gx(C) (C-'a')
void ade(int u1,int v1){v[cnte]=v1,nxt[cnte]=fir[u1],fir[u1]=cnte++;}
void ext(int i)
{
	int p=lst,val=gx(s[i]),np=++cntnd;dis[np]=i,lst=np,pos[n-i+1]=np;
	for(;p&&!ch[p][val];p=fa[p])ch[p][val]=np;
	if(!p)fa[np]=rt;
	else
	{
		int q=ch[p][val];
		if(dis[p]+1==dis[q])fa[np]=q;
		else
		{
			int nq=++cntnd;
			dis[nq]=dis[p]+1;
			memcpy(ch[nq],ch[q],sizeof(ch[q]));
			fa[nq]=fa[q],fa[q]=fa[np]=nq;
			for(;p&&ch[p][val]==q;p=fa[p])ch[p][val]=nq;
		} 
	}
}
void reset()
{
	rep(i,1,n)c[i]=0;
	rep(i,1,cntnd)rep(j,0,25)ch[i][j]=0;
	rep(i,1,cntnd)rep(j,0,19)anc[i][j]=0;
	rep(i,1,cntnd){fa[i]=0;}
	dwn(i,(tot<<1),1)dfn[i]=in[i]=maxd[i]=0,fir[i]=-1,w[i]=0;
	rt=lst=cntnd=hd=1;tl=cnte=yes=ans=tim=tp=tot=0; 
}
int getp(int l,int r)
{
	int len=r-l+1,u=pos[l];
	dwn(i,19,0)if(anc[u][i]&&dis[anc[u][i]]>=len)u=anc[u][i];
	return u;
} 
void tar(int u)
{
	dfn[u]=low[u]=++tim;
	ins[u]=1,stk[++tp]=u;
	view(u,k)
	{
		if(!dfn[v[k]])tar(v[k]),low[u]=min(low[u],low[v[k]]);
		else if(ins[v[k]])low[u]=min(low[u],dfn[v[k]]);
	}
	if(low[u]==dfn[u])
	{
		int num=0;
		while(1)
		{
			num++;
			ins[stk[tp]]=0;
			if(stk[tp--]==u)break;
		}
		if(num>1)yes=1;
	}
}
int getpos(int u,int len)
{
	int l=0,r=val[u].size()-1,res=r;
	while(l<=r)
	{
		int mid=(l+r)>>1;
		if(val[u][mid]>=len)res=min(res,mid),r=mid-1;
		else l=mid+1;
	}
	return id[u][res];
}
int main()
{
	T=read();
	memset(fir,-1,sizeof(fir));rt=lst=++cntnd;hd=1,tl=0;
	while(T--)
	{
		scanf("%s",s+1);
		n=strlen(s+1);
		reverse(s+1,s+n+1);
		rep(i,1,n)ext(i);
		rep(i,1,cntnd)c[dis[i]]++;
		rep(i,1,n)c[i]+=c[i-1];
		rep(i,1,cntnd)ord[c[dis[i]]--]=i;
		rep(j,1,cntnd)
		{
			int i=ord[j];
			anc[i][0]=fa[i];
			rep(k,1,19)anc[i][k]=anc[anc[i][k-1]][k-1];
		}
		na=read();
		rep(i,1,na){int l=read(),r=read();posa[i]=getp(l,r),val[posa[i]].pb(r-l+1),ltha[i]=r-l+1;}
		nb=read();
		rep(i,1,nb){int l=read(),r=read();posb[i]=getp(l,r),lthb[i]=r-l+1;}
		rep(i,1,cntnd)sort(val[i].begin(),val[i].end());
		rep(j,1,cntnd)
		{
			int lim=val[ord[j]].size(),Lst=0;if(fa[ord[j]]){Lst=id[fa[ord[j]]][id[fa[ord[j]]].size()-1];} 
			rep(k,0,lim-1){++tot,w[tot]=val[ord[j]][k];if(fa[ord[j]])ade(Lst,tot),++in[tot];Lst=tot,id[ord[j]].pb(tot);}
			if(!lim||val[ord[j]][lim-1]!=dis[ord[j]]){++tot;if(fa[ord[j]])ade(Lst,tot),++in[tot];id[ord[j]].pb(tot);val[ord[j]].pb(dis[ord[j]]);}
		}
		m=read();
		rep(i,1,tot)ade(i,i+tot),++in[i+tot],w[i+tot]=w[i],w[i]=0;
		rep(i,1,m)
		{
			int x=read(),y=read();
			x=getpos(posa[x],ltha[x]),y=getpos(posb[y],lthb[y]);
			ade(x+tot,y),in[y]++;
		}
		{rep(i,1,cntnd){val[i].clear(),val[i].shrink_to_fit();id[i].clear(),id[i].shrink_to_fit();}}
		dwn(i,(tot<<1),1)if(!dfn[i])tar(i);
		if(yes)puts("-1");
		else
		{
			dwn(i,(tot<<1),1)if(!in[i])Q[++tl]=i;
			while(hd<=tl)
			{
				int u=Q[hd++];maxd[u]+=w[u],ans=max(ans,maxd[u]);
				view(u,k)
				{
					maxd[v[k]]=max(maxd[v[k]],maxd[u]),in[v[k]]--;
					if(!in[v[k]])Q[++tl]=v[k];
				}
			}
			write(ans);
		}
		reset();
	}
	return 0;
}

一些感想

据说猎人公会会长是云猎人。。。我佛了
sb猎人公会nmsl