【BZOJ】1068: [SCOI2007]压缩（dp）

http://www.lydsy.com/JudgeOnline/problem.php?id=1068

发现如果只设一维的话无法转移

那么我们开第二维，发现对于前i个来说，如果确定了M在哪里，第i个是用R还是不用就能确定了（如果用R那么在中间一定变成了缓冲串）

那么可以转移了

设d[i,j]表示前i个串，最近的一个M在i的前边一个格子，的最短长度，有

d[1,1]=1

d[i,i]=min{d[i-1,j]}+2 //即用一次M又补上i，所以+2

d[i,j]=d[pos,j]+1，其中pos=(i+j-1)/2，且是整数（就是确定了M找前一个R），且s[j, pos]=s[pos+1, i]

d[i,j]=min{d[i,j], d[k, j]+i-k}

答案就是min{d[n,i], 1<=i<=n}

因为数据小所以直接暴力处理，也就是n^3的，n^2的话应该挺好优化的

#include <cstdio>
#include <cstring>
#include <cmath>
#include <string>
#include <iostream>
#include <algorithm>
#include <queue>
#include <set>
#include <map>
using namespace std;
typedef long long ll;
#define pii pair<int, int>
#define mkpii make_pair<int, int>
#define pdi pair<double, int>
#define mkpdi make_pair<double, int>
#define pli pair<ll, int>
#define mkpli make_pair<ll, int>
#define rep(i, n) for(int i=0; i<(n); ++i)
#define for1(i,a,n) for(int i=(a);i<=(n);++i)
#define for2(i,a,n) for(int i=(a);i<(n);++i)
#define for3(i,a,n) for(int i=(a);i>=(n);--i)
#define for4(i,a,n) for(int i=(a);i>(n);--i)
#define CC(i,a) memset(i,a,sizeof(i))
#define read(a) a=getint()
#define print(a) printf("%d", a)
#define dbg(x) cout << (#x) << " = " << (x) << endl
#define error(x) (!(x)?puts("error"):0)
#define printarr2(a, b, c) for1(_, 1, b) { for1(__, 1, c) cout << a[_][__]; cout << endl; }
#define printarr1(a, b) for1(_, 1, b) cout << a[_] << ‘\t‘; cout << endl
inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<‘0‘||c>‘9‘; c=getchar()) if(c==‘-‘) k=-1; for(; c>=‘0‘&&c<=‘9‘; c=getchar()) r=r*10+c-‘0‘; return k*r; }
inline const int max(const int &a, const int &b) { return a>b?a:b; }
inline const int min(const int &a, const int &b) { return a<b?a:b; }

const int N=52, oo=0x3f3f3f3f;
int d[N][N], n;
char s[N];
bool check(int i, int j) {
	int pos=(i+j-1), k=1;
	if(pos&1) return 0;
	pos>>=1;
	if(d[pos][j]>=oo) return 0;
	for1(x, j, pos) { if(s[x]!=s[pos+k]) return 0; ++k; }
	return 1;
}
int main() {
	scanf("%s", s+1);
	n=strlen(s+1);
	CC(d, oo);
	d[1][1]=1;
	for1(i, 2, n) {
		for1(j, 1, i-1) d[i][i]=min(d[i-1][j], d[i][i]);
		d[i][i]+=2;
		for1(j, 1, i-1) {
			if(check(i, j)) d[i][j]=d[(i+j-1)>>1][j]+1;
			for1(k, j, i-1) d[i][j]=min(d[i][j], d[k][j]+i-k);
		}
	}
	int ans=oo;
	for1(i, 1, n) ans=min(d[n][i], ans);
	print(ans);
	return 0;
}

Description

给一个由小写字母组成的字符串，我们可以用一种简单的方法来压缩其中的重复信息。压缩后的字符串除了小写字母外还可以（但不必）包含大写字母R与M，其中M标记重复串的开始，R重复从上一个M（如果当前位置左边没有M，则从串的开始算起）开始的解压结果（称为缓冲串）。 bcdcdcdcd可以压缩为bMcdRR，下面是解压缩的过程：

另一个例子是abcabcdabcabcdxyxyz可以被压缩为abcRdRMxyRz。

Input

输入仅一行，包含待压缩字符串，仅包含小写字母，长度为n。

Output

输出仅一行，即压缩后字符串的最短长度。

Sample Input

bcdcdcdcdxcdcdcdcd

Sample Output

12

HINT

在第一个例子中，解为aaaRa，在第二个例子中，解为bMcdRRxMcdRR。

【限制】

100%的数据满足：1<=n<=50 100%的数据满足：1<=n<=50

Source