思路:用f[i][j]用来表示[ i, j ]区间中合并之后的最小长度,那么显然f[1][n]就是答案。
dp有两种转移方式:
1.两个区间的答案合并:
2.考虑[ i, j ]这个区间内部合并(可能是两个合并,也可能是3个、4个....相同的循环节合并)
这部分稍微难一些,我们回想一下kmp里学过的,如果len%(len-ne[r])==0,那么len-ne[r]就是最小的循环节长度!这个性质非常重要,在这里就能够用到,如果条件满足我们可以写出如下转移:
知道了思路,代码就不难写了:
#include <bits/stdc++.h>
using namespace std;
#define FOR(i, a, b) for (int i = (a); i <= (b); ++i)
#define pii pair<int,int>
#define int long long
const int N = 105, inf = 0x3f3f3f3f;
char s[N]; int n, f[N][N];
int ne[N];
bool check(int l,int r,int len){
FOR(i,1,r)
if(s[i]!=s[(i-l)%len+l]) return 0;
return 1;
}
int dfs(int l,int r){
if(l==r) return 1;
if(f[l][r]!=-1) return f[l][r]; //记忆化
int&res=f[l][r]=inf;
//1.两部分加起来取最小
FOR(i,l,r-1) res=min(res,dfs(l,i)+dfs(i+1,r));
//2.考虑当前区间可以合并
// 1)先跑个kmp,求出周期长度
int len=r-l+1, dif=l-1; //和普通kmp的“字符位置差值”是l-1
ne[1]=0;
for(int i=2,j=0; i<=len; i++){
while(j && s[i+dif]!=s[j+1+dif]) j=ne[j];
if(s[i+dif]==s[j+1+dif]) j++;
ne[i] = j;
}
// 2)计算出长度和最小周期长度,如果最小周期长度能被(r-l+1)整除,那就是可以合并
int k=len-ne[r-dif]; //记得这里不是ne[r]而是ne[r-dif],因为要把dif值减回来
if(len%k==0) f[l][r]=min(f[l][r], dfs(l,l+k-1));
return res;
}
void solve(){
n=strlen(s+1);
memset(f,-1,sizeof(f));
cout<<dfs(1,n)<<'\n';
}
signed main(){
ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);
while(cin>>(s+1) && s[1]!='*') solve();
}