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
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 YZY
freopen("yzy.txt","r",stdin);
#endif
cin >> 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;
}