[JSOI2007]字符加密Cipher SA

[JSOI2007]字符加密Cipher

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 7859  Solved: 3410
[Submit][Status][Discuss]

Description

　　喜欢钻研问题的JS同学，最近又迷上了对加密方法的思考。一天，他突然想出了一种他认为是终极的加密办法

：把需要加密的信息排成一圈，显然，它们有很多种不同的读法。例如下图，可以读作：

JSOI07 SOI07J OI07JS I07JSO 07JSOI 7JSOI0把它们按照字符串的大小排序：07JSOI 7JSOI0 I07JSO JSOI07

 OI07JS SOI07J读出最后一列字符：I0O7SJ，就是加密后的字符串（其实这个加密手段实在很容易破解，鉴于这是

突然想出来的，那就^^）。但是，如果想加密的字符串实在太长，你能写一个程序完成这个任务吗？

Input

　　输入文件包含一行，欲加密的字符串。注意字符串的内容不一定是字母、数字，也可以是符号等。

Output

　　输出一行，为加密后的字符串。

Sample Input

JSOI07

Sample Output

I0O7SJ

HINT

对于100%的数据字符串的长度不超过100000。

 

题解：先将于2倍的然后构造后缀数组即可，

满足条件的是SA中<=n，的那些+n-1个点。

 

 #include<cstring>
 #include<cmath>
 #include<cstdio>
 #include<algorithm>
 #include<iostream>

 #define N 200007
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(ch<''||ch>''){if (ch=='-')f=-;ch=getchar();}
     while(ch>=''&&ch<=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return x*f;
 }

 int n,k;
 int s[N],sa[N],rk[N*],A[N],zhi[N];
 int a[N],b[N],cnta[N],cntb[N],tsa[N],height[N];
 char ch[N];

 void Get_SA()
 {
     for (int i=;i<=n;i++)cnta[i]=;
     for (int i=;i<=n;i++)cnta[s[i]]++;
     for (int i=;i<=n;i++)cnta[i]+=cnta[i-];
     for (int i=n;i>=;i--)sa[cnta[s[i]]--]=i;
     rk[sa[]]=;
     for (int i=;i<=n;i++)rk[sa[i]]=rk[sa[i-]]+(s[sa[i]]!=s[sa[i-]]);
     for (int i=;rk[sa[n]]!=n;i<<=)
     {
         for (int j=;j<=n;j++)a[j]=rk[j],b[j]=rk[j+i];
         for (int j=;j<=n;j++)cnta[j]=cntb[j]=;
         for (int j=;j<=n;j++)cnta[a[j]]++,cntb[b[j]]++;
         for (int j=;j<=n;j++)cnta[j]+=cnta[j-],cntb[j]+=cntb[j-];
         for (int j=n;j>=;j--)tsa[cntb[b[j]]--]=j;
         for (int j=n;j>=;j--)sa[cnta[a[tsa[j]]]--]=tsa[j];
         rk[sa[]]=;
         for (int j=;j<=n;j++)
             rk[sa[j]]=rk[sa[j-]]+(a[sa[j]]!=a[sa[j-]]||b[sa[j]]!=b[sa[j-]]);
     }
 }
 int main()
 {
     scanf("%s",ch+);
     int len;
     len=n=strlen(ch+);
     for (int i=;i<=n;i++)A[i]=s[i]=(int)(ch[i]);
     sort(A+,A+len+);
     len=unique(A+,A+len+)-A-;
     for (int i=;i<=len;i++)zhi[i]=A[i];
     for (int i=;i<=n;i++)s[i+n]=s[i]=lower_bound(A+,A+len+,s[i])-A;
     s[n+n]=;
     int num=n;n=*n-;
     Get_SA();
     for (int i=;i<=n;i++)
         if (sa[i]<=num) putchar(char(zhi[s[sa[i]+num-]]));
 }