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