【题解】PERIOD - Period [POJ1961] [SP263]
【题解】PERIOD - Period [POJ1961] [SP263]
在进入这道题之前,我们需要了解 kmp 算法
不知道的童鞋可以去看一下Silent_EAG(一个可爱的女孩纸)的讲解。
关于 kmp 算法中 next 数组的周期性质
引理:
对于某一字符串 \(S[1\)~\(i ]\),在它众多的\(next[ i ]\)的“候选项”中,如果存在某一个\(next[ i ]\),使得: \(i\) % \(( i - next[ i ] ) == 0\) ,那么 \(S[ 1\)~ \(( i -next[ i ] ) ]\) 可以为 \(S[ 1\) ~ \(i ]\) 的循环元而 \(i / ( i - next[ i ] )\) 即是它的循环次数 \(K\)。
证明如下:
如图,\(next[ i ] = j\),由定义得红色部分两个子串完全相同。
那么有\(S[ 1\)~\(j ] = S[ m\)~\(i ]\) \(( m = i - next[ i ] )\)。
如果我们在两个子串的前面框选一个长度为 m 的小子串(橙色部分)。
可以得到:\(S[ 1\)~\(m ] = S[ m\)~\(b ]\)。
如果在紧挨着之前框选的子串后面再框选一个长度为 m 的小子串(绿色部分),同样的道理,
可以得到:\(S[ m\)~\(b ] = S[ b\)~\(c ]\)
又因为:\(S[ 1\)~\(m ] = S[ m\)~\(b ]\)
所以:\(S[ 1\)~\(m ] = S[ m\)~\(b ] = S[ b\)~\(c ]\)
如果一直这样框选下去,无限推进,总会有一个尽头。
当满足 \(i\) % \(m == 0\) 时,刚好可以分出 \(K\) 个这样的小子串,
且形成循环\(( K = i / m )\)。
因此我们需要从 \(1\)~\(N\) 扫一遍,判断如果 \(next[ i ]\)
可以整除 \(i\) ,即满足 \(i\) % \(( i - next[ ] ) == 0\) ,那么就可以
肯定\(S[ 1\)~ \(( i - next[ i ] ) ]\)是 \(S[ 1\)~\(i ]\) 的最小循环元,而
\(i / ( i - next[ i ] )\) 即是它的最大循环次数 \(K\) ,直接依次输
出这些 \(i\)和 \(K\) 即可。
那么为什么只判断 \(next[ i ]\) 而不判断 \(next[?]\)呢?
(注:\(next[i]\)和\(next[?]\)表述的意义不同,为方便描
述,这里定义\(next[?]\)为\(next[i]\)的$“候选项”中的某一个)
实际上由这道题可以总结出很多结论:
结论一:
若\(i-next[i]\)可整除\(i\),则\(s[1\)~\(i]\)具有长度为\(i-next[i]\)
的循环元,即\(s[1\)~\(i-next[i]]\)。(前面的一大堆文
字和图片已经给出了这个结论的证明,同时结论一
也是后面推导其他结论的理论基础)
结论二:
若\(i-next[?]\)可整除\(i\),则\(s[1\)~\(i]\)具有长度为\(i-next[?]\)
的循环元,即\(s[1\)~\(i-next[?]]\)。
(用与结论一同样的证明方法可以推导出结论二)
(由此处可知,\(i-next[?]\)想用作循环元要满足的
条件是:\(i-next[?]\)可整除\(i\))。
结论三:
任意一个循环元的长度必定是最小循环元长度的倍数
结论四:
如果\(i-next[i]\)不可整除\(i\),\(s[1\)~\(i-next[?]]\)一定
不能作为\(s[1\)~\(i]\)的循环元。
关于结论四的证明和扩展:
①.如果\(s[1\)~\(i-next[i]]\)不能作为\(s[1\)~\(i]\)循环元,那么
\(s[1\)~\(i-next[?]]\)也都一定不能作为\(s[1\)~\(i]\)的循环元。
(即结论四)
②.如果\(i-next[i]\)不可整除\(i\),\(s[1\)~\(i]\)一定不存在循环元。
③.如果\(i-next[i]\)不可整除\(i\),\(i-next[?]\)也都一定不可整除\(i\)。
④.如果\(s[1\)~\(m]\)是\(s[1\)~\(i]\)的循环元,\(next[i]\)一定为\(i-m\)(\(i-m\)一定为
\(next[i]\))。(在算法竞赛进阶指南上有这么一句话:如果\(s[1\)~\(m]\)为\(s[1\)~\(i]\)的循
环元,\(i-m\)一定是\(next[i]\)的“候选项”,与此处意义略有不同)
⑤.若\(m=i-next[i]\),\(j=next[?]\),\(next[j]=j-m\)。(无论\(m\)可否整除\(i\))
(由④扩展而来)
一些题外话:
关于③的证明,有一个很有趣的想法。
有两个数\(a\),\(c\)和一个数的集合\(b\),且\(b\)与\(a\)有一定的关系(限制)。
已知\(a\)不可整除\(c\),证明\(x(x∈b\))不可整除\(c\)(目前尚未成功)。
虽然表面上看起来并没有什么用,但这种思想把图形匹配转
化为了代数证明。
如果有大佬感兴趣可以思考一下。。。
附:来自某李姓Math大佬。
②③④比较好理解,这个⑤是个什么意思呢?
其实不难懂,通俗点说就是\(i-next[?]\)一定是在\(m\)的倍数处
\((m,2m,3m...)\),如果有循环,也可以说是\(i-next[?]\)一定在循环
节点上,或者说是一定在我们先前图片中框选的黑色块的边界相
邻处,不可能在某个黑色块的中间(如图红色为不可能的情况)
注意一下这个等式:\(i-next[j]=i-j+m\)
可以化简为:\(next[j]=j-m\)
那么可以发现每个\(next[?]\)和\(next[next[?]]\)之间刚好相差m,
只是要由⑤推导①的话,用化简前的样子似乎会更好懂一些。
假如④⑤得证,那么它们和①有什么关系呢?
如果\(i-next[?]\)一定是在\(m\)的倍数处\((m=i-next[i])\),
因为当\(m\)不可整除\(i\)时,\(m\)的倍数也不可整除\(i\),
所以\(i-next[?]\)均不满足作为\(s[1\)~\(i]\)循环元的条件(前面已
提到过“条件”具体指什么)。
因此,⑤\(→\)①得证。
如何证明④或者⑤?
如图,\(j=next[i]\),\(m=i-next[i]\)
先按照与之前相同的方法先将\(s[1\)~\(i]\)划分成\(K\)个黑色块
\(j0=next[j]\),\(n=i-next[j]\),假设n不在m的倍数处,如图红色。
同样的,框选出红色块。
然后再作一些辅助线。接下来就开始推理。
设\(v=j-j0\)。
先看左边:\(s[1\)~\(1+v]=s[m\)~\(m+v]\),\(s[1+v\)~\(1+2v]=s[m+v\)~$m+2v] $
再看右边: \(s[1\)~\(1+v]=s[m+v\)~\(m+2v]\)
综合可得:\(s[1\)~\(1+v]=s[m\)~\(m+v]=s[m+v\)~\(m+2v]=s[1+v\)~\(1+2v]\)
无限的推进,再推进,辅助线划分出的长度为v的区域全部相等,直至边界。而此时的边界出现了两种情况:
⑴v可整除i。
此时刚好将\(s[1\)~\(i]\)分成了若干个完全相同的长度为\(v\)的小块,明显形成了循环元\(s[1\)~\(v]\),那么\(next[i]\)至少应为\(i-v\),这与之前的\(next[i]=j\)相矛盾。
⑵v不可整除i。
观察下列图片,你发现了什么?
将蓝圈处放大,发现了一种交叉相等的情况(如图绿色处)。
再把它压扁,并取几个新的名字\(1',m',j',i'\)。此时它变得和初始
的情况一模一样,于是经过相同的操作后,再一次使出了无限
推进,假如每次的\(v'\)都不可整除\(m'\),那么就一路推了边界:\(v'=1\)。
\(1\)可以整除任何数,于是\(s[1\)~\(i]\)形成了长度为\(1\)的循环元,矛盾。
当n不在m的倍数处时,一定会出现矛盾,所以假设不成立。
因此④得证。同理⑤也得证。
完结附代码:
#include<cstdio>
int t,i,j,n,nex[1000005];char a[1000005];
int main(){
while(scanf("%d",&n),n){
scanf("%s",a+1);
printf("Test case #%d\n",++t);
for(i=2,j=0;i<=n;i++){//最基本的 next[] 数组求法
while(j&&a[i]!=a[j+1])j=nex[j];
if(a[i]==a[j+1])j++;nex[i]=j;
}
for(i=2;i<=n;i++)//由于1~1只有一个字母,只能是它本身构成长度为1的循环,所以从2开始枚举
if(i%(i-nex[i])==0&&nex[i])//判断时还要注意nex[i]是否为0
printf("%d %d\n",i,i/(i-nex[i]));
//如果i含有长度大于1的最小循环元,输出i的长度(即i)以及最大循环次数K(即i-nex[i])
printf("\n");//记得输出一个空行
}
}
写了这么多证明,结果最后代码简单得不要不要的。。。。。
最新文章
- linux下解压
- 解决对含有第三方jar包的项目打包出现java.lang.NoClassDefFoundError问题
- linux常用命令 (mac ),积少成多
- git 给远程库 添加多个url地址
- sqoop导入数据到hive
- Html笔记(七)表单
- chmod -R o+rX /data
- OC封装的TLV数据格式解析库
- php一些函数及方法...
- c#简单易用的短信发送服务 悠逸企业短信服务
- MySQL优化原理
- 【Java学习笔记之二十】final关键字在Java继承中的用法小结
- ABP之模块系统
- 马凯军201771010116《面向对象与程序设计Java》第十二周学习总结
- Linux用户、用户组、文件权限学习笔记
- Spring循环依赖
- UIPanGestureRecognizer translateInView, locationInView
- mongodb及mongoclient在win7下的编译和使用
- solr报错 ERROR SolrDispatchFilter null:ClientAbortException: java.net.SocketException: Broken pipe 原因是nginx截断了请求
- js urlencode
热门文章
- SpringBoot+logback实现按业务输出日志到不同的文件
- Asp.Net SignalR 使用记录 技术回炉重造-总纲 动态类型dynamic转换为特定类型T的方案 通过对象方法获取委托_C#反射获取委托_ .net core入门-跨域访问配置
- SpringCloud学习第三章-springcloud 父项目创建
- cocoapods升级
- JavaScript 关于金额、数字验证的正则表达式
- request有get,post,put,delete等方法大全
- 201871010118-唐敬博《面向对象程序设计(Java)》第四周学习总结
- django 解析上传xls文件
- 禁用wordpress模板默认样式
- application platform as a service (aPaaS)