HDU5769 Substring（后缀数组）

链接：http://acm.hdu.edu.cn/showproblem.php?pid=5769

#include <iostream>
#include <stdio.h>
#include <string.h>
#include <string>
#include <algorithm>
using namespace std;
const int MAXN = 1e5+;
int t1[MAXN],t2[MAXN],c[MAXN];
bool cmp(int *r,int a,int b,int l)
{
    return r[a] == r[b] && r[a+l] == r[b+l];
}
void da(int str[],int sa[],int rankk[],int height[],int n,int m)
{
    n++;
    int i, j, p, *x = t1, *y = t2;
    //第一轮基数排序，如果s的最大值很大，可改为快速排序
    for(i = ; i < m; i++)
        c[i] = ;
    for(i = ; i < n; i++)
        c[x[i] = str[i]]++;
    for(i = ; i < m; i++)
        c[i] += c[i-];
    for(i = n-; i >= ; i--)
        sa[--c[x[i]]] = i;
    for(j = ; j <= n; j <<= )
    {
        p = ;
        //直接利用sa数组排序第二关键字
        for(i = n-j; i < n; i++)
            y[p++] = i;//后面的j个数第二关键字为空的最小
        for(i = ; i < n; i++)
            if(sa[i] >= j)
                y[p++] = sa[i] - j;
        //这样数组y保存的就是按照第二关键字排序的结果
        //基数排序第一关键字
        for(i = ; i < m; i++)c[i] = ;
        for(i = ; i < n; i++)
            c[x[y[i]]]++;
        for(i = ; i < m; i++)
            c[i] += c[i-];
        for(i = n-; i >= ; i--)
            sa[--c[x[y[i]]]] = y[i];
        //根据sa和x数组计算新的x数组
        swap(x,y);
        p = ;
        x[sa[]] = ;
        for(i = ; i < n; i++)
            x[sa[i]] = cmp(y,sa[i-],sa[i],j)?p-:p++;
        if(p >= n)break;
        m = p;//下次基数排序的最大值
    }
    int k = ;
    n--;
    for(i = ; i <= n; i++)rankk[sa[i]] = i;
    for(i = ; i < n; i++)
    {
        if(k)k--;
        j = sa[rankk[i]-];
        while(str[i+k] == str[j+k])k++;
        height[rankk[i]] = k;
    }
}
int rankk[MAXN],height[MAXN];
int RMQ[MAXN];
int mm[MAXN];
int best[][MAXN];
char str[MAXN];
int r[MAXN];
int sa[MAXN];
int len;
int hxpos[MAXN];
int xpos[MAXN];
int main()
{
    int t,cas;
    __int64 sum;
    char x[];
    int i,n,cnt;
    scanf("%d",&t);
    for(cas=; cas<=t; cas++)
    {
        sum=;
        scanf("%s",x);
        scanf("%s",str);
        len=strlen(str);
        memset(xpos,-,sizeof(xpos));
        cnt=;
        for(i=; i<len; i++)
        {
            if(str[i]==x[])
                xpos[cnt++]=i;
        }
        for( i = ; i < len; i++)r[i] = str[i]-'a'+;
        r[len] = ;
        da(r,sa,rankk,height,len,);
        cnt=;
        for(i=; i<len; i++)
        {
            while(xpos[cnt]<i&&xpos[cnt]!=-)
                cnt++;
            if(xpos[cnt]==-)
                hxpos[rankk[i]]=len-i+;
            else
                hxpos[rankk[i]]=xpos[cnt]-i+;
        }
        for(i=; i<len; i++)
        {
            sum+=(__int64)(len-sa[i+]-max(height[i+],hxpos[i+]-));
        }
        printf("Case #%d: %I64d\n",cas,sum);
    }
    return ;
}