后缀数组 (hihocoder重复旋律系列)

最新推荐文章于 2019-09-17 11:40:47 发布
最新推荐文章于 2019-09-17 11:40:47 发布
阅读量279 收藏
点赞数
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/hlxng/article/details/75050328
版权

重复旋律1

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#define RG register
#define IL inline
#define pi acos(-1.0)
#define ll long long
using namespace std;

int gi() {
  char ch=getchar(); int x=0;
  while(ch<'0' || ch>'9') ch=getchar();
  while(ch>='0' && ch<='9') {x=10*x+ch-'0';ch=getchar();}
  return x;
}

const int maxn = 20010;

int n,ck,inf=1<<30,ans;
int t1[maxn],t2[maxn],c[maxn],s[maxn],sa[maxn];
int ra[maxn],he[maxn];
int q[maxn];

void get_sa(int m) {
  int i,*x=t1,*y=t2;
  for(i=0; i<m; i++) c[i]=0;
  for(i=0; i<n; i++) c[x[i]=s[i]]++;
  for(i=1; i<m; i++) c[i]+=c[i-1];
  for(i=n-1; i>=0; i--) sa[--c[x[i]]]=i;
  for(int k=1; k<=n; k<<=1) {
    int p=0;
    for(i=n-k; i<n; i++) y[p++]=i;
    for(i=0; i<n; i++) if(sa[i]>=k) y[p++]=sa[i]-k;
    for(i=0; i<m; i++) c[i]=0;
    for(i=0; i<n; i++) c[x[y[i]]]++;
    for(i=1; i<m; i++) c[i]+=c[i-1];
    for(i=n-1; i>=0; i--) sa[--c[x[y[i]]]]=y[i];//这时的sa是以第二关键字的排名的排列的
    swap(x,y);//把x数组赋给y数组
    p=1,x[sa[0]]=0;
    for(i=1; i<n; i++)  
      x[sa[i]]=y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
    if(p>=n) break;
    m=p;
  }
}

void get_he() {
  int i,j,k=0;
  for(i=0; i<=n; i++) ra[sa[i]]=i;
  for(i=0; i<n; i++) {
    if(k) k--;
    j=sa[ra[i]-1];//rank[n]=0,不会出现负下标
    while(s[i+k] == s[j+k]) k++;
    he[ra[i]]=k;
  }
}

int cal() {//单调队列O(n)维护最小值(单增队列)
  int hd=1,ta=0,res=0;
  for(int i=1; i<ck; i++) {
    while(hd<=ta && he[i]<he[q[ta]]) ta--;
    q[++ta]=i;
  }//先将队列扩展到[1,k-1]
  for(int i=ck; i<=n; i++) {//每次逐次往后移动一个1格,求[i-k+1,i]的最值,即为队头
    while(hd<=ta && he[i]<he[q[ta]]) ta--;
    q[++ta]=i;
    while(hd<=ta && q[hd]<=i-ck+1) hd++;
    res=max(res,he[q[hd]]);
  }
  return res;
}

/*
  要点:
  1'先扩展到k-1
  2'维持单调
  3'区间求最值
*/

int main() {
  n=gi(),ck=gi();
  for(int i=0; i<n; i++) s[i]=gi();
  n++;
  get_sa(100);
  n--;
  get_he();
  int ans=cal();
  printf("%d", ans);
  return 0;
}

重复旋律2

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#define RG register
#define IL inline
#define pi acos(-1.0)
#define ll long long 
using namespace std;

int gi() {
  char ch=getchar(); int x=0;
  while(ch<'0' || ch>'9') ch=getchar();
  while(ch>='0' && ch<='9') {x=10*x+ch-'0';ch=getchar();}
  return x;
}

const int maxn = 100010;

int n,m=1000,ans;
int s[maxn],t1[maxn],t2[maxn],c[maxn];
int sa[maxn],he[maxn],ra[maxn];

