51nod除了在做OJ之外,还开展了很多副业。婚姻介绍所就是其中之一。
对于一个客户,我们可以使用一个字符串来描述该客户的特质。
假设现在我们有两个客户A和B。
A的特质字符串为:abcdefg
B的特质字符串为:abcxyz
则A和B的匹配度f(A, B)为A和B的最长公共前缀的长度,即len(‘abc’) = 3
由于最近51nod经费紧张,所以夹克大老爷设计了一种压缩算法以节约内存。
所有用户的特质字符串都被存储在了一个长为n的字符串S中。(n <= 1000)用户的特质使用一个整数p表示,表示该用户的特质字符串为S[p…n - 1]。
现给定字符串S,与q次查询<ai, bi>(ai, bi分别为合法的用户特质整数)。请输出q次查询分别对应的客户匹配度。
现给定字符串长度n,与字符串S。接下来是整数q,代表接下来有q次查询。
下面q行有两个整数ai, bi。代表查询特质为ai与bi的用户的匹配度。
1 <= n <= 1000
1 <= q <= 10^6
输入数据全部合法。
每一行输出一个用户匹配度整数。
题意:求给定x位置的后缀与y位置的后缀的最长公共前缀
思路:后缀x和后缀y的最长公共前缀
min(height[rak[x]+1], height[rak[x] + 2]…height[rak[y]])
使用RMQ预处理出height数组,然后O(1)时间复杂度查询最长公共前缀
#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e6 + 5;
char s[maxn];
struct node {
int x, y, id;
}a[maxn], b[maxn];
int n, m, p;
int sa[maxn], rak[maxn], tx[maxn];
int lg[maxn];
void rsort() {
for (int i = 1; i <= m; i++) {
tx[i] = 0;
}
for (int i = 1; i <= n; i++) {
tx[a[i].y]++;
}
for (int i = 1; i <= m; i++) {
tx[i] += tx[i - 1];
}
for (int i = 1; i <= n; i++) {
b[tx[a[i].y]--] = a[i];
}
for (int i = 1; i <= m; i++) {
tx[i] = 0;
}
for (int i = 1; i <= n; i++) {
tx[b[i].x]++;
}
for (int i = 1; i <= m; i++) {
tx[i] += tx[i - 1];
}
for (int i = n; i >= 1; i--) {
a[tx[b[i].x]--] = b[i];
}
}
void solve() {
rsort();
p = 0;
for (int i = 1; i <= n; i++) {
if (a[i].x == a[i - 1].x && a[i].y == a[i - 1].y) {
rak[a[i].id] = p;
} else {
rak[a[i].id] = ++p;
}
}
for (int i = 1; i <= n; i++) {
a[i].x = rak[i];
a[i].id = sa[rak[i]] = i;
a[i].y = 0;
}
m = p;
}
void ssort() {
m = 127;
for (int i = 1; i <= n; i++) {
a[i].x = a[i].y = s[i];
a[i].id = i;
}
solve();
for (int j = 1; j <= n; j <<= 1) {
for (int i = 1; i + j <= n; i++) {
a[i].y = a[i + j].x;
}
solve();
if (p == n) {
break;
}
}
}
int height[maxn];
inline void get_Height(){//求Height数组
int k=0;
for(int i=1;i<=n;++i){
if(k)k--;
int j=sa[rak[i]-1];
while(s[i+k]==s[j+k])++k;
height[rak[i]]=k;
}
}
int dp[maxn][20];
void RMQ() {
lg[0] = -1;
for (int i = 1; i <= n; i++) {
dp[i][0] = height[i];
lg[i] = lg[i >> 1] + 1;
}
for (int j = 1; j <= lg[n]; j++) {
for (int i = 1; i + (1 << j) - 1 <= n; i++) {
dp[i][j] = min(dp[i][j - 1], dp[i + (1 << (j - 1))][j - 1]);
}
}
}
inline int read(){
int k=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){k=k*10+ch-'0';ch=getchar();}
return k*f;
}
int main() {
n = read();
scanf("%s", s + 1);
n = strlen(s + 1);
ssort();
get_Height();
m = read();
int x, y, l, r;
RMQ();
while (m--) {
x = read(); y = read();
x++; y++;
if (x == y) {
printf("%d\n", n - x + 1);
} else {
l = rak[x];
r = rak[y];
if (l > r) {
swap(l, r);
}
l++;
int k = lg[r - l + 1];
printf("%d\n", min(dp[l][k], dp[r - (1 << k) + 1][k]));
}
}
return 0;
}