【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=0;i<=m;i++) Rs[i]=0;
    for(int i=1;i<=cl;i++) Rs[rk[i]]++;
    for(int i=1;i<=m;i++) Rs[i]+=Rs[i-1];
    for(int i=cl;i>=1;i--) sa[Rs[rk[i]]--]=i;
    
    int ln=1,p=0;
    while(p<cl)
    {
       int k=0;
        for(int i=cl-ln+1;i<=cl;i++) y[++k]=i;
        for(int i=1;i<=cl;i++) if(sa[i]>ln) y[++k]=sa[i]-ln;
        for(int i=1;i<=cl;i++) wr[i]=rk[y[i]];
        
        for(int i=0;i<=m;i++) Rs[i]=0;
        for(int i=1;i<=cl;i++) Rs[wr[i]]++;
        for(int i=1;i<=m;i++) Rs[i]+=Rs[i-1];
        for(int i=cl;i>=1;i--) sa[Rs[wr[i]]--]=y[i];
        
        for(int i=1;i<=cl;i++) wr[i]=rk[i];
        for(int i=cl+1;i<=cl+ln;i++) wr[i]=0;
        p=1;rk[sa[1]]=1;
        for(int i=2;i<=cl;i++)
        {
            if(wr[sa[i]]!=wr[sa[i-1]]||wr[sa[i]+ln]!=wr[sa[i-1]+ln]) p++;
            rk[sa[i]]=p;
        }
        m=p,ln*=2;
    }
    sa[0]=rk[0]=0;
}

void get_he()
{
    int kk=0;
    for(int i=1;i<=cl;i++)
    {
        int j=sa[rk[i]-1];
        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=0,minn=INF;
    ans+=cl;
    for(int i=rk[1];i>=2;i--)
    {
        minn=mymin(minn,height[i]);
        ans=(ans+minn)%Mod;
    }minn=INF;
    for(int i=rk[1]+1;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=0;i<cl;i++) c[i+1]=a[i]-'a'+1;
        get_sa(30);
        get_he();
        ffind();
    }
    return 0;
}

 

2016-07-17 15:50:30