exkmp(Z函数) 笔记
exkmp 用于求解这样的问题:
求文本串 \(T\) 的每一个后缀与模式串 \(M\) 的匹配长度(即最长公共前缀长度)。特别的,取 \(M=T\),得到的这个长度被称为 \(Z\) 函数。“函数”只是一个叫法,它本质上是个数组...为了好听,后面叫他“\(Z\) 数组” (互联网上的确有人这么叫)
符号(字符串)
\(|S|\) 表示 \(S\) 的长度
\(S[l:r]\) 表示 \(S\) 从 \(l\) 到 \(r\) 的子串。如果 \(l\) 空着,默认为 \(1\);同理 \(r\) 默认为 \(|S|\)。
也就是 \(S[:x]\) 表示 \(S\) 到 \(x\) 结束的前缀,\(S[x:]\) 表示 \(S\) 从 \(x\) 开始的后缀。
\(LCP(S_1,S_2)\) 表示 \(S_1,S_2\) 的 最长公共前缀 (Longest Common Prefix
)
算法讲解
设 \(p_i=LCP(T_i,M)\)
定义从 \(l\) 开始的匹配区间为 \([l,l+p_l-1]\) (设 \(l+p_l-1=r\))
我们枚举处理。假设现在已经求好了 \([1,i-1]\) 的 \(p\) 数组,要求 \(p_i\)。记录一个 最靠后 的匹配区间 \([l,r]\) (\(l<i\),以 \(r\) 靠后为第一关键字,\(l\) 靠后为第二关键字),考虑直接从 \([l,r]\) 中继承点答案来,那很显然一个前提就是 \(i\le r\) (你 \(i\) 在 \(r\) 外面继承啥)
显然,\(p_i\ge LCP(T[i:r],M)\) (因为 \(T[i:r]\) 是 \(T[i:]\) 前缀)
由定义, \([l,r]\) 是最长匹配长度,可知 \(T[l:r]=M[1:r-l+1]\)。
然后现在假如 \(l<i\le r\),那么显然 \(T[i:r]=M[i-l+1:r-l+1]\)
那么 \(LCP(T[i:r],M)=LCP(M[i-l+1:r-l+1],M)\)
简单想一下,\(LCP(A[l:r],A)=min(LCP(A[l:],A),r-l+1)\)
我们要求 \([l,r]\) 子串与整个串的 \(LCP\),可以先求以 \(l\) 开头的整个后缀的与整个串的 \(LCP\),然后和区间长度取 \(min\)。这显然正确。
然后有:
\(LCP(M[i-l+1:r-l+1],M)=min(LCP(M[i-l+1:],M),(r-l+1)-(i-l+1+1))\)
右边的 \(-l+1\) 两个抵消了,就变成 \(r-i+1\)
然后前面是 \(LCP(M[i-l+1:],M)\) 。这不就是 \(M\) 的 \(Z\) 数组的第 \(i-l+1\) 个位置吗!(还记得 \(Z\) 数组的定义吗?)
觉得看字母理解不了的看图(自己画的)(纯鼠标):
红色的部分就是我们推出来的匹配部分。然后现在我们把 \(M\) 移到 \(i\) 开头的位置来匹配,就相当于把 \(M[i-l+1:r-l+1]\) 这一段(红色)移到 \(M\) 的开头处匹配。这一段匹配的长度就是 \(min(Z_{i-l+1},r-i+1)\)。
假设我们现在能求这个 \(Z\) 数组,那么我们已经知道 \(p_i\) 的最小值了 ,就是 \(min(Z_{i-l+1},r-i+1)\) 。从这个位置开始暴力即可。这样就不用每次从 \(1\) 开始匹配了。
求完 \(p_i\) 之后,记得用 \([i,i+p_i-1]\) 更新 \([l,r]\)。
update 2021.01.20:
时间是线性的。原因是,每次成功匹配之后,能直接继承答案的区间便多了一位,复杂度便是线性的。
证明由神仙学长 tzxydby 提供
如何求 Z 数组
我们发现 \(Z\) 数组就是自己和自己匹配的过程。然后我们把上面过程中 \(M\) 换成 \(T\) 即可。
所以我们还是记录一个最靠后的匹配区间 \([l,r]\),然后 \(p_i\) 就相当于 \(Z_i\) 了。
易得:
\(Z_i=min(LCP(M[i-l+1:],M),r-i+1)=min(LCP(T[i-l+1:],T),r-i+1)=min(Z_{i-l+1},r-i+1)\)
求完 \(Z_i\) 之后,记得用 \([i,i+Z_i-1]\) 来更新 \([l,r]\)。
一样,也是从这里开始暴力即可。同理,时间复杂度依然是线性的。
模板
#include <bits/stdc++.h>
using namespace std;
#define N 20000007
#define F(i,l,r) for(int i=l;i<=r;++i)
#define D(i,r,l) for(int i=r;i>=l;--i)
#define Fs(i,l,r,c) for(int i=l;i<=r;c)
#define Ds(i,r,l,c) for(int i=r;i>=l;c)
#define MEM(x,a) memset(x,a,sizeof(x))
#define FK(x) MEM(x,0)
#define Tra(i,u) for(int i=G.Start(u),v=G.To(i);~i;i=G.Next(i),v=G.To(i))
#define p_b push_back
#define sz(a) ((int)a.size())
#define all(a) a.begin(),a.end()
#define iter(a,p) (a.begin()+p)
#define Flandre_Scarlet int
#define IsMyWife main
char _c;
int I()
{
int x=0; int f=1;
while(_c<'0' or _c>'9') f=(_c=='-')?-1:1,_c=getchar();
while(_c>='0' and _c<='9') x=(x<<1)+(x<<3)+(_c^48),_c=getchar();
return (x=(f==1)?x:-x);
}
void Rd(int cnt,...)
{
va_list args; va_start(args,cnt);
F(i,1,cnt) {int* x=va_arg(args,int*);(*x)=I();}
va_end(args);
}
char a[N],b[N];
void Input()
{
scanf("%s%s",a+1,b+1);
}
int z[N];
void Z(char s[]) // 求 Z 函数
{
int n=strlen(s+1);
z[1]=n; F(i,2,n) z[i]=0;
// Z[1]=n 特判,同时也是递推边界
int l=0,r=0;
F(i,2,n)
{
if (i<=r) z[i]=min(z[i-l+1],r-i+1); // 推理出下界
while(i+z[i]<=n and s[i+z[i]]==s[z[i]+1]) ++z[i]; // 暴力
if (i+z[i]-1>=r) l=i,r=i+z[i]-1; // 更新最靠后的匹配位置
}
}
int p[N];
void ExKmp(char s[],char t[])
{
int n=strlen(s+1);
Z(t);
int l=0,r=0;
F(i,1,n)
{
if (i<=r) p[i]=min(z[i-l+1],r-i+1); // 推理出下界
while(i+p[i]<=n and s[i+p[i]]==t[p[i]+1]) ++p[i]; // 暴力
if (i+p[i]-1>r) l=i,r=i+p[i]-1; // 更新最靠后的匹配位置
}
}
void Soviet()
{
ExKmp(a,b);
int n=strlen(a+1),m=strlen(b+1);
long long ans=0;
F(i,1,m) ans^=1ll*i*(z[i]+1);
printf("%lld\n",ans);
ans=0;
F(i,1,n) ans^=1ll*i*(p[i]+1);
printf("%lld\n",ans);
}
Flandre_Scarlet IsMyWife()
{
Input();
Soviet();
getchar();
return 0;
}
exkmp(Z函数) 笔记的更多相关文章
- luogu P5410 模板 扩展 KMP Z函数 模板
LINK:P5410 模板 扩展 KMP Z 函数 画了10min学习了一下. 不算很难 思想就是利用前面的最长匹配来更新后面的东西. 复杂度是线性的 如果不要求线性可能直接上SA更舒服一点? 不管了 ...
- MySQL函数笔记
MySQL函数笔记 日期函数 SELECT t1.xcjyrq, t1.* FROM view_sbxx t1 WHERE t1.syzt ; SELECT t1.xcjyrq, t1.* FROM ...
- JavaScript基础——JavaScript函数(笔记)
avaScript 函数(笔记) JavaScript 是函数式编程语言,在JavaScript脚本中可以随处看到函数,函数构成了JavaScript源代码的主体. 一.定义函数 定义函数的方法有两种 ...
- STL之vector常用函数笔记
STL之vector常用函数笔记 学会一些常用的vector就足够去刷acm的题了 ps:for(auto x:b) cout<<x<<" ";是基于范围的 ...
- numpy函数笔记(持续更新)
numpy函数笔记 np.isin用法 np.isin(a,b) 用于判定a中的元素在b中是否出现过,如果出现过返回True,否则返回False,最终结果为一个形状和a一模一样的数组.(注意:这里的a ...
- 题解-洛谷P5410 【模板】扩展 KMP(Z 函数)
题面 洛谷P5410 [模板]扩展 KMP(Z 函数) 给定两个字符串 \(a,b\),要求出两个数组:\(b\) 的 \(z\) 函数数组 \(z\).\(b\) 与 \(a\) 的每一个后缀的 L ...
- Atcoder Regular Contest 058 D - 文字列大好きいろはちゃん / Iroha Loves Strings(单调栈+Z 函数)
洛谷题面传送门 & Atcoder 题面传送门 神仙题. mol 一发现场(bushi)独立切掉此题的 ycx %%%%%%% 首先咱们可以想到一个非常 naive 的 DP,\(dp_{i, ...
- 前缀函数与Z函数介绍
字符串算法果然玄学=_= 参考资料: OI Wiki:前缀函数与KMP算法 OI Wiki:Z函数(扩展KMP) 0. 约定 字符串的下标从 \(0\) 开始.\(|s|\) 表示字符串 \(s\) ...
- KMP&Z函数详解
KMP 一些简单的定义: 真前缀:不是整个字符串的前缀 真后缀:不是整个字符串的后缀 当然不可能这么简单的,来个重要的定义 前缀函数: 给定一个长度为\(n\)的字符串\(s\),其 \(前缀函数\) ...
随机推荐
- 搭建一个简单的springMVC框架
//新建一个简单的maven项目,选择war包 //web.xml配置 <?xml version="1.0" encoding="UTF-8"?> ...
- Idea mybatis maper接口与mapper.xml文件关联 会根据接口中的方法点在mxl中生成相应sql方法
- MM-RFQ询价报价
(1).询价报价单事务码:ME41/ME42/ME43 需要的主数据:采购组织.供应商.采购组,物料 (2)ME47:维护供应商的报价.可以用项目明细的条件对供应商的报价进行详细设置. (3)供应商价 ...
- 远程控制卡 使用ipmitools设置ipmi
远程控制卡 使用ipmitools设置ipmi 使用DELL的远程控制卡可以方便的管理服务器 在CentOS中可以使用ipmitools管理 IPMI( Intelligent Platform Ma ...
- 聊一聊这个总下载量36039K的XSS-NPM库,是如何工作的?
上篇文章这一次,彻底理解XSS攻击讲解了XSS攻击的类型和预防方式,本篇文章我们来看这个36039K的XSS-NPM库(你没有看错就是3603W次, 36039K次,36,039,651次,数据来自h ...
- 并发编程之JMM&Volatile(一)
并发 很多程序员应该对并发一词并不陌生,并发如同一把双刃剑,如果使用得当,可以帮助我们更好的压榨硬件的性能,反之,也会产生一些难以排查的问题.这里,先简单介绍下并发的几个基本概念. 进程与线程 进程: ...
- Elasticsearch java api操作(二)(Java High Level Rest Client)
一.说明: 一.Elasticsearch提供了两个JAVA REST Client版本: 1.java low level rest client: 低级别的rest客户端,通过http与集群交互, ...
- java操作hive和beeline的使用
一.java操作hive 1.启动服务:hiveserver2,让hive开启与外部连接的服务 nohup hiveserver2 1>/dev/null 2>/dev/null & ...
- Java 中泛型的实现原理
泛型是 Java 开发中常用的技术,了解泛型的几种形式和实现泛型的基本原理,有助于写出更优质的代码.本文总结了 Java 泛型的三种形式以及泛型实现原理. 泛型 泛型的本质是对类型进行参数化,在代码逻 ...
- 【老孟Flutter】为什么 build 方法放在 State 中而不是在 StatefulWidget 中
老孟导读:此篇文章是生命周期相关文章的番外篇,在查看源码的过程中发现了这一有趣的问题,欢迎大家一起探讨. Flutter 中Stateful 组件的生命周期:http://laomengit.com/ ...