Hdu 4117 GRE Words （后缀数组+dp）

题目大意：

求出最多能记住的单词的权值和，要求最大。

记住的规则就是上一个单词是这个单词的子串。

思路分析：

首先得声明这题是数据水了才能用sa做的。

sa的复杂度最多可以达到 Orz（sumlen * sumlen） ...

所以我们sa处理的就是这个串是否是下一个串的子串，如果是就转移方程。

dp[i] = max (dp[i] , dp[j] + val[i])...

#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
#define maxn 350005
#define inf 0x3f3f3f3f
using namespace std;

int str[maxn];
int sa[maxn],t1[maxn],t2[maxn],c[maxn],n;
int bel[maxn];
int st[maxn],sumlen[maxn],val[maxn];
int dp[maxn];
void suffix(int m)
{
    int *x=t1,*y=t2;
    for(int i=0;i<m;i++)c[i]=0;
    for(int i=0;i<n;i++)c[x[i]=str[i]]++;
    for(int i=1;i<m;i++)c[i]+=c[i-1];
    for(int i=n-1;i>=0;i--)sa[--c[x[i]]]=i;
    for(int k=1;k<=n;k<<=1)
    {
        int p=0;
        for(int i=n-k;i<n;i++)y[p++]=i;
        for(int i=0;i<n;i++)if(sa[i]>=k)y[p++]=sa[i]-k;
        for(int i=0;i<m;i++)c[i]=0;
        for(int i=0;i<n;i++)c[x[y[i]]]++;
        for(int i=0;i<m;i++)c[i]+=c[i-1];
        for(int i=n-1;i>=0;i--)sa[--c[x[y[i]]]]=y[i];
        swap(x,y);
        p=1;x[sa[0]]=0;
        for(int i=1;i<n;i++)
        x[sa[i]]=y[sa[i-1]]==y[sa[i]]&&y[sa[i-1]+k]==y[sa[i]+k]?p-1:p++;
        if(p>=n)break;
        m=p;
    }
}
int rank[maxn],height[maxn];
void getheight()
{
    int k=0;
    for(int i=0;i<n;i++)rank[sa[i]]=i;
    for(int i=0;i<n;i++)
    {
        if(k)k--;
        if(!rank[i])continue;
        int j=sa[rank[i]-1];
        while(str[i+k]==str[j+k])k++;
        height[rank[i]]=k;
    }
}
char tmp[maxn];
int main()
{
    int T,cas=1;
    for(scanf("%d",&T);T--;)
    {
        int N;
        scanf("%d",&N);
        n=0;
        for(int i=1;i<=N;i++)
        {
            scanf("%s%d",tmp,&val[i]);
            dp[i]=val[i];
            int len=strlen(tmp);
            st[i]=n;
            sumlen[i]=len;
            for(int j=0;j<len;j++)
            {
                str[n]=tmp[j];
                bel[n++]=i;
            }
            bel[n]=i;
            str[n++]=i+128;
        }
        str[n-1]=0;
        suffix(30000);
        getheight();

        for(int i=1;i<=N;i++)
        {
            int mlen=inf;
            for(int j=rank[st[i]]+1;j<n;j++)
            {
                mlen=min(mlen,height[j]);
                if(mlen<sumlen[i])break;
                if(bel[sa[j]]>i)
                {
                    dp[bel[sa[j]]]=max(dp[bel[sa[j]]],val[bel[sa[j]]]+dp[i]);
                }
            }

            mlen=inf;
            for(int j=rank[st[i]]-1;j>0;j--)
            {
                mlen=min(mlen,height[j+1]);
                if(mlen<sumlen[i])break;
                if(bel[sa[j]]>i)
                {
                    dp[bel[sa[j]]]=max(dp[bel[sa[j]]],val[bel[sa[j]]]+dp[i]);
                }
            }
        }
        printf("Case #%d: %d\n",cas++,max(0,*max_element(dp+1,dp+1+N)));
    }
    return 0;
}
/*
1
6
abbab 800
a 1
ab 2
abb 3
baba 5
abbab 8
*/