HDU6583:Typewriter(dp+后缀自动机)
Posted heyuhhh
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了HDU6583:Typewriter(dp+后缀自动机)相关的知识,希望对你有一定的参考价值。
题意:
给出\(p,q\),现在要你生成一个字符串\(s\)。
你可以进行两种操作:一种是花费\(p\)的代价随意在后面添加一个字符,另一种是花费\(q\)的代价可以随意赋值前面的一个子串。
现在问最小代价是多少。
思路:
考虑\(dp\),那么就有转移方程:\(dp[i]=min\dp[i-1]+p,dp[j]+q\\),即直接对两种操作取min即可。
注意到该\(dp\)方程有一个性质:其值肯定为单调不降的。因为如果有\(dp[i-1]>dp[i]\),那此时我们跟着复制过来肯定更好,不能复制就花\(p\)在后面增加,那值岂不是也变大了?
所以转移方程中的\(j\)应为最远的一个\(j\)。
之后思路如下:
- 考虑维护离当前\(i\)最远的一个\(j\),使得\(s[1,\cdots,j-1]\)的子串中含有\(s[j,\cdots,i]\)。
- 注意到\(j\)是单调不减的,考虑用后缀自动机来维护这样一个位置。
- 当向\(i+1\)转移时,若目前后缀自动机中无法成功向\(s[i+1]\)转移,那么将第\(j\)个字符加入后缀自动机并且\(j++\),直到成功转移。
- 注意我们始终要保证后缀自动机中的状态长度尽可能小,这样我们才能保证\(j\)尽可能远,详见\(withdraw\)操作。
挺好的一道题,需要透彻分析问题的性质。
代码如下:
#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 4e5 + 5;
int p, q;
char s[N];
struct SAM
struct node
int ch[26];
int len, fa;
node()memset(ch, 0, sizeof(ch)), len = 0;
dian[N];
int last, tot, now;
void init(int n)
last = tot = now = 1;
for(int i = 1; i <= 2 * n; i++)
for(int j = 0; j < 26; j++) dian[i].ch[j] = 0;
dian[i].len = 0;
void add(int c)
int p = last;
int np = last = ++tot;
dian[np].len = dian[p].len + 1;
for(; p && !dian[p].ch[c]; p = dian[p].fa) dian[p].ch[c] = np;
if(!p) dian[np].fa = 1;
else
int q = dian[p].ch[c];
if(dian[q].len == dian[p].len + 1) dian[np].fa = q;
else
int nq = ++tot; dian[nq] = dian[q];
dian[nq].len = dian[p].len + 1;
dian[q].fa = dian[np].fa = nq;
for(; p && dian[p].ch[c] == q; p = dian[p].fa) dian[p].ch[c] = nq;
void withdraw(int lens)
while(now && dian[dian[now].fa].len >= lens) now = dian[now].fa;
if(now == 0) now = 1;
void trans(int t, int lens)
now = dian[now].ch[t];
withdraw(lens);
bool match(int t)
return dian[now].ch[t];
A;
ll dp[N];
int main()
ios::sync_with_stdio(false); cin.tie(0);
while(cin >> s + 1)
int n = strlen(s + 1);
cin >> p >> q;
A.init(n);
int l = 2, r = 1;
A.add(s[1] - 'a'); dp[1] = p;
for(int i = 2; i <= n; i++)
++r; int tmp = s[i] - 'a';
dp[i] = dp[i - 1] + p;
while((!A.match(tmp) || r - l + 1 > i / 2) && l <= r)
A.add(s[l++] - 'a');
A.withdraw(r - l);
A.trans(tmp, r - l + 1);
dp[i] = min(dp[i], dp[l - 1] + q);
cout << dp[n] << '\n';
return 0;
以上是关于HDU6583:Typewriter(dp+后缀自动机)的主要内容,如果未能解决你的问题,请参考以下文章
HDU3247 Resource Archiver(AC自动机+BFS+DP)