void get_sa() {
  int i,*x=t1,*y=t2;
  for(i=0; i<m; i++) c[i]=0;
  for(i=0; i<n; i++) c[x[i]=s[i]]++;
  for(i=1; i<m; i++) c[i]+=c[i-1];
  for(i=n-1; i>=0; i--) sa[--c[x[i]]]=i;
  for(int k=1; k<=n; k<<=1) {
    int p=0;
    for(i=n-k; i<n; i++) y[p++]=i;
    for(i=0; i<n; i++)  if(sa[i]>=k) y[p++]=sa[i]-k;
    for(i=0; i<m; i++) c[i]=0;
    for(i=0; i<n; i++) c[x[y[i]]]++;
    for(i=1; i<m; i++) c[i]+=c[i-1];
    for(i=n-1; i>=0; i--) sa[--c[x[y[i]]]]=y[i];
    swap(x,y);
    p=1,x[sa[0]]=0;
    for(i=1; i<n; i++)
      x[sa[i]]=y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
    if(p>=n) break;
    m=p;
  }
}

void get_he() {
  int i,j,k=0;
  for(i=1; i<=n; i++) ra[sa[i]]=i;
  for(i=0; i<n; i++) {
    if(k) k--;
    j=sa[ra[i]-1];//排名在i前面的字符下标j 
    while(s[i+k]==s[j+k]) k++;
    he[ra[i]]=k;
  }
}

bool check(int k) {
  int minsa=1<<30,maxsa=0;
  for(int i=1; i<=n; i++) {
    if(he[i]<k) {
      maxsa=sa[i];
      minsa=sa[i];
    }
    else {
      maxsa=max(sa[i],maxsa);
      minsa=min(sa[i],minsa);
      if(maxsa-minsa>=k) return true;
    }
  }
  return false;
}

int main() {
  n=gi();
  for(int i=0; i<n; i++) s[i]=gi();
  m++,n++,get_sa(),n--;
  get_he();
  int l=0,r=n,mid;
  while(l<=r) {
    mid=(l+r)>>1;
    if(check(mid)) ans=max(ans,mid),l=mid+1;
    else r=mid-1;
  }
  printf("%d", ans);
  return 0;
}

重复旋律3 

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#define RG register
#define IL inline
#define pi acos(-1.0)
#define ll long long 
using namespace std;

int gi() {
  char ch=getchar(); int x=0;
  while(ch<'0' || ch>'9') ch=getchar();
  while(ch>='0' && ch<='9') {x=10*x+ch-'0';ch=getchar();}
  return x;
}

const int maxn = 100100;

int n,m=101,ls1,ls2;
char s1[maxn],s2[maxn];
int s[maxn*2],c[maxn*3],t1[maxn*2],t2[maxn*2];
int sa[maxn*2],ra[maxn*2],he[maxn*2];

void get_sa() {
  int i,*x=t1,*y=t2;
  for(i=0; i<m; i++) c[i]=0;
  for(i=0; i<n; i++) c[x[i]=s[i]]++;
  for(i=1; i<m; i++) c[i]+=c[i-1];
  for(i=n-1; i>=0; i--) sa[--c[x[i]]]=i;
  for(int k=1; k<=n; k<<=1) {
    int p=0;
    for(i=n-k; i<n; i++) y[p++]=i;
    for(i=0; i<n; i++) if(sa[i]>=k) y[p++]=sa[i]-k;
    for(i=0; i<m; i++) c[i]=0;
    for(i=0; i<n; i++) c[x[y[i]]]++;
    for(i=1; i<=m; i++) c[i]+=c[i-1];
    for(i=n-1; i>=0; i--) sa[--c[x[y[i]]]]=y[i];
    swap(x,y);
    p=1,x[sa[0]]=0; 
    for(i=1; i<n; i++)
      x[sa[i]]=y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
    if(p>=n) break;
    m=p;
  }
}

void get_he() {
  int i,j,k=0;
  for(i=0; i<=n; i++) ra[sa[i]]=i;
  for(i=0; i<n; i++) {
    if(k) k--;
    j=sa[ra[i]-1];
    while(s[j+k]==s[i+k]) k++;
    he[ra[i]]=k;
  }
}

