POJ - 3581 Sequence(后缀数组)
題目鏈接:點擊查看
題目大意:給出一個數列,保證第一個數一定大于其余的所有數,現在需要讓我們將其分為連續的三段,每一段獨立反轉后保證字典序最小,題目要求輸出反轉后字典序最小的數列
題目分析:因為第一個數一定大于其余的所有數,所以將整個數列反轉后,求排名最小的后綴數組就是第一段的答案了,因為這個性質保證了反轉后的數組中最后一個數是最大的,肯定不可能是以最后一個數字為開頭的后綴,只可能是以前面的數字為開頭的后綴,包含了最后一個數字,所以保證了答案的正確性,那么剩下第二段和第三段該如何分割呢,網上都說了正解的方法,但卻沒說明白為什么,正解就是在反轉后的數列中,去掉已經確定好的第一段后,將剩余的數列復制一遍,然后求后綴數組,通過sa來確定排名最小的后綴就是答案了,為什么可以這樣做呢,因為在去掉第一段后的開頭已經不能保證是最大的了,也就不滿足上面所說的性質了,但是仔細考慮一下,我們需要將剩下的字符串分成兩段并分別反轉,那么意味著反轉后復制一下,在前半段的每個后綴都代表了以每個字母為斷點反轉后的答案了,這個時候就可以直接比較這些后綴的大小,選出排名最小的那個后綴就是本題的答案了
注意,數字的范圍題目沒給出,我看網上大部分代碼都是將原模板中的桶排序替換為了sort,我感覺不妥,為什么不直接離散化呢,一開始我也是固執的不想用離散化,但是自己main函數里寫的已經驗證過很多次是沒有問題的了,交上去還是WA,沒辦法向離散化妥協,改完離散化后立馬AC,不得不說不能亂用別人的模板,出了錯誤都不知道是哪里出的
代碼:
#include<iostream> #include<cstdio> #include<cstdlib> #include<string> #include<ctime> #include<climits> #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=2e5+100;int sa[N]; //SA數組,表示將S的n個后綴從小到大排序后把排好序的 //的后綴的開頭位置順次放入SA中 int t1[N],t2[N],c[N],a[N];int s[N],len;vector<int>v;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;} }int main() { // freopen("input.txt","r",stdin); // ios::sync_with_stdio(false);int n,cut1,cut2;//找出兩個斷點后最后統一輸出scanf("%d",&n);len=n;for(int i=0;i<n;i++){scanf("%d",a+i);v.push_back(a[i]);}sort(v.begin(),v.end());v.erase(unique(v.begin(),v.end()),v.end());for(int i=0;i<n;i++)a[i]=lower_bound(v.begin(),v.end(),a[i])-v.begin()+1;for(int i=n-1;i>=0;i--)s[i]=a[n-i-1];s[len]=0;build_sa(s,len+1,v.size()+1);int pos=0;for(int i=1;i<=len;i++){if(sa[i]>=2){cut1=len-sa[i];break;}}len=0;for(int i=n-1;i>=cut1;i--)s[len++]=a[i];for(int i=n-1;i>=cut1;i--)s[len++]=a[i];s[len]=0;build_sa(s,len+1,v.size()+1);for(int i=1;i<=len;i++){if(sa[i]>=1&&sa[i]<len/2){cut2=cut1+len/2-sa[i];break;}}reverse(a,a+cut1);reverse(a+cut1,a+cut2);reverse(a+cut2,a+n);for(int i=0;i<n;i++)printf("%d\n",v[a[i]-1]);return 0; }?
總結
以上是生活随笔為你收集整理的POJ - 3581 Sequence(后缀数组)的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: POJ - 1226 Substring
- 下一篇: (转)计算几何模板 - kuangbin