字符串整理
Posted hangzz
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了字符串整理相关的知识,希望对你有一定的参考价值。
字符串整理(Manachar,KMP,扩展KMP,ACAM,SA,SAM,最小表示法)
基础
- 字符集(sum):一个字符集是一个建立了全序关系的集合,即任意属于(sum)的元素可以比较,字符集中的元素叫做字符
- 字符串:一个字符串(S)将n个字符顺次排列组成,(n)为(S)的长度,计作(|S|),此文使用的字符串均从1下标开始
- 子串:字符串(S)的子串(S[i...j],i leq j),表示(S)串中从(i)到(j)这一段,也就是顺次排列(S[i],S[i+1],...,S[j])形成的字符串。
- 子序列:字符串(S)的子序列是从(S)中将若干元素提取出来并不改变相对位置形成的序列,即(S[p_1],S[p_2],...,S[p_t]),(1 leq p_1 leq p_1 leq ... leq p_t leq |S|)。
- 后缀:是指从某个位置(i)开始到整个串末尾结束的一个特殊子串。字符串(S)的从(i)开头的后缀表示为(Suffix(S,i)),也就是(Suffix(S,i)=S[i...|S|]),(suf(S,k)表示S后k个字符组成的后缀)。
- 前缀:是指从串首开始到某个位置(i)结束的一个特殊子串。字符串(S)的从(i)结尾的前缀表示为(Preffix(S,i)),也就是(Preffix(S,i)=S[1...i]),(pre(S,k)表示S前k个字符组成的前缀)。
- (Borde)r:(0 le r le |S|),(pre(S,r)=suf(S,r)),(pre(S,r))为s的(Border)
(Border有)传递性,(Border)的(Border)还是原串的(Border) - 周期:(0 le p le |S|),(S[i]=S[i+p],forall i in{1,2,...,|S|-p}),p为s的周期
(pre(S,k))是(S)的(Border)(Leftrightarrow) (|S|-k)是(S)的周期 - 自动机
KMP
- (nwxt[i])表示(s[1..i])的最大的(Border)
- 根据(Border)的传递性,我们记录s最大的(Border)后(mb(S)),s的所有(Border)为(mb(S),mb(mb(S)),...)
- 求得时候运用border传递性与以求信息
- 可以(O(N))进行单串的匹配
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1000001;
char t[MAXN], s[MAXN];
int n, m;
int nxt[MAXN], k, ans, sh[MAXN];
int main() {
scanf("%s", t + 1);
scanf("%s", s + 1);
n = strlen(s + 1); m = strlen(t + 1);
nxt[0] = -1; nxt[1] = 0;
for (int i = 2; i <= n; i++) {
k = nxt[i - 1];
while(k != -1 && s[k + 1] != s[i]) {
k = nxt[k];
}
nxt[i] = k + 1;
}
ans = 0; k = 0;
for (int i = 1; i <= m; i++) {
while(k != -1 && s[k + 1] != t[i]) {
k = nxt[k];
}
k = k + 1;
if(k == n) ans++, sh[ans] = i - n + 1;
}
for (int i = 1; i <= ans; i++) printf("%d
", sh[i]);
for (int i = 1; i <= n; i++) printf("%d ", nxt[i]);
return 0;
}
扩展KMP
- 对于字符串(S),(Z_i)表示(lcp(Suffix(S,i),S))
- 在求解时充分利用已有信息(与manachar有相似之处)
- 做i的时候,记录1..i-1中匹配到的lcp到达最右边的位置
#include <bits/stdc++.h>
const int N = 20000007;
int n, m, z[N], ext[N];
char s[N], t[N];
long long ans;
int min(int x, int y) {
return x < y ? x : y;
}
inline void getz() {
z[1] = n;
ans ^= (1ll * 1 * (z[1] + 1));
for (register int i = 2, l = 0, r = 0; i <= n; i++) {//r 维护当前匹配过的子串中达到的最右边的位置
if(i <= r) z[i] = min(z[i - l + 1], r - i + 1);
while(i + z[i] <= n && s[i + z[i]] == s[z[i] + 1]) ++z[i];
if(i + z[i] - 1 > r) l = i, r = i + z[i] - 1;
ans ^= (1ll * i * (z[i] + 1));
}
}
inline void exkmp() {
for (register int i = 1, l = 0, r = 0; i <= m; i++) {
if(i <= r) ext[i] = min(z[i - l + 1], r - i + 1);
while(i + ext[i] <= m && t[i + ext[i]] == s[ext[i] + 1]) ++ext[i];
if(i + ext[i] - 1 > r) l = i, r = i + ext[i] - 1;
ans ^= (1ll * i * (ext[i] + 1));
}
}
int main() {
std::ios::sync_with_stdio(0);
std::cin.tie(0);
std::cin >> (t + 1);
std::cin >> (s + 1);
n = strlen(s + 1);
m = strlen(t + 1);
ans = 0;
getz(); std::cout << ans << std::endl;
ans = 0;
exkmp(); std::cout << ans << std::endl;
return 0;
}
Manachar
#include <bits/stdc++.h>
using namespace std;
int const maxn = 11000000;
char s[maxn << 1];
int cnt;
int hw[maxn << 1], maxright, mid;
int ans;
inline void read() {
char c = getchar();
s[0] = ‘#‘; s[cnt = 1] = ‘#‘;
while(c > ‘z‘ || c < ‘a‘) c = getchar();
while(c <= ‘z‘ && c >= ‘a‘) s[++cnt] = c, s[++cnt] = ‘#‘, c = getchar();
}
void manacher() {
for (int i = 1; i <= cnt; i++) {
if(i <= maxright) hw[i] = min(hw[(mid << 1) - i], maxright - i + 1);
else hw[i] = 1;
while(s[hw[i] + i] == s[i - hw[i]]) hw[i]++;
if(hw[i] + i - 1 >= maxright) maxright = hw[i] + i - 1, mid = i;
ans = max(ans, hw[i]);
}
}
int main() {
read();
manacher();
printf("%d", ans - 1);
return 0;
}
Trie
#include <bits/stdc++.h>
using namespace std;
const int MAXN = 1e5;
int Trie[MAXN][30], tot = 1;
int en[MAXN];
void Insert(char* str) {
int len = strlen(str), p = 1;
for (int i = 0; i < len; i++) {
int ch = str[i] - ‘a‘;
if(Trie[p][ch] == 0) Trie[p][ch] = ++tot;
p = Trie[p][ch];
}
en[p] = 1;
return;
}
bool search(char* str) {
int len = strlen(str), p = 1;
for (int i = 0; i < len; i++) {
int ch = str[i] - ‘a‘;
p = Trie[p][ch];
if(p == 0) return 0;
}
return en[p];
}
int main() {
return 0;
}
ACAM
#include <bits/stdc++.h>
using namespace std;
int n;
char s[5000010];
struct tr
{
int fail, ch[27];
#define fail(x) trie[x].fail
#define ch(x, y) trie[x].ch[y]
} trie[2000001];
int tot = 1, root = 1;
int match[2000001];
int head[2000001], nxt[2000001], v[2000001], cnt;
int size[2000001];
inline void add(int x, int y)
{
nxt[++cnt] = head[x]; head[x] = cnt; v[cnt] = y;
}
inline int son(int x, int y)
{
return ch(x, y) ? ch(x, y) : son(fail(x), y);
}
inline void insert(char str[], int id)
{
int len = strlen(str + 1);
int p = root;
for (int i = 1; i <= len; i++)
{
int c = str[i] - ‘a‘ + 1;
if(!ch(p, c)) ch(p, c) = ++tot;
p = ch(p, c);
}
match[id] = p;
}
inline void getfail()
{
queue<int> q;
for (int i = 1; i <= 26; i++) ch(0, i) = 1;
q.push(1);
while(!q.empty())
{
int u = q.front();
q.pop();
for (int i = 1; i <= 26; i++)
{
if(ch(u, i))
{
fail(ch(u, i)) = son(fail(u), i);
q.push(ch(u, i));
}
}
}
}
inline void buildfail()
{
for (int i = 2; i <= tot; i++)
{
add(fail(i), i);
}
}
inline void getin(char str[])
{
int p = 1;
int len = strlen(s + 1);
for (int i = 1; i <= len; i++)
{
p = son(p, s[i] - ‘a‘ + 1);
++size[p];
}
}
inline void getans(int p)
{
for (int i = head[p]; i; i = nxt[i])
{
getans(v[i]);
size[p] += size[v[i]];
}
}
int main()
{
scanf("%d", &n);
for (int i = 1; i <= n; i++)
{
scanf("%s", s + 1);
insert(s, i);
}
getfail();
buildfail();
scanf("%s", s + 1);
getin(s);
getans(1);
for (int i = 1; i <= n; i++) printf("%d
", size[match[i]]);
return 0;
}
SA
#include<bits/stdc++.h>
const int N=1000010;
char s[N];
int n,m,num;
int sa[N],rk[N],bac[N],y[N],tmp[N];//sa[i]数组表示后缀排名为i的位置,rk[i]表示后缀[i..n]的排名
int height[N];//height[i]表示rk为i的后缀与rk为i-1的后缀的LCP
void getsa(){
//rk[i]表示位置i的第一关键字(排名)
//二元组(rk[i],rk[i+k])
//初始化基数排序(单字符,其实是单元)
for(int i=1;i<=n;i++)bac[rk[i]=s[i]]++;//bac[i]表示第一关键字小于i的个数
for(int i=2;i<=m;i++)bac[i]+=bac[i-1];//故就先逐个统计,再求前缀和,相当于用桶计数
for(int i=n;i>=1;i--)sa[bac[rk[i]]--]=i;//在第二关键字有序的情况下,对位置按照第一关键字的bac数组,逐个附排名
for(int k=1;k<=n;k<<=1){
num=0;//排名的数量,初始化为单个字符的ASCLL码上限,m=ascll(‘z‘)=122
//y[i]表示第二关键字排名为i的二元组的第一关键字位置(i)
for(int i=n-k+1;i<=n;i++)y[++num]=i;//没有第二关键字的排名最靠前
for(int i=1;i<=n;i++)if(sa[i]>k)y[++num]=sa[i]-k;//利用已有的sa_(k/2)数组,rk靠前的且有第一关键字的,将第一关键字位置sa_(k/2)[i]-k放入y
//在上一轮中,rk_k数组已经求好了,倍增求解,
//接下来利用基数排序,求出sa_k
for(int i=1;i<=m;i++)bac[i]=0;
for(int i=1;i<=n;i++)bac[rk[i]]++;//bac[i]表示第一关键字小于i的个数
for(int i=2;i<=m;i++)bac[i]+=bac[i-1];//故就先逐个统计,再求前缀和,相当于用桶计数
for(int i=n;i>=1;i--)sa[bac[rk[y[i]]]--]=y[i];//在第二关键字有序的情况下,对位置按照第一关键字的bac数组,逐个附排名
memcpy(tmp,rk,sizeof(tmp));//用tmp存一下rk_(k/2),在求rk_k时仍会用到
//其实就是利用sa_k和rk_(k/2)数组求rk_k
rk[sa[1]]=1;num=1;
for(int i=2;i<=n;i++){
if(tmp[sa[i]]==tmp[sa[i-1]]&&tmp[sa[i]+k]==tmp[sa[i-1]+k])rk[sa[i]]=num;
else rk[sa[i]]=++num;
}
if(num==n)break;//已经有了n种排名
m=num;//m表示排名的数量(在桶排中为值域)
}
}
void geth(){
//rk[i]=rk1,sa[rk1-1]=x;
//rk[i-1]=rk2,as[rk2-1]=y;
//height[rk1]=lcp(s[i..],s[x...]);
//height[rk2]=lcp[s[i-1..],s[y...]);
//lcp(s[i..],s[x..])>=lcp(s[i-1...],s[y...])
for(int i=1;i<=n;i++){
int x=sa[rk[i]-1];
int k=std::max(0,height[rk[i-1]]-1);
while(s[i+k]==s[x+k])++k;
height[rk[i]]=k;
}
}
int main(){
scanf("%s",s+1);n=strlen(s+1);m=122;
getsa();for(int i=1;i<=n;i++)printf("%d ",sa[i]);
printf("
");
geth();for(int i=2;i<=n;i++)printf("%d ",height[i]);
return 0;
}
SAM
#include<bits/stdc++.h>
#define fi first
#define se second
#define pb push_back
#define mp make_pair
#define SZ(x) ((int)x.size())
#define ALL(x) x.begin(),x.end()
#define U(i,u) for(register int i=head[u];i;i=nxt[i])
#define rep(i,a,b) for(register int i=(a);i<=(b);++i)
#define per(i,a,b) for(register int i=(a);i>=(b);--i)
using namespace std;
typedef long double ld;
typedef long long ll;
typedef unsigned int ui;
typedef pair<int,int> PII;
typedef vector<int> VI;
template<class T> inline void read(T &x){
x=0; char c=getchar(); int f=1;
while (!isdigit(c)) {if (c==‘-‘) f=-1; c=getchar();}
while (isdigit(c)) {x=x*10+c-‘0‘; c=getchar();} x*=f;
}
template<class T> inline void cmin(T &x, T y){x=x<y?x:y;}
template<class T> inline void cmax(T &x, T y){x=x>y?x:y;}
const int N=2000010;
char s[N];
int bac[N],a[N],cnt[N],n;
ll ans;
struct SAM{
int size,last;
int len[N<<1],link[N<<1];
int to[N<<1][30];
void extend(int c){
int p,cur=++size;
len[cur]=len[last]+1;
for(p=last;p&&!to[p][c];p=link[p])to[p][c]=cur;//情况1
if(!p)link[cur]=1;
else{
int q=to[p][c];
//2-A类
if(len[q]==len[p]+1)link[cur]=q;
//2-B类
else{
int cl=++size;
len[cl]=len[p]+1;
memcpy(to[cl],to[q],sizeof(to[cl]));link[cl]=link[q];
while(p&&to[p][c]==q){
to[p][c]=cl;
p=link[p];
}
link[cur]=link[q]=cl;
}
}
last=cur;
cnt[cur]=1;
}
void build(){
scanf("%s",s+1);last=size=1;n=strlen(s+1);
rep(i,1,n)extend(s[i]-‘a‘);
}
void solve(){
rep(i,1,size)bac[len[i]]++;
rep(i,2,size)bac[i]+=bac[i-1];
rep(i,1,size)a[bac[len[i]]--]=i;//对后缀长度排序
per(i,size,1){
int p=a[i];//排名为i的状态为p
cnt[link[p]]+=cnt[p];
if(cnt[p]>1)cmax(ans,1ll*len[p]*cnt[p]);
}
printf("%lld
",ans);
}
}sam;
int main(){
sam.build();
sam.solve();
return 0;
}
以上是关于字符串整理的主要内容,如果未能解决你的问题,请参考以下文章