int solve() {
  int ans=0;
  for(int i=1; i<=n; i++) {
    if((sa[i]<ls1 && sa[i-1]>ls1) || (sa[i]>ls1 && sa[i-1]<ls1)) {
      ans=max(ans,he[i]);
    } 
  }
  return ans;
}

int main() {
  scanf("%s%s", s1,s2);
  ls1=strlen(s1),ls2=strlen(s2);
  for(int i=0; i<ls1; i++) s[i]=s1[i]-'a'+1;
  s[ls1]=100;
  for(int i=1; i<=ls2; i++)
    s[i+ls1]=s2[i-1]-'a'+1;
  n=ls1+ls2+1;
  n++;
  get_sa();
  n--;
  get_he();
  printf("%d", solve());
  return 0;
}

重复旋律4 

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#define RG register
#define IL inline
#define pi acos(-1.0)
#define ll long long 
using namespace std;

int gi() {
  char ch=getchar(); int x=0;
  while(ch<'0' || ch>'9') ch=getchar();
  while(ch>='0' && ch<='9') {x=10*x+ch-'0';ch=getchar();}
  return x;
}

const int maxn = 100010;

int n,m=27,ans;
char s1[maxn];
int s[maxn],c[maxn],t1[maxn],t2[maxn];
int sa[maxn],ra[maxn],he[maxn];
int dp[maxn][22];

void get_sa() {
  int i,*x=t1,*y=t2;
  for(i=0; i<m; i++) c[i]=0;
  for(i=0; i<n; i++) c[x[i]=s[i]]++;
  for(i=1; i<m; i++) c[i]+=c[i-1];
  for(i=n-1; i>=0; i--) sa[--c[x[i]]]=i;
  for(int k=1; k<=n; k<<=1) {
    int p=0;
    for(i=n-k; i<n; i++) y[p++]=i;
    for(i=0; i<n; i++) if(sa[i]>=k) y[p++]=sa[i]-k;
    for(i=0; i<m; i++) c[i]=0;
    for(i=0; i<n; i++) c[x[y[i]]]++;
    for(i=1; i<m; i++) c[i]+=c[i-1];
    for(i=n-1; i>=0; i--) sa[--c[x[y[i]]]]=y[i];
    swap(x,y);
    p=1,x[sa[0]]=0;
    for(i=1; i<n; i++)
      x[sa[i]]=y[sa[i-1]]==y[sa[i]] && y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
    if(p>=n) break;
    m=p;		 
  }
}

void get_he() {
  int i,j,k=0;
  for(i=1; i<=n; i++) ra[sa[i]]=i;
  for(i=0; i<n; i++) {
    if(k) k--;
    j=sa[ra[i]-1];
    while(s[j+k]==s[i+k]) k++;
    he[ra[i]]=k;
  }
}

void prepare() {
  memset(dp,127/3,sizeof(dp));
  for(int i=1; i<=n; i++)
    dp[i][0]=he[i];
  for(int j=1; j<=20; j++)
    for(int i=1; i+(1<<j-1)<=n; i++)
      dp[i][j]=min(dp[i][j-1],dp[i+(1<<j-1)][j-1]);
}

int lcp(int l, int r) {
  l=ra[l],r=ra[r];//注意取排名
  if(l>r) swap(l,r);
  l++;
  int len=r-l+1;
  int k=(int)(log(len)/log(2));
  return min(dp[l][k],dp[r-(1<<k)+1][k]);
}

int solve() {//纯抄题解
  for(int L=1; L<=n; L++) {
    for(int i=1; i+L<=n; i+=L) {
      int R=lcp(i,i+L);
      ans=max(ans,R/L+1);
      if(i>=L-R%L) ans=max(ans,lcp(i-L+R%L,i+R%L)/L+1);
    }
  }
  return ans;
}

