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 20using 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;}