bzoj2746 旅行问题

2746: [HEOI2012]旅行问题

Time Limit: 30 Sec  Memory Limit: 256 MB
[Submit][Status][Discuss]

Description

yz是Z国的领导人，他规定每个地区的名字只能为26个小写拉丁字母的一个。由于地 区数有可能超过26个，便产生了一个问题，如何辨别名字相同的地区？于是yz规定，一个 地区的描述必须包含它的所有上级，且上级按次序排列。于是，一个地区的描述是一个字符 串。比如说，一个地区的名字为c，它的上级为b，b的上级为a，a没有上级，那么这个地 区就描述为abc。显然，这个描述同时包含了c的上级b和b的上级a的描述，分别为ab和a。 值得注意的是，每个地区最多有一个上级，同一上级的地区之间名字不同，没有上级的 地区之间名字不同。现在，yz对外公布了n个地区的描述，这些描述中包含了Z国所有地区的描述，并让 你处理来访者的旅行问题。现有m对人访问这个国家，对于每对人，第一个人喜欢第i个描述中的第j个地区，设 这个地区描述为s1，第二个人喜欢第k个描述中的第l个地区，设这个地区描述为s2。他们为了统一行程，决定访问描述为s的地区（显然他们只关心地区的名字，并非是地区本身）， 设s的长度为t，s需要满足以下条件： 
1：t<=j, t<=l; 
1：s[1..t] = s1[j-t+1 … j], s[1..t] = s2[l-t+1 … l]；（即s为s1中1到k位 与s2中1到l位的公共后缀） 
2：t最大化。 
为了不使输出过大，你只需把这个字符串按照如下生成的26进制数转成10进制后mod 1000000007后输出: 
a->0 
b->1 
. 
. 
. 
z->25 
比如地区cab被编码成2 *    26? + 0 * 26? + 1 * 26? = 1353。 

Input

第一行给定一个整数n 
第2…n+1行:每i+1行给定一个字符串a[i],表示第i个描述。 
接下来一行一个整数m 
接下来m行:每行给定四个整数i,j,k,l，字母含义与题目描述一致。 

Output

共m行，每行一个整数，表示答案字符串的编码。 

Sample Input

2
aabb babb
2
1 3 2 3
1 4 2 4

Sample Output

1
1
【样例说明】
询问1中的公共后缀有ab和b，但是没有ab这个地区，只有b地区，所以只能选择b这个 地区；
询问2中的公共后缀有abb、bb和b，但是没有abb和bb这两个地区，只有b地区，所以 只能选择b这个地区。

HINT

【数据范围】

 设这个国家地区总数数为tot（注意：输入的字符串总长度可能超过tot！） 对于30%的数据，满足tot，m，n<=100； 

对于50%的数据，满足tot，m，n<=1000； 

对于80%的数据，满足tot，m，n<=100000； 

对于100%的数据，满足tot，m，n<=1000000； 

保证输入文件不超过20MB。 

Analysis

　　呃。。。AC自动机恶心题

　　求的那个串就是在trie树上做LCA（我用的倍增），最恶心的在于卡空间

Code

//AC自动机
#include <cstdio>
#include <algorithm>
#define maxn 1000010
#define maxk 20
using namespace std;
int f[maxn][21], trie[maxn][26], fail[maxn], w[maxn], tot=1, N, M, deep[maxn],
	len[maxn], q[maxn], acaf[maxn], ans[maxn];
char s[5000000];
struct quiry{int a, b, num, pos;}qu[maxn<<1], *pt;
bool cmp1(quiry q1, quiry q2){return q1.a == q2.a ? q1.b < q2.b : q1.a < q2.a;}
bool cmp2(quiry q1, quiry q2){return q1.num < q2.num;}
void insert(char *&s, int num)
{
	int pos, to;
	for(pos=1;*s;s++)
	{
		to = trie[pos][*s-97] ? trie[pos][*s-97] : trie[pos][*s-97]=++tot;
		w[to]=*s;
		len[to]=len[pos]+1;
		acaf[to]=pos;
		ans[to]=((long long)(*s-97)+(long long)ans[pos]*26)%1000000007;
		pos=to;
		while(pt->a==num and pt->b==len[pos])pt->pos=pos,pt++;
	}
}
void acabuild()
{
	int u, f, i, l=1 ,r=1;
	q[r++]=1;
	deep[1]=1;
	for(i=2;i<=tot;i++)fail[i]=1;
	while(l<r)
	{
		u=q[l++];
		for(i=0;i<26;i++)
		{
			if(trie[u][i])
			{
				for(f=fail[u];f and !trie[f][i];f=fail[f]);
				if(f)fail[trie[u][i]]=trie[f][i];
				else fail[trie[u][i]]=1;
				deep[trie[u][i]]=deep[fail[trie[u][i]]]+1;
				q[r++]=trie[u][i];
			}
		}
	}
}
void calcf()
{
	int i, k;
	for(i=2;i<=tot;i++)f[i][0]=fail[i];
	for(k=1;k<=maxk;k++)
		for(i=2;i<=tot;i++)
			f[i][k]=f[f[i][k-1]][k-1];
}
int lca(int x, int y)
{
	int k, t;
	if(deep[x]>deep[y])t=x,x=y,y=t;
	for(k=maxk;deep[y]>deep[x];k--)
		if(deep[f[y][k]]>=deep[x])y=f[y][k];
	if(x==y)return x;
	for(k=maxk;k>=0;k--)
		if(f[x][k]!=f[y][k])x=f[x][k],y=f[y][k];
	return f[x][0];
}
int solve(int num)
{
	int x, y;
	x=qu[num].pos;
	y=qu[M+num].pos;
	printf("%d\n",ans[lca(x,y)]);
}
void init()
{
	int i;
	char *p;
	scanf("%d",&N);
	for(i=1,p=s-1;i<=N;i++)for(scanf("%s",++p);*p;p++);
	scanf("%d",&M);
	for(i=1;i<=M;i++)
		scanf("%d%d%d%d",&qu[i].a,&qu[i].b,&qu[M+i].a,&qu[M+i].b),qu[i].num=i, qu[M+i].num=M+i;
	sort(qu+1,qu+M+M+1,cmp1);
	pt=qu+1;
	for(p=s,i=1;i<=N;i++,p++)insert(p,i);
}
int main()
{
	init();
	acabuild();
	calcf();
	sort(qu+1,qu+M+M+1,cmp2);
	for(int i=1;i<=M;i++)solve(i);
	return 0;
}