BZOJ 3097 Hash Killer I
3097: Hash Killer I
Description
这天天气不错,hzhwcmhf神犇给VFleaKing出了一道题: 给你一个长度为N的字符串S,求有多少个不同的长度为L的子串。 子串的定义是S[l]、S[l + 1]、... S[r]这样连续的一段。 两个字符串被认为是不同的当且仅当某个位置上的字符不同。
VFleaKing一看觉得这不是Hash的裸题么!于是果断写了哈希 + 排序。 而hzhwcmhf神犇心里自然知道,这题就是后缀数组的height中 < L的个数 + 1,就是后缀自动机上代表的长度区间包含L的结点个数,就是后缀树深度为L的结点的数量。 但是hzhwcmhf神犇看了看VFleaKing的做法表示非常汗。于是想卡掉他。
VFleaKing使用的是字典序哈希,其代码大致如下: u64 val = 0; for (int i = 0; i < l; i++) val = val * base + s[i] - 'a'; u64是无符号int64,范围是[0, 2^64)。VFleaKing让val自然溢出。 base是一个常量,VFleaKing会根据心情决定其值。 VFleaKing还求出来了base ^ l,即base的l次方,这样就能方便地求出所有长度为L的子串的哈希值。 然后VFleaKing给哈希值排序,去重,求出有多少个不同的哈希值,把这个数作为结果。 其算法的C++代码如下:
typedef unsigned long long u64;
const int MaxN = 100000;
inline int hash_handle(const char *s, const int &n, const int &l, const int &base) { u64 hash_pow_l = 1; for (int i = 1; i <= l; i++) hash_pow_l *= base;
int li_n = 0; static u64 li[MaxN];
u64 val = 0; for (int i = 0; i < l; i++) val = val * base + s[i] - 'a'; li[li_n++] = val; for (int i = l; i < n; i++) { val = val * base + s[i] - 'a'; val -= (s[i - l] - 'a') * hash_pow_l; li[li_n++] = val; }
sort(li, li + li_n); li_n = unique(li, li + li_n) - li; return li_n; }
hzhwcmhf当然知道怎么卡啦!但是他想考考你。
Input
没有输入。
Output
你需要输出一组数据使得VFleaKing的代码WA掉。我们会使用Special Judge检查你的结果的正确性。 输出文件共两行。 第一行两个用空格隔开的数n、l。 第二行是一个长度为n的字符串。只能包含'a'~'z'。 需要保证1 <= n <= 10^5, 1 <= l <= n, 不符合以上格式会WA。 不要有多余字符,很可能导致你WA。
Sample Input
Sample Output
HINT
orz 波兰人 & fotile96 & sillycross
Source
VFleaKing & hzhwcmhf
这是一个很强大的系列。
3097: Hash Killer I base看心情,unsigned long long自然溢出 Submit: 951 Solved: 356
3098: Hash Killer II base看心情,MOD1000000007 Submit: 1597 Solved: 838
3099: Hash Killer III base看心情,MOD双组随机质数 Submit: 605 Solved: 0
昨天做了3098,当时直接把N、L开大,输出100000个随机数,就卡了过去。本人坚信与3097是相似的,于是重复交了一波,全部WA了。本人表示很伤心。
最后发现(以下来自http://blog.csdn.net/regina8023/article/details/43112899):
非常神奇的构造题。
首先明白两点:
1.卡hash的关键在于构造两个不同的串对应的hash值相同。
2.爆u64相当于对2^64这个数取模。
如果base是偶数,那么a.........aaa(>64个a)与ba.......aa(a的数量为前面那么串a的数量-1),这两个串长度相同,hash值相同,显然串是不同的,这样就卡掉了。
如果base是奇数,就比较麻烦了。
看vfk的做法吧:
如果base是奇数的话,现在只考虑a、b两个字母。
a \ b表示a能整除b。(orz 具体数学)
设数学上的函数not(S)表示把字符串S中每个位置的'a'变成'b',把'b'变成'a'后形成的字符串。比如not("ababaa") = "bababb"
strA . strB代表字符串串联。如"娃" . "哈哈" = "娃哈哈"
|str|表示字符串str的长度。
设字符串序列{orzstr[i]},orzstr[1] = "a", orzstr[i] = orzstr[i - 1] . not(orzstr[i - 1])
那么|orzstr[i]| = |orzstr[i - 1]| * 2。显然这是等比数列,得到:|orzstr[i]| = |orzstr[1]| . 2 ^ (i - 1) = 2 ^ (i - 1)
设hash(str)为str的哈希值。
则:
hash(orzstr[i]) = hash(orzstr[i - 1]) * base ^ |not(orzstr[i - 1])| + hash(not(orzstr[i - 1]))
= hash(orzstr[i - 1]) * base ^ (2 ^ (i - 2)) + hash(not(orzstr[i - 1]))
hash(not(orzstr[i])) = hash(not(orzstr[i - 1])) * base ^ (2 ^ (i - 2)) + hash(orzstr[i - 1])
两式相减:
hash(orzstr[i]) - hash(not(orzstr[i]))
= (hash(orzstr[i - 1]) * base ^ (2 ^ (i - 2)) + hash(not(orzstr[i - 1]))) - (hash(not(orzstr[i - 1])) * base ^ (2 ^ (i - 2)) + hash(orzstr[i - 1]))
= (hash(orzstr[i - 1]) - hash(not(orzstr[i - 1]))) * (base ^ (2 ^ (i - 2)) - 1)
这让我们发现,hash(orzstr[i]) - hash(not(orzstr[i]))似乎是个神奇的东西。而我们的目的实际上是要找两个字符串strA, strB使得
hash(strA) % 2^64 = hash(strB) % 2^64
相当与
2^64 \ hash(strA) - hash(strB)
设数列{f[i]},f[i] = hash(orzstr[i]) - hash(not(orzstr[i]))
这样就有:
f[i] = f[i - 1] * (base ^ (2 ^ (i - 2)) - 1)
还是有点不爽啊……我们再设数列{g[i]},g[i] = base ^ (2 ^ (i - 1)) - 1
于是能写成:
f[i] = f[i - 1] * g[i - 1]
则f[i] = f[1] * g[1] * g[2] * ... * g[i - 1]
然后发现一个神奇的事情?
base是奇数,则base的任意正整数次方也一定是奇数。所以对于任意的i必有g[i]为偶数,所以2 ^ (i - 1) \ f[i]
问题是不是结束了呢……发现没有……这样的话我们要使2 ^ 64 \ f[i],至少得让i = 65……然后发现|orzstr[65]|是个天文数字。
发现我们刚才那样分析太坑爹了……
i > 1时有:
g[i] = base ^ (2 ^ (i - 1)) - 1 = (base ^ (2 ^ (i - 2)) - 1) * (base ^ (2 ^ (i - 2)) + 1) = g[i - 1] * 一个偶数
而g[1]显然是偶数吧……
那么4 \ g[2],8 \ g[3]...
也就是说2 ^ i \ g[i]
所以f[i] 实际上有:
(2 ^ 1) * (2 ^ 2) * (2 ^ 3) * ... * (2 ^ (i - 1)) \ f[i]
2 ^ (i * (i - 1) / 2) \ f[i]
当i取12时,就有66个2了哟!
这就是卡base为奇数时的方法。orzstr[12]和not(orzstr[12])即为所求。
而读入中base既有奇数又有偶数,直接在奇数构造的字符串后面加64个a就可以了。
这样看来,忽然发现,2^64远远大于1000000007,所以卡起来不是那么简单。那么,对于3099,想要AC谈何容易。
1 /************************************************************** 2 Problem: 3097 3 User: Doggu 4 Language: C++ 5 Result: Accepted 6 Time:0 ms 7 Memory:916 kb 8 ****************************************************************/ 9 10 #include <cstdio> 11 #include <algorithm> 12 char s[100000]; 13 int main() { 14 printf("%d %d\n",(1<<11)+65,(1<<10)); 15 int now=1;s[1]='a'; 16 for( int i = 1; i <= 11; i++ ) { 17 for( int j = 1; j <= now; j++ ) 18 s[j+now]=s[j]=='a'?'b':'a'; 19 now<<=1; 20 } 21 for( int i = 1; i <= now; i++ ) putchar(s[i]); 22 for( int i = 1; i <= 65; i++ ) putchar('a'); 23 putchar('\n'); 24 return 0; 25 } 26 27
构造
转载于:https://www.cnblogs.com/Doggu/p/bzoj3097.html
BZOJ 3097 Hash Killer I相关推荐
- 字符串——BZOJ 3097: Hash Killer I【构造题,思维题】
BZOJ 3097: Hash Killer I[构造题,思维题] 题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=3097 题意:让你出一组 ...
- BZOJ 3098: Hash Killer II(新生必做的水题)
3098: Hash Killer II Time Limit: 5 Sec Memory Limit: 128 MBSec Special Judge Submit: 1555 Solved: ...
- BZOJ 3098(Hash Killer II-生日攻击)
3098: Hash Killer II Time Limit: 5 Sec Memory Limit: 128 MBSec Special Judge Submit: 66 Solved ...
- 【BZOJ 3098】 Hash Killer II
Description 这天天气不错,hzhwcmhf神犇给VFleaKing出了一道题: 给你一个长度为N的字符串S,求有多少个不同的长度为L的子串. 子串的定义是S[l].S[l + 1].... ...
- BZOJ3098. Hash Killer II(生日攻击)
Description 这天天气不错,hzhwcmhf神犇给VFleaKing出了一道题: 给你一个长度为N的字符串S,求有多少个不同的长度为L的子串. 子串的定义是S[l].S[l + 1].- S ...
- hash进阶:使用字符串hash乱搞的姿势
前言 此文主要介绍hash的各种乱搞方法,hash入门请参照我之前这篇文章 不好意思hash真的可以为所欲为 在开头先放一下题表(其实就是我题解中的hash题目qwq) 查询子串hash值 必备的入门 ...
- 2019.4.summary
2019.4.1 BZOJ1061: [Noi2008]志愿者招募 真心有点难QAQ https://www.byvoid.com/zhs/blog/noi-2008-employee 看void爷的 ...
- 中高级数论 [欧拉函数线性筛,二次剩余]
欧拉函数线性筛 对于素数ppp, φ(p∗i)={p−1i=1p∗φ(i)p∣i(p−1)∗φ(i)p∤i\varphi (p*i)= \begin{cases} p-1& i=1\\ p*\ ...
- [BZOJ 3207] 花神的嘲讽计划Ⅰ【Hash + 可持久化线段树】
题目链接:BZOJ - 3207 题目分析 先使用Hash,把每个长度为 k 的序列转为一个整数,然后题目就转化为了询问某个区间内有没有整数 x . 这一步可以使用可持久化线段树来做,虽然感觉可以有更 ...
- BZOJ.4598.[SDOI2016]模式字符串(点分治 Hash)
LOJ BZOJ 洛谷 点分治.考虑如何计算过\(rt\)的答案. 记\(pre[i]\)表示(之前的)子树内循环匹配了\(S\)的前缀\(i\)的路径有多少,\(suf[i]\)表示(之前的)子树内 ...
最新文章
- The Road to learn React书籍学习笔记(第三章)
- minui 向div放html,c# – 如何添加css类到html通用控件div?
- [RQNOJ313]波浪数
- 解决问题 inner element must either be a resource reference or empty.
- FD.io/VPP — VNF 应用场景
- kswapd0 挖矿_bioset linux_linux bioset 进程 腾讯云
- android cm 老罗,Android之父打造了一款全面屏旗舰 罗永浩如此评价
- (转) 淘淘商城系列——解决KindEditor上传图片浏览器兼容性问题
- baq在聊天中啥意思_职场中的“老实人”如何实现逆袭,得到领导的重用?
- vue 利用axios请求接口下载excel
- CSS选择器、优先级与匹配原理
- 删除Windows网络连接方法
- 小程序样式写了没有用,或许你就差一行代码
- Chrome谷歌浏览器插件的本地存放位置
- 【vue+HT+flyTo】HT图扑软件中的flyTo应用
- java递归处理单位人员组织机构树
- HDU oj wod sticks
- 【Nginx】Nginx主机域名配置
- 文件恢复 文件改名后如何恢复原来名称分享操作步骤
- opencart 3.0 版本数据库数据表字典(详细篇)
热门文章
- 记一次wkhtmltopdf填坑经历
- MC2D v0.0.6 Source Code
- SE96X、CM500-I27条码扫描引擎使用经验
- Ubuntu中ls之后文件的颜色含义
- 实战NLP beginner任务一:基于机器学习的文本分类
- html网页表格相同行自动合并,jQuery实现HTML表格单元格的合并功能
- 摄影_光圈、快门、曝光度(ISO)
- 【WIN10】PS/2 标准键盘鼠标带感叹号(代码10,39)(适用于红米G)
- (转帖)SpringBoot自定义Starter
- usb万能驱动win7_突破封锁,倒行逆施强行安装Win7系统,繁琐变简洁