UOJ #35. 后缀排序(后缀数组模板题)
2016-11-06 16:17
295 查看
测试一下模板
#include<bits/stdc++.h>
#include<string>
using namespace std;
const int maxn =2e5+10;
struct Suffix
{
int r[maxn];
int sa[maxn],rank[maxn],height[maxn];
int t[maxn],t2[maxn],c[maxn],n,m;
void init(string s)
{
n=s.size();
for(int i=0;i<n;i++) r[i]=(int)s[i];
m=128;//
}
int cmp(int *r,int a,int b,int l) {return r[a]==r[b]&&r[a+l]==r[b+l];}
void build()
{
int i,k,p,*x=t,*y=t2;
r[n++]=0;
for (i=0; i<m; i++) c[i]=0;
for (i=0; i<n; i++) c[x[i]=r[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 (k=1,p=1; k<n; k*=2,m=p)
{
for (p=0,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]]=cmp(y,sa[i-1],sa[i],k)?p-1:p++;
}
n--;
}
void LCP()
{
int i,j,k=0;
for (i=1; i<=n; i++) rank[sa[i]]=i;
for (i=0; i<n; i++)
{
if (k) k--;
j=sa[rank[i]-1];
while (r[i+k]==r[j+k]) k++;
height[rank[i]-1]=k;
}
}
int LCS(string s1,string s2)
{
int len=s1.size();
s1=s1+"$"+s2;
init(s1);
build();
LCP();
int ans=0;
for(int i=2;i<=n;i++)
if((sa[i-1]<len)!=(sa[i]<len)) ans=max(ans,height[i]);
return ans;
}
};
Suffix SA;
string str;
int main()
{
cin>>str;
int len=str.size();
SA.init(str);
SA.build();
for(int i=1;i<=len;i++)
cout<< SA.sa[i]+1<<(i==len?'\n':' ');
SA.LCP();
for(int i=1;i<len;i++)
cout<< SA.height[i]<<(i==len-1?'\n':' ');
return 0;
}
#include<bits/stdc++.h>
#include<string>
using namespace std;
const int maxn =2e5+10;
struct Suffix
{
int r[maxn];
int sa[maxn],rank[maxn],height[maxn];
int t[maxn],t2[maxn],c[maxn],n,m;
void init(string s)
{
n=s.size();
for(int i=0;i<n;i++) r[i]=(int)s[i];
m=128;//
}
int cmp(int *r,int a,int b,int l) {return r[a]==r[b]&&r[a+l]==r[b+l];}
void build()
{
int i,k,p,*x=t,*y=t2;
r[n++]=0;
for (i=0; i<m; i++) c[i]=0;
for (i=0; i<n; i++) c[x[i]=r[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 (k=1,p=1; k<n; k*=2,m=p)
{
for (p=0,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]]=cmp(y,sa[i-1],sa[i],k)?p-1:p++;
}
n--;
}
void LCP()
{
int i,j,k=0;
for (i=1; i<=n; i++) rank[sa[i]]=i;
for (i=0; i<n; i++)
{
if (k) k--;
j=sa[rank[i]-1];
while (r[i+k]==r[j+k]) k++;
height[rank[i]-1]=k;
}
}
int LCS(string s1,string s2)
{
int len=s1.size();
s1=s1+"$"+s2;
init(s1);
build();
LCP();
int ans=0;
for(int i=2;i<=n;i++)
if((sa[i-1]<len)!=(sa[i]<len)) ans=max(ans,height[i]);
return ans;
}
};
Suffix SA;
string str;
int main()
{
cin>>str;
int len=str.size();
SA.init(str);
SA.build();
for(int i=1;i<=len;i++)
cout<< SA.sa[i]+1<<(i==len?'\n':' ');
SA.LCP();
for(int i=1;i<len;i++)
cout<< SA.height[i]<<(i==len-1?'\n':' ');
return 0;
}
相关文章推荐
- UOJ #35. 后缀排序 后缀数组模板题
- Uoj#35. 后缀排序
- Uoj #35. 后缀排序(后缀数组)
- 【UOJ 35】 后缀排序|后缀数组 *2
- 【UOJ #35】后缀排序 后缀数组模板
- 【UOJ 35】 后缀排序|后缀数组 *3
- luogu #3809 【模板】后缀排序(后缀数组)
- UOJ #35. BZOJ 1031 后缀排序 后缀数组模板
- uoj #35. 后缀排序
- UOJ.35.[模板]后缀排序(后缀数组 倍增)
- [省选前题目整理][UOJ 35]后缀排序(后缀数组)
- 【后缀数组】uoj#35. 后缀排序
- uoj#35. 后缀排序 后缀数组
- UOJ #35. 后缀排序[后缀数组详细整理]
- 洛谷.3809.[模板]后缀排序(后缀数组 倍增)
- 后缀数组模板-boj477.新来的小妹妹 & boj477. 田田背课文
- 【后缀数组】后缀数组模板
- 【后缀数组】洛谷P3809模板题
- Codevs 1500 后缀排序(后缀数组)
- 后缀数组——附1(含模板)