BZOJ4650: [Noi2016]优秀的拆分(hash 调和级数)
时间:2022-06-16
本文章向大家介绍BZOJ4650: [Noi2016]优秀的拆分(hash 调和级数),主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
题意
Sol
NOI的题都这么良心么。。
先交个(n^4)暴力 => 75
hash优化一下 => 90
然后(90)到(100)分之间至少差了(10)难度台阶= =
(90)分的暴力hash就比较trival了。
考虑怎么优化。 显然我们只要找出所有形如(AA)的字符串就行了,设(pre[i])表示以(i)为端点,向前的所有(AA)的数量,(suf[i])表示以(i)为端点,向后的所有(AA)的数量
这样最终答案就是(sum_{i = 1}^{N - 1} pre[i] * suf[i + 1])
那怎么求(pre)呢?((suf)同理,只要把原串翻转一下就和pre一样了)
首先枚举一个长度(len),然后每隔(len)个点打一个标记。
比如(abcabca)在长度为(len = 3)的时候是这样的(abc|abc|a)
对于相邻的两个标记,我们二分找出他们的(LCS)和(LIS),然后考虑在第一个标记左端点的所有点的贡献,一个显然的结论是:(其实也不是很显然,自己举几个例子试试吧)
在([(i - pre + 1, min(i, i - pre + 1 + (pre + suf - len )) + 1]))内的点会产生贡献(可能不是很严格,但是可以A。。)
然后暴力加就可以A了(因为数据水。。)
实际上可以直接差分一波
时间复杂度:(O(nlog^2n))
mmp居然卡unsigned 自然溢出
#include<bits/stdc++.h>
#define ull unsigned long long
#define LL long long
using namespace std;
const int MAXN = 1e5 + 10, mod = 998244353;
inline LL read() {
char c = getchar(); int x = 0, f = 1;
while(c < '0' || c > '9') {if(c == '-') f = -1; c = getchar();}
while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
return x * f;
}
int N, sum1[MAXN], sum2[MAXN];
ull ha[MAXN], base = 123, po[MAXN];
char s[MAXN];
ull gethash(int l, int r) {
if(l <= 0) return ha[r];
return (ha[r] - ha[l - 1] * po[r - l + 1] % mod + mod) % mod;
}
void Get(char *a, int *sum) {
po[0] = 1; ha[0] = 0; memset(ha, 0, sizeof(ha));
for(int i = 1; i <= N; i++) ha[i] = (ha[i - 1] * base % mod + a[i]) % mod, po[i] = base * po[i - 1] % mod;
for(int len = 1; len <= N / 2 + 1; len++) {
for(int i = len; i < N; i += len) {
int ll = i, rr = i + len, l = 0, r = 0, pre = 0, suf = 0, ans = 0;
if(rr > N) continue;
l = 1, r = len; ans = 0;
while(l <= r) {
int mid = l + r >> 1;
if(gethash(ll - mid + 1, ll) == gethash(rr - mid + 1, rr)) l = mid + 1, ans = mid;
else r = mid - 1;
}
pre = ans;
l = 1, r = len; ans = 0;
while(l <= r) {
int mid = l + r >> 1;
if(gethash(ll + 1, ll + mid) == gethash(rr + 1, rr + mid)) l = mid + 1, ans = mid;
else r = mid - 1;
}
suf = ans;
if(pre + suf < len) continue;
//for(int j = i - pre + 1; j <= min(i, i - pre + 1 + (pre + suf - len )); j++) sum[j]++;
sum[i - pre + 1]++;
sum[min(i, i - pre + 1 + (pre + suf - len )) + 1]--;
}
}
for(int i = 1; i <= N; i++) sum[i] += sum[i - 1];
}
void solve() {
memset(sum1, 0, sizeof(sum1));
memset(sum2, 0, sizeof(sum2));
scanf("%s", s + 1); N = strlen(s + 1);
Get(s, sum1);
reverse(s + 1, s + N + 1);
Get(s, sum2);
reverse(sum2 + 1, sum2 + N + 1);
LL ans = 0;
for(int i = 1; i <= N - 1; i++) ans += sum1[i + 1] * sum2[i];
cout << ans << endl;
}
signed main() {
//freopen("a.in", "r", stdin);
for(int T = read(); T; T--, solve());
return 0;
}
- WordPress强迫症技巧:让文章(ID)地址完美连续(障眼法)
- iOS内存管理:从MRC到ARC实践
- MySQL错误修复:Table xx is marked as crashed and last (automatic?) repair failed
- PHP跨站脚本攻击(XSS)漏洞修复方法(一)
- Windows下获取网络连线实际名称,加强IP类设置脚本的兼容性
- Android APP 快速 Pad 化实现
- PHP彩蛋还是漏洞?expose_php彩蛋的触发和屏蔽方法
- 深入源码探索 ReactNative 通信机制
- PHP跨站脚本攻击(XSS)漏洞修复思路(二)
- WordPress发布文章自动同步到新浪微博(带特色图片)
- go http 服务器编程(1)
- Linux系统内存监控、性能诊断工具vmstat命令详解
- go http 服务器编程(2)
- 利用placeholder属性来添加输入框默认文字提示,提高用户体验
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法
- 谱聚类的python实现
- 基于TypeScript封装Axios笔记(三)
- hadoop伪分布式之启动HDFS并运行MR程序(WordCount)
- 细品Reids的HyperLogLog数据结构
- 布隆过滤器与缓存击穿
- spring之通过注解方式配置Bean(二)
- 【python-leetcode103-树的宽度遍历】二叉树的锯齿形层次遍历
- 图卷积网络(GCN)python实现
- hadoop伪分布式之配置yarn并运行MR程序(WordCount)
- 【python-leetcode637-树的宽度遍历】二叉树的层平均值
- graphSAGE的python实现
- 【python-leetcode111-树的宽度遍历】二叉树的最小深度
- 基于TypeScript封装Axios笔记(四)
- spring之泛型依赖注入
- 【python-leetcode113-树的深度遍历】路径总和Ⅱ