2746: [HEOI2012]旅行问题

2746: [HEOI2012]旅行问题

Time Limit: 30 Sec  Memory Limit: 256 MB
Submit: 718  Solved: 243
[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。 

Source

[Submit][Status][Discuss]



AC自动机。。

把所有串串扔进自动机构造fail tree

每次询问输出对应节点的LCA就是

6666666666666666666666666666666

一开始数组开太大了sad。。。

因为是用后悔边建的树，所以深度要在递推完f数组后才能算(第二次写的时候不幸出错)

#include<iostream>#include<cstdio>#include<cstring>#include<queue>#include<vector>using namespace std;const int maxn = 1E6 + 10;typedef long long LL;const LL mo = 1000000007LL;int n,m,cnt = 0,ch[maxn][26],f[maxn][20],L[maxn];LL ans[maxn];char c[maxn];vector <int> pos[maxn];queue <int> q;int LCA(int a,int b){if (L[a] < L[b]) swap(a,b);int log;for (log = 0; L[a] - (1<<log) >= 0; log++); log--;for (int j = log; j >= 0; j--)if (L[a] - (1<<j) >= L[b])a = f[a][j];if (a == b) return ans[a];for (int j = log; j >= 0; j--)if (f[a][j] != f[b][j]) {a = f[a][j]; b = f[b][j];}return ans[f[a][0]];}int main(){#ifdef YZYfreopen("yzy.txt","r",stdin);#endifcin >> n;for (int i = 1; i <= n; i++) {scanf("%s",c);int len = strlen(c);int x = 0;for (int j = 0; j < len; j++) {int po = c[j] - 'a';if (!ch[x][po]) {ch[x][po] = ++cnt;LL ADD = c[j] - 'a';ans[cnt] = (ans[x]*26LL + ADD)%mo;}x = ch[x][po];pos[i].push_back(x);}}for (int i = 0; i < 26; i++)if (ch[0][i]) {L[ch[0][i]] = 1;f[ch[0][i]][0] = 0;q.push(ch[0][i]);}while (!q.empty()) {int k = q.front(); q.pop();for (int i = 0; i < 26; i++) {int u = ch[k][i];if (!u) {ch[k][i] = ch[f[k][0]][i];continue;}int v = f[k][0];while (v && !ch[v][i]) v = f[v][0];L[u] = L[ch[v][i]] + 1;f[u][0] = ch[v][i];q.push(u);}}for (int j = 1; j < 20; j++)for (int i = 1; i <= cnt; i++)f[i][j] = f[f[i][j-1]][j-1];cin >> m;while (m--) {int a,b,x,y;scanf("%d%d%d%d",&a,&b,&x,&y);printf("%d\n",LCA(pos[a][b-1],pos[x][y-1]));}return 0;}