
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 OKR−A 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(n≤500 000),表示 S S S的长度。
第二行 n n n个小写英文字母,表示字符串 S S S。 第三行一个正整数 q ( q ≤ 2 000 000 ) q(q≤2\ 000\ 000) q(q≤2 000 000),表示询问次数。 下面 q q q行每行两个正整数 a , b ( 1 ≤ a ≤ b ≤ n ) a,b(1≤a≤b≤n) a,b(1≤a≤b≤n),表示询问字符串 S [ a … b ] S[a…b] S[a…b]的最短循环节长度。输出
给出一个由小写英文字母组成的字符串 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][j−len]==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;}
发表评论
最新留言
路过,博主的博客真漂亮。。
[***.116.15.85]2025年04月12日 06时33分27秒
关于作者

喝酒易醉,品茶养心,人生如梦,品茶悟道,何以解忧?唯有杜康!
-- 愿君每日到此一游!
推荐文章
zmq的send
2021-05-09
阿里钉钉面试题
2021-05-09
C++中找资源或者函数的方法
2021-05-09
delete对象时会自动调用类的析构函数
2021-05-09
SQL联表的方式(逗号, Left Join, Right Join)
2021-05-09
字符串初始化时的注意点
2021-05-09
POD类型
2021-05-09
const与常量,傻傻分不清楚~
2021-05-09
Head First设计模式——迭代器模式
2021-05-09
MongoDB版本及存储引擎区别
2021-05-09
shell echo单行和多行文字定向写入到文件中
2021-05-09
Java高性能编程之CAS与ABA及解决方法
2021-05-09
HTML5新特性
2021-05-09
cmp命令
2021-05-09
Linux 磁盘管理(df fu fdisk mkfs mount)
2021-05-09
老Python总结的字典相关知识
2021-05-09
jQuery的事件绑定与触发 - 学习笔记
2021-05-09
Linux上TCP的几个内核参数调优
2021-05-09
记一次讲故事机器人的开发-我有故事,让机器人来读
2021-05-09
seo 回忆录百度基本概念(一)
2021-05-09