HDU - 3518 Boring counting(后缀数组)
生活随笔
收集整理的這篇文章主要介紹了
HDU - 3518 Boring counting(后缀数组)
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
題目鏈接:點擊查看
題目大意:給出一個字符串,問有多少個子串出現了兩次以上,計算時彼此不能覆蓋
題目分析:因為數據范圍比較小,我們可以直接暴力,枚舉子串的長度,利用height數組的性質將后綴分為不同的組,在每個組中嘗試尋找答案即可
代碼:
#include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cstring> #include<algorithm> #include<stack> #include<queue> #include<map> #include<set> #include<cmath> #include<sstream> using namespace std;typedef long long LL;const int inf=0x3f3f3f3f;const int N=1e3+100;char str[N];int sa[N]; //SA數組,表示將S的n個后綴從小到大排序后把排好序的 //的后綴的開頭位置順次放入SA中 int t1[N],t2[N],c[N];int rk[N],height[N],len;int s[N];void build_sa(int s[],int n,int m)//n為添加0后的總長 {int i,j,p,*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(j=1;j<=n;j<<=1) {p=0;for(i=n-j;i<n;i++) y[p++]=i;for(i=0;i<n;i++) if(sa[i]>=j) y[p++]=sa[i]-j;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]+j]==y[sa[i]+j]?p-1:p++;if(p>=n) break;m=p;} }void get_height(int s[],int n)//n為添加0后的總長 {int i,j,k=0;for(i=0;i<=n;i++)rk[sa[i]]=i;for(i=0;i<n;i++) {if(k) k--;j=sa[rk[i]-1];while(s[i+k]==s[j+k]) k++;height[rk[i]]=k;} }void solve(int base=128) {build_sa(s,len+1,base);get_height(s,len); }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);while(scanf("%s",str)!=EOF&&str[0]!='#'){len=strlen(str);for(int i=0;i<len;i++)s[i]=str[i]-'a'+1;s[len]=0;solve(30);int ans=0;for(int i=1;i<=len/2;i++){int mmax=-inf,mmin=inf;for(int j=2;j<=len;j++){if(height[j]>=i){mmax=max(mmax,max(sa[j],sa[j-1]));mmin=min(mmin,min(sa[j],sa[j-1]));}else{if(mmax-mmin>=i)ans++;mmax=-inf;mmin=inf;}}if(mmax-mmin>=i)ans++;}printf("%d\n",ans);}return 0; }?
總結
以上是生活随笔為你收集整理的HDU - 3518 Boring counting(后缀数组)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: SPOJ - DISUBSTR Dist
- 下一篇: HDU - 6194 string st