【十二省联考2019】【LOJ3049】【洛谷P5248】【BZOJ5496】字符串问题（SAM优化建图）（拓扑排序）

LOJ传送门

洛谷传送门

BZOJ传送门

---

解析：

瞎写一通成功拿下BZOJ倒数第一（在线丢人 ）

常数实在是太大了没办法。。。

首先我们很显然发现要求的是DAG最长链，这是一个图论问题。只要把图建出来一切都好说。

如果建出图有环肯定就咕咕咕了。

显然不能暴力建图，考虑怎么优化前缀那部分的连边。

建立反串后缀自动机，那么前缀关系就转化成了后缀。

树上倍增找到每个串的对应节点，用一个vector存下来。

将SAM上每个节点拆成入点和出点，中间放的是这个节点存的串的信息。所有父亲的出点向儿子的入点连边。显然父亲表示串全部都是儿子串的后缀。

现在考虑怎么处理同一个点的连边。

直接按照长度为第一关键字，是否为$a$类串为第二关键字（$a$串如果和某一个$b$相同，显然应该放在$b$的出来方向）对所有串排序。

然后把串给连一连就行了。这样建立的图显然是我们需要求的图

---

代码：

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define cs const

namespace IO{
	inline int getint(){
		re char c;
		while(!isdigit(c=gc()));re int num=c^48;
		while(isdigit(c=gc()))num=((num+(num<<2))<<1)+(c^48);
		return num;
	}
}
using namespace IO;

using std::cout;
using std::cerr;

typedef std::pair<int,int> pii;
#define mp std::make_pair
#define fi first
#define se second

cs int N=8e5+5;

int n,na,nb;

namespace Graph{
	cs int N=::N<<1;
	std::vector<int> G[N];
	int deg[N],n;
	
	inline void addedge(int u,int v){
		G[u].push_back(v);++deg[v];
		n=std::max(n,std::max(u,v));
	}
	
	ll val[N],len[N];
	std::queue<int> q;
	inline void solve(){
		for(int i=1;i<=n;++i)if(!deg[i])q.push(i);
		int tot=0;
		while(!q.empty()){
			int u=q.front();q.pop();
			len[u]+=val[u];++tot;
			for(int re e=G[u].size()-1;~e;--e){
				int v=G[u][e];
				len[v]=std::max(len[v],len[u]);
				if(!--deg[v])q.push(v);
			}
		}
		if(tot!=n)cout<<"-1\n";
		else cout<<*std::max_element(len+1,len+n+1)<<"\n";
	}
	
	inline void clear(){
		while(n){
			G[n].clear();
			val[n]=len[n]=deg[n]=0;
			n--;
		}
	}
}

namespace SAM{
	cs int N=::N<<1;
	cs int B=20;
	std::vector<pii> nd[N];
	int len[N],fa[N],son[N][26],f[N][B+1];
	int pos[N],id[N];
	int last=1,now=1,tot;
	
	inline void push_back(int c,int id){
		c-='a';
		int cur=++now,p=last;
		len[cur]=len[last]+1;pos[id]=cur;
		for(;p&&!son[p][c];p=fa[p])son[p][c]=cur;
		if(!p)fa[cur]=1;
		else if(len[son[p][c]]==len[p]+1)fa[cur]=son[p][c];
		else {
			int clone=++now,q=son[p][c];
			len[clone]=len[p]+1;
			memcpy(son[clone],son[q],sizeof son[q]);
			fa[clone]=fa[q];
			fa[q]=fa[cur]=clone;
			for(;p&&son[p][c]==q;p=fa[p])son[p][c]=clone;
		}
		last=cur;
	}
	
	inline void init(){
		for(int re i=2;i<=now;++i)f[i][0]=fa[i],Graph::addedge(fa[i]+now,i);
		for(int re i=1;i<=B;++i)
		for(int re j=1;j<=now;++j)f[j][i]=f[f[j][i-1]][i-1];
	}
	
	inline void locate(int l,int r,int id){
		int u=pos[l];
		for(int re i=B;~i;--i)if(len[f[u][i]]>=r-l+1)u=f[u][i];
		nd[u].push_back(mp(r-l+1,id));
	}
	
	inline bool cmp(cs pii &a,cs pii &b){
		return a.fi!=b.fi?a.fi<b.fi:a.se>b.se;
	}
	
	inline void build(){
		for(int re i=1;i<=now;++i){
			std::vector<pii> &vec=nd[i];
			if(!vec.size())Graph::addedge(i,i+now);
			else {
				sort(vec.begin(),vec.end(),cmp);
				for(int re j=0;j<(int)vec.size();++j){
					int x=vec[j].se;
					id[x]=x+now*2,x=id[x];
					if(j==0)Graph::addedge(i,x);
					if(j==(int)vec.size()-1)Graph::addedge(x,i+now);
					else Graph::addedge(x,vec[j+1].se+now*2);
					if(x<=now*2+na){
						++tot;
						int u=now*2+na+nb+tot;
						Graph::addedge(x,u);
						Graph::val[u]=vec[j].fi,id[x-now*2]=u;
					}
				}
			}
		}
	}
	
	inline void clear(){
		while(now){
			nd[now].clear();
			len[now]=fa[now]=0;
			memset(son[now],0,sizeof son[now]);
			memset(f[now],0,sizeof f[now]);
			--now;
		}
		last=now=1;tot=0;
	}
}

char s[N];
inline void solve(){
	std::cin>>(s+1);n=strlen(s+1);
	for(int re i=n;i;--i)SAM::push_back(s[i],i);
	SAM::init();
	na=getint();
	for(int re i=1,l,r;i<=na;++i){
		l=getint(),r=getint();
		SAM::locate(l,r,i);
	}
	nb=getint();
	for(int re i=1,l,r;i<=nb;++i){
		l=getint(),r=getint();
		SAM::locate(l,r,na+i);
	}
	SAM::build();
	int m=getint();
	while(m--){
		int x=getint(),y=getint();
		Graph::addedge(SAM::id[x],SAM::id[y+na]);
	}
	Graph::solve();
}

inline void init(){
	memset(SAM::id+1,0,sizeof(int)*(na+nb));
	memset(SAM::pos+1,0,sizeof(int)*n);
	SAM::clear();
	Graph::clear();
}

signed main(){
	int T=getint();
	while(T--)init(),solve();
	return 0;
}