int main() {
  scanf("%s", s1);
  n=strlen(s1);
  for(int i=0; i<n; i++) s[i]=s1[i]-'a'+1;
  n++,get_sa(),n--;
  get_he();
  prepare();
  printf("%d\n", solve());
  return 0;
}



HLX_BLOG
关注
后缀数组
itegel的专栏
11-18 6904
<br />转自：<br />后缀数组（Suffix Array）——理论和思想<br />　　这一篇和下一篇博文我准备写一下我在参加ACM/ICPC期间曾经研究过的后缀数组。关于后缀数组，网上有很多英文资料，但是很多现在的研究结果都是受1991年Udi Manber & Gene Myers的《Suffix arrays: a new method for on-line string searches》中所提的方法的启发，采用倍增的思想。当然，现在有国外学者提到的三分+分治的线性构造方法除外。<br /
hihoCoder #1407 : 后缀数组二·重复旋律2
你真的介意真的在乎真的想要真的渴望的话你一定能够做到
11-08 436
题目分析 如提示的那样，二分长度，然后进行判断即可。。 #include <cmath> #include <cstdio> #include <cstring> #include <iostream> #include <algorithm> using namespace std; const int maxn = 1e5+100;int s[maxn]; int sa[maxn],t[ma
hihocoder 后缀数组
1999
09-17 145
后缀数组一·重复旋律 时间限制:5000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi平时的一大兴趣爱好就是演奏钢琴。我们知道一个音乐旋律被表示为长度为 N 的数构成的数列。 小Hi在练习过很多曲子以后发现很多作品自身包含一样的旋律。旋律是一段连续的数列，相似的旋律在原数列可重叠。比如在1 2 3 2 3 2 1 中 2 3 2 出现了两次。 小Hi想知道一段旋律中...
11.罗穗骞《后缀数组——处理字符串的有力工具》.zip
08-08
后缀数组的构造是通过一系列算法完成的，其中提到的“倍增算法”和“DC3算法”是两种常用的构建方法。倍增算法利用了后缀数组的一些性质，通过逐步增加比较长度来构建，其时间复杂度可以达到O(n log^2 n)，n为字符串...
后缀数组1
08-03
后缀数组是一种用于字符串处理的强大数据结构，它存储了一个字符串的所有后缀按字典序排序后的起始位置。在本文中，作者许智磊主要探讨了后缀数组的基本概念、构建算法以及其在模式匹配和寻找最长回文子串等问题上的...
后缀数组——处理字符串的有力工具.pdf
08-10
后缀数组对于处理单一字符串中的重复子串问题、子串的数量计算、回文子串以及连续重复子串等都极为有效。通过后缀数组，我们可以快速找到字符串中最长的重复子串、计算出不相同子串的个数，甚至检测字符串中的最长...
后缀数组 后缀数组.pdf
09-17
字符串S的后缀是指从任意位置i开始到字符串末尾的部分，对于字符串S，可以构造出一系列后缀，这些后缀根据字典序进行比较。后缀数组SA是一种数据结构，它通过保存字符串所有后缀的起始位置，并按字典序排序这些后缀...
Java后缀数组之求sa数组的实例代码
08-27
"Java后缀数组之求sa数组的实例代码" 后缀数组是一个字符串所有后缀的集合，然后我们根据后缀数组的一些性质就可以实现各种需求。Java后缀数组之求sa数组的实例代码将会对后缀数组的一些基本概念进行详细解释，并...
[知识点]后缀数组
weixin_33795806的博客
08-19 146
// 本文部分内容参照刘汝佳《算法竞赛入门经典训练指南》，特此说明。 [20190129更新！] 终于！时隔多年对这篇文章重新整理了一下，感谢大家提出的建议与意见。 1、前言 　　趁着这几天上午，把后缀数组大致看完了。这个东西本身的概念可能没太大理解问题，但是它所延伸出来的知识很复杂，很多，还有它的两个兄弟——后缀树，后缀自动机，编起来都不是盖的。 2、概念 　　前面曾...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值