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


你可能感兴趣的:(2746: [HEOI2012]旅行问题)