BZOJ 2795: [Poi2012]A Horrible Poem( hash )

...字符串hash.

假如长度x是一个循环节, 那么对于任意n(x | n)也是一个循环节.

设当前询问区间[l, r]长度为len = ∏piai, 最终答案ans = ∏piai' ,我们只需枚举len的质因数来确定ai'即可

#include<bits/stdc++.h>

using namespace std;

typedef unsigned long long ull;

const int maxn = 500009;
const ull P = 1000173169;

char S[maxn];
int check[maxn], prime[maxn], N = 0, n;
ull H[maxn], K[maxn];

void init() {
memset(check, 0, sizeof check);
for(int i = 2; i <= n; i++) {
if(!check[i])
check[i] = prime[N++] = i;
for(int j = 0; j < N && i * prime[j] <= n; j++) {
check[i * prime[j]] = prime[j];
if(i % prime[j] == 0) break;
}
}
H
= 0; K[0] = 1;
for(int i = n - 1; ~i; i--)
H[i] = H[i + 1] * P + S[i];
for(int i = 1; i <= n; i++)
K[i] = K[i - 1] * P;
}

inline bool jud(int a, int b, int L) {
return H[a] - H[a + L] * K[L] == H[b] - H[b + L] * K[L];
}

void Read() {
cin >> n;
char* t = S;
for(; !islower(*t); *t = getchar());
for(int i = 1; i < n; i++)
*++t = getchar();
}

int main() {

Read();
init();

int q; scanf("%d", &q);
while(q--) {
int l, r; scanf("%d%d", &l, &r); l--; r--;
int ans = r - l + 1, t = ans;
while(t > 1) {
int x = check[t];
while(ans % x == 0 && jud(l, l + ans / x, r - l  + 1 - ans / x)) ans /= x;
for(; t % x == 0; t /= x);
}
printf("%d\n", ans);
}

return 0;
}

　　

2795: [Poi2012]A Horrible Poem

Time Limit: 50 Sec Memory Limit: 128 MB
Submit: 202 Solved: 115
[Submit][Status][Discuss]

Description

给出一个由小写英文字母组成的字符串S，再给出q个询问，要求回答S某个子串的最短循环节。
如果字符串B是字符串A的循环节，那么A可以由B重复若干次得到。

Input

第一行一个正整数n (n<=500,000)，表示S的长度。
第二行n个小写英文字母，表示字符串S。
第三行一个正整数q (q<=2,000,000)，表示询问个数。
下面q行每行两个正整数a,b (1<=a<=b<=n)，表示询问字符串S[a..b]的最短循环节长度。

Output

依次输出q行正整数，第i行的正整数对应第i个询问的答案。

Sample Input

8

aaabcabc

3

1 3

3 8

4 8

Sample Output

1

3

5

HINT

Source

鸣谢 jiangzoi&oimaster