【HDOJ6583】Typewriter（SAM，DP）

题意：给定一个由小写字母组成的字符串，每次可以花费p在串后加上任意一个字母，花费q在串后复制一个当前串的子串，问生成字符串的最小花费

n<=2e5,1<=p,q<2^31

思路：

 SAM上每个结点表示的串长度为[st[F[p]]+1,st[p]]，如果长度不符合当前i,j的要求，比如复制的时候长度小于当前的一半或者长度大于[j+1,i]就暴力往上跳parent

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef pair<int,int> PII;
typedef pair<ll,ll> Pll;
typedef vector<int> VI;
typedef vector<PII> VII;
typedef pair<ll,int>P;
#define N  500010
#define M  210000
#define fi first
#define se second
#define MP make_pair
#define pi acos(-1)
#define mem(a,b) memset(a,b,sizeof(a))
#define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
#define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
#define lowbit(x) x&(-x)
#define Rand (rand()*(1<<16)+rand())
#define id(x) ((x)<=B?(x):m-n/(x)+1)
#define ls p<<1
#define rs p<<1|1

const ll MOD=1e9+7,inv2=(MOD+1)/2;
      double eps=1e-6;
      ll INF=1ll<<62;
      ll inf=5e13;
      int dx[4]={-1,1,0,0};
      int dy[4]={0,0,-1,1};

char s[N];
int ch[N][26],st[N],F[N],p,np,q,nq,cnt,t;
ll dp[N],A,B;

int read()
{
   int v=0,f=1;
   char c=getchar();
   while(c<48||57<c) {if(c=='-') f=-1; c=getchar();}
   while(48<=c&&c<=57) v=(v<<3)+v+v+c-48,c=getchar();
   return v*f;
}

struct sam
{
    void extend(int x)
    {
        p=np;
        st[np=++cnt]=st[p]+1;
        while(p&&!ch[p][x])
        {
            ch[p][x]=np;
            p=F[p];
        }

        if(!p) F[np]=1;
          else if(st[p]+1==st[q=ch[p][x]]) F[np]=q;
            else
            {
                  st[nq=++cnt]=st[p]+1;
                  memcpy(ch[nq],ch[q],sizeof ch[q]);
                  F[nq]=F[q];
                  F[q]=F[np]=nq;
                  while(p&&ch[p][x]==q)
                  {
                      ch[p][x]=nq;
                      p=F[p];
                  }
            }
    }

    void cancel(int len)
    {
        while(t&&st[F[t]]>=len) t=F[t];
        if(!t) t=1;
    }

    void trans(int len,int x)
    {
        t=ch[t][x];
        if(!t) t=1;
        cancel(len);
    }

}sam;

void solve()
{
    int n=strlen(s);
    //rep(i,1,n) dp[i]=INF;

    cnt=np=t=1;
    sam.extend(s[0]-'a');
    dp[0]=A;
    int j=1;
    rep(i,1,n-1)
    {
        dp[i]=dp[i-1]+A;
        int x=s[i]-'a';
        while((!ch[t][x]||i-j+1>(i+1)/2)&&j<=i)
        {
            sam.extend(s[j++]-'a');
            sam.cancel(i-j);

        }
        sam.trans(i-j+1,x);
        if(j<=i) dp[i]=min(dp[i],dp[j-1]+B);
    }

    rep(i,1,cnt)
    {
        rep(j,0,25) ch[i][j]=0;
        F[i]=st[i]=0;
    }
    printf("%lld
",dp[n-1]);
}

int main()
{
    //freopen("1.in","r",stdin);
    //freopen("1.out","w",stdout);
    while(scanf("%s%lld%lld",s,&A,&B)!=EOF)
    {
        solve();
    }
    return 0;
}