[SCOI2007]压缩(动态规划,区间dp,字符串哈希)

[SCOI2007]压缩

状态:设\(dp[i][j]\)表示前i个字符,最后一个\(M\)放置在\(j\)位置之后的最短字串长度.

转移有三类,用刷表法来实现.

第一种是直接往压缩串后面填字符,这样就是:

\[dp[i+1][j]=min(dp[i+1][j],dp[i][j]+1) \]

另外一种就是往字串里添加\(R\),要满足相邻两个字符串是匹配的,可以用字符串哈希来快速匹,另外下面的\(k\)和\(sz\)要倍增往后跳(因为重复的串长在成倍增加,题目的样例解释的很清楚了).

\[dp[k][j]=min(dp[k][j],dp[k-sz][j]+1) \]

当然还要往后填\(M\)

\[dp[i][i]=min(dp[i][i],dp[i][j]) \]

时间复杂度\(O(n^2logn)\)

代码很好懂.

#include<bits/stdc++.h>
#define ull unsigned long long
#define maxn 55
using namespace std;
const int base=233;
ull ha[maxn],pw[maxn];
char s[maxn];
int n,ans,dp[maxn][maxn];
ull gethash(int l,int r){return ha[r]-ha[l-1]*pw[r-l+1];}
int main()
{
	scanf("%s",s+1);n=strlen(s+1);pw[0]=1;
	for(int i=1;i<=n;i++)ha[i]=ha[i-1]*base+s[i];
	for(int i=1;i<=n;i++)pw[i]=pw[i-1]*base;
	memset(dp,0x3f,sizeof(dp));dp[1][0]=1;
	//dp[i][j]表示前i个字符,上一个M放在j位置之后的最短加密字串
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<i;j++)dp[i][i]=min(dp[i][i],dp[i][j]+1);
		for(int j=0;j<=i;j++)
		{
			dp[i+1][j]=min(dp[i+1][j],dp[i][j]+1);//往后直接添加字母
			if(i==j)continue;
			int sz=i-j;ull haha=gethash(j+1,i);
			for(int k=i+sz;k<=n;k+=sz)//倍增加R
			{
				if(haha==gethash(k-sz+1,k))//哈希匹配字符串
					dp[k][j]=min(dp[k][j],dp[k-sz][j]+1);
				else break;
				haha=haha*pw[sz]+haha;sz<<=1;
			}
		}
	}
	ans=1e9;
	for(int j=0;j<n;j++)ans=min(ans,dp[n][j]);
	cout<<ans<<endl;
	return 0;
}