OKR-A Horrible Poem
发布日期:2021-05-07 06:54:53 浏览次数:21 分类:精选文章

本文共 2619 字,大约阅读时间需要 8 分钟。

O K R − A   H o r r i b l e   P o e m OKR-A\ Horrible\ Poem OKRA Horrible Poem

题目链接:

题目大意

给出一个由小写英文字母组成的字符串 S S S,再给出 q q q个询问,要求回答 S S S某个子串的最短循环节。

如果字符串 B B B是字符串 A A A的循环节,那么 A A A可以由 B B B 重复若干次得到。

输入

第一行一个正整数 n ( n ≤ 500   000 ) n(n≤500\ 000) n(n500 000),表示 S S S的长度。

第二行 n n n个小写英文字母,表示字符串 S S S
第三行一个正整数 q ( q ≤ 2   000   000 ) q(q≤2\ 000\ 000) q(q2 000 000),表示询问次数。
下面 q q q行每行两个正整数 a , b ( 1 ≤ a ≤ b ≤ n ) a,b(1≤a≤b≤n) a,b(1abn),表示询问字符串 S [ a … b ] S[a…b] S[ab]的最短循环节长度。

输出

给出一个由小写英文字母组成的字符串 S S S,再给出 q q q个询问,要求回答 S S S某个子串的最短循环节。

如果字符串 B B B是字符串 A A A的循环节,那么 A A A可以由 B B B重复若干次得到。

样例输入

8aaabcabc31 33 84 8

样例输出

135

思路

这道题用 h a s h hash hash来做。

因为要求字符串的最短循环节,就要判断字符串是否相同,自然就会用到 h a s h hash hash

既然知道了怎么判断两个字符串是否相同,就要判断某个小字符串是否是某个字符串的循环节。当然,我们不能直接一个一个地枚举那个字符串里面的小字符串(因为会超时),那怎么办呢?

有一种方法:设 a [ i ] [ j ] a[i][j] a[i][j] a a a字符串从第 i i i为到第 j j j位,如果我们要判断 b [ k ] [ l ] b[k][l] b[k][l]这个小字符串(长度为 l e n len len)是它的循环节的话,我们只要判断 a [ i ] [ j − l e n ] = = a [ i + l e n ] [ j ] a[i][j−len]==a[i+len][j] a[i][jlen]==a[i+len][j]。原因的话自己模拟看看就知道了。

在枚举循环节的长度的时候,我们不一定要枚举全部的长度,我们只要在循环节循环次数是字符串的长度的约数的情况下才需要枚举这个长度。

代码

#include
#define rr register#define min(x, y) (x) < (y) ? (x) : (y) using namespace std;int n, m, x, y, p = 131, P[500001], hash[500001], b[500001][27], ans;char a[500001];int geth(int l, int r) { //得到字符串的hash值 return hash[r] - hash[l - 1] * P[r - l + 1];}void find (int length) { //如果是判断出是循环节的话就判断是否比当前最优解还要优 if (geth(x + length, y) == geth(x, y - length)) ans = min(ans, length);}int gcd(int l, int r) { //求gcd if (r == 0) return l; return gcd(r, l % r);}int read() { //快读 int an = 0; char c = getchar(); while (c > '9' || c < '0') c = getchar(); while (c <= '9' && c >= '0') { an = an * 10 + c - 48; c = getchar(); } return an;}int main() { n = read();//读入 scanf("%s", a + 1);//读入 P[0] = 1;//初始化 for (rr int i = 1; i <= n; i++) P[i] = P[i - 1] * p;//预处理hash的进制 for (rr int i = 1; i <= n; i++) hash[i] = hash[i - 1] * p + a[i] - 'a' + 1; //转换hash并记录前缀和 for (rr int i = 1; i <= n; i++) { for (int j = 1; j <= 26; j++) b[i][j] = b[i - 1][j];//前缀和出单词出现的位置与频率 b[i][a[i] - 'a' + 1]++;//记录单词的出现 } scanf("%d", &m);//读入 for (rr int i = 1; i <= m; i++) { x = read();//读入 y = read();//读入 int len = y - x + 1;//求出长度 ans = 2147483647;//初始化 for (int i = 1; i <= 26; i++) len = gcd(len, b[y][i] - b[x - 1][i]);//对小序列中每个字母出现次数取gcd for (int i = 1; i * i <= len; i++) { //枚举循环次数 if (len % i) continue;//循环次数须为序列长度的约数 find((y - x + 1) / i);//计算循环次数为i次的结果,更新答案 find((y - x + 1) / (len / i));//计算循环次数为len/i次的结果,更新答案 } printf("%d\n", ans);//输出 } return 0;}
上一篇:变幻数
下一篇:Mathematics

发表评论

最新留言

路过,博主的博客真漂亮。。
[***.116.15.85]2025年04月12日 06时33分27秒