2019 杭电多校1 - 1006 Typewriter
题目链接:hdu 6583
题解:
记
f
[
i
]
f[i]
f[i] 为输出前
i
i
i 个字符的最小代价,
f
[
i
]
f[i]
f[i]是非递减的(反证法)。
对于
i
i
i 从小到大处理,维护使得
s
[
j
+
1...
i
]
∈
s
[
0...
j
]
s[j+1 ... i] ∈ s[0 ... j]
s[j+1...i]∈s[0...j] 的最小的
j
j
j ,则
f
[
i
]
=
m
i
n
(
f
[
i
−
1
]
+
p
,
f
[
j
]
+
q
)
f[i] = min(f[i-1]+p , f[j]+q)
f[i]=min(f[i−1]+p,f[j]+q)。
用 SAM 维护
s
[
0...
j
]
s[0...j]
s[0...j],若
s
[
0...
j
]
s[0...j]
s[0...j] 中包含
s
[
j
+
1...
i
]
s[j+1...i]
s[j+1...i],即加入第
i
i
i 个字符仍然能复制,就不需要做任何处理,不需要在SAM中加入
s
[
j
+
1
]
s[j+1]
s[j+1]。否则,重复地将第
j
+
1
j+1
j+1 个字符加入后缀自动机并
j
=
j
+
1
j = j + 1
j=j+1,相应维护
s
[
j
+
1...
i
]
s[j+1...i]
s[j+1...i] 在后缀自动机上新的匹配位置,直到
s
[
j
+
1...
i
]
∈
s
[
0...
j
]
s[j+1... i] ∈ s[0...j]
s[j+1...i]∈s[0...j]。
代码:
#include <bits/stdc++.h>
#define LL long long
#define LD long double
#define ULL unsigned long long
#define UI unsigned int
#define PII pair<int,int>
#define MPII(x,y) pair<int,int>{x,y}
#define _for(i,j,k) for(int i=j;i<=k;i++)
#define for_(i,j,k) for(int i=j;i>=k;i--)
#define efor(i,u) for(int i=head[u];i;i=net[i])
#define lowbit(x) (x&-x)
#define ls(x) x<<1
#define rs(x) x<<1|1
#define inf 0x3fffffff3f3f3f3fll
using namespace std;
const int maxn = 2e6 + 5;
const int M = 1e9 + 7;
inline int mad(int a,int b){return (a+=b)>=M?a-M:a;}
int len[maxn],link[maxn],cnt,last,nxt[maxn][27];
string s;
void extend(int c){
int np=++cnt,p=last;
len[np]=len[p]+1;
while(p&&!nxt[p][c]){
nxt[p][c]=np;
p=link[p];
}
if(!p){
link[np]=1;
}
else{
int q=nxt[p][c];
if(len[p]+1==len[q]){
link[np]=q;
}
else{
int nq=++cnt;
len[nq]=len[p]+1;
_for(i,0,26) nxt[nq][i]=nxt[q][i];
link[nq]=link[q];
while(p&&nxt[p][c]==q){
nxt[p][c]=nq;
p=link[p];
}
link[np]=link[q]=nq;
}
}
last=np;
}
LL dp[maxn];
int main(){
ios::sync_with_stdio(false);
cin.tie(0);cout.tie(0);
while(cin>>s){
int p,q,lens=s.length();
cin>>p>>q;
cnt=last=1;
link[1]=len[1]=0;
dp[0]=p;
extend(s[0]-'a');
int cur = nxt[1][s[0]-'a'],j=0;
_for(i,1,lens-1){
dp[i]=dp[i-1]+p;
while(1){
while(cur!=1&&len[link[cur]]>=(i-j-1)) cur=link[cur];
if(nxt[cur][s[i]-'a']){
cur = nxt[cur][s[i]-'a'];
break;
}
else extend(s[++j]-'a');
}
dp[i]=min(dp[i],dp[j]+q);
}
_for(i,1,cnt){
_for(j,0,26) nxt[i][j]=0;
}
cout<<dp[lens-1]<<"\n";
}
return 0;
}