【HDU4552】 怪盗基德的挑战书（后缀数组）

怪盗基德的挑战书

Problem Description

　　“在树最美丽的那天，当时间老人再次把大钟平均分开时，我会降临在灯火之城的金字塔前，带走那最珍贵的笑容。”这是怪盗基德盗取巴黎卢浮宫的《蒙娜丽莎的微笑》这幅画时，挑战书上的内容。
　　但这次，怪盗基德的挑战书上出现了一串串小写字母“aaab sdfeeddd...”。柯南以小学生的眼睛，超凡高中生的头脑，快速统计各种字母频率，字符串长度，并结合挑战书出现的时间等信息，试图分析怪盗基德的意图。最后，他将线索锁定在字符串的循环次数上。并且进一步推理发现，从字符串的第一位开始，到第i位，形成该字符串的子串(c1, c2, c3 ... ci )。对于某一子串ci在该字符串中出现的次数记为ki，则全部子串的循环次数总和AIM = k1 + k2 + ... + ki + ... + kn，柯南发现，AIM恰好对应一个ASCII码！所以，只要把挑战书上的字符串转变成数字，再找到对应的ASCII码，就可以破解这份挑战书了！
　　现在，你的任务就是把字符串转变成对应数字，因为ASCII码以及扩展ASCII码全部只有256个，所以，本题只要把结果对256取余即可。

Input

输入有多组测试数据；
每组测试数据只有一个字符串，由各种小写字母组成，中间无空格。
字符串的长度为L(0 < L <= 100000)。

Output

请计算并输出字符串的AIM值，每组数据输出一行。

Sample Input

aaa

abab

 

Sample Output

6
6

 

 

【题意】

　　给一个串s，求s的每个前缀出现次数之和。

 

【分析】

　　用后缀数组求的话，就是求出每个后缀和最长的后缀的公共前缀长度就可以了。

　　就是rank[1]的位置往两边找。

代码如下：

 #include<cstdio>
 #include<cstdlib>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<queue>
 using namespace std;
 #define INF 0xfffffff
 #define Maxl 100010
 #define Mod 256
 
 int k,la;
 char a[Maxl],b[Maxl];
 int c[Maxl];
 int cl;
 
 int sa[Maxl],rk[Maxl],Rs[Maxl],wr[Maxl],y[Maxl];
 //sa -> 排名第几的是谁
 //rk -> i的排名
 //Rs数值小于等于i的有多少个
 //y -> 第二关键字排名第几的是谁（类似sa）
 int height[Maxl];
 
 int mymin(int x,int y) {return x<y?x:y;}
 
 void get_sa(int m)
 {
    memcpy(rk,c,sizeof(rk));
     for(int i=;i<=m;i++) Rs[i]=;
     for(int i=;i<=cl;i++) Rs[rk[i]]++;
     for(int i=;i<=m;i++) Rs[i]+=Rs[i-];
     for(int i=cl;i>=;i--) sa[Rs[rk[i]]--]=i;
 
     int ln=,p=;
     while(p<cl)
     {
        int k=;
         for(int i=cl-ln+;i<=cl;i++) y[++k]=i;
         for(int i=;i<=cl;i++) if(sa[i]>ln) y[++k]=sa[i]-ln;
         for(int i=;i<=cl;i++) wr[i]=rk[y[i]];
 
         for(int i=;i<=m;i++) Rs[i]=;
         for(int i=;i<=cl;i++) Rs[wr[i]]++;
         for(int i=;i<=m;i++) Rs[i]+=Rs[i-];
         for(int i=cl;i>=;i--) sa[Rs[wr[i]]--]=y[i];
 
         for(int i=;i<=cl;i++) wr[i]=rk[i];
         for(int i=cl+;i<=cl+ln;i++) wr[i]=;
         p=;rk[sa[]]=;
         for(int i=;i<=cl;i++)
         {
             if(wr[sa[i]]!=wr[sa[i-]]||wr[sa[i]+ln]!=wr[sa[i-]+ln]) p++;
             rk[sa[i]]=p;
         }
         m=p,ln*=;
     }
     sa[]=rk[]=;
 }
 
 void get_he()
 {
     int kk=;
     for(int i=;i<=cl;i++)
     {
         int j=sa[rk[i]-];
         if(kk) kk--;
         while(c[i+kk]==c[j+kk]&&i+kk<=cl&&j+kk<=cl) kk++;
         height[rk[i]]=kk;
     }
 }
 
 void ffind()
 {
     int ans=,minn=INF;
     ans+=cl;
     for(int i=rk[];i>=;i--)
     {
         minn=mymin(minn,height[i]);
         ans=(ans+minn)%Mod;
     }minn=INF;
     for(int i=rk[]+;i<=cl;i++)
     {
         minn=mymin(minn,height[i]);
         ans=(ans+minn)%Mod;
     }
     printf("%d\n",ans);
 }
 
 int main()
 {
     while(scanf("%s",a)!=EOF)
     {
         cl=strlen(a);
         for(int i=;i<cl;i++) c[i+]=a[i]-'a'+;
         get_sa();
         get_he();
         ffind();
     }
     return ;
 }

[HDU4522]

2016-07-17 15:50:30