Codeforces Beta Round #46 (Div. 2) E. Common ancestor

题意

    定义变换：a_i->b_ic_i ：表示把字符串中的a_i 字符变成b_ic_i 字符。再定义两个字符串 s1 s2 的公共祖先 s3：s1 s2 能够由 s3 经过一些变换分别得到。现在给你两个长度不超过 50 的字符串，问你他们的公共祖先中长度最短的是多少，输出这个最短长度。

做法分析

    动态规划，如果我们知道了每个字符串中从第 i 个位置到第 j 个位置能否变成某个特定的字符 ch，那么我们就可以按照第一个字符串的前 i1 个和第二个字符串的前 i2 找最短的公共祖先来划分阶段，显然这样划分的阶段满足最优子结构和无后效性，定义的状态就是：第一个字符串的前 a 个字符和第二个字符串的前 b 个字符的最短公共祖先的长度，于是可以写出如下的状态转移方程：

    f[a][b]=max{ f[i][j]+1, f[a][b] } 其中，i<a, j<b, 第一个字符串的 [i+1, a] 和第二个字符串的 [j+1, b] 经过压缩（逆变换）后能够得到相同的字符

    现在就变成怎么判断第每个字符串的 [L, R] 区间内的子串能否经过压缩得到一个特定的字符 ch

    这也可以用动态规划来做

    我们以区间的长度来划分阶段（不能以起点和终点的长度来划分阶段），定义的状态就是 g[i, i+L-1][ch] ，表示 [i, i+L-1] 能否缩成 ch 这个字符。有点类似传递闭包，状态转移方程如下：

    g[i, L+i-1][ch] = g[i, L+i-1][ch] || (g[i, k][b] && g[k+1, i+L-1][c] && ch->bc) ，其中 ch->bc 表示 bc 能够缩成 ch

    这样，这道题目就解决了。

    搞了一下午，DP 做的还是太少了，没有划分阶段的经验，还需努力进行大量的练习

AC通道

Codeforces Beta Round #46 (Div. 2) E. Common ancestor

参考代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>

using namespace std;

const int INT_INF=0x3fffffff;

bool compress[2][60][60][30];
char s[2][60], substitution[60][10];
int f[60][60], n, len[2];

void init(int id)
{
	len[id]=(int)strlen(s[id]);
	for(int i=0; i<len[id]; i++)
		compress[id][i][i][s[id][i]-'a']=1;
	for(int L=1; L<=len[id]; L++)
		for(int st=0, en=st+L-1; en<len[id]; st++, en++)
		{
			for(int mid=st; mid<en; mid++)
				for(int i=0; i<n; i++)
					if(compress[id][st][mid][substitution[i][3]-'a'] && compress[id][mid+1][en][substitution[i][4]-'a'])
						compress[id][st][en][substitution[i][0]-'a']=1;
		}
}

int main()
{
	scanf("%s%s%d", s[0], s[1], &n);
	for(int i=0; i<n; i++)
		scanf("%s", substitution[i]);
	memset(compress, 0, sizeof(compress));
	init(0), init(1);
	for(int i=0; i<=len[0]; i++)
		for(int j=0; j<=len[1]; j++)
			f[i][j]=INT_INF;
	f[0][0]=0;
	for(int mid0=0; mid0<len[0]; mid0++)
		for(int mid1=0; mid1<len[1]; mid1++)
			for(int st0=mid0; st0<len[0]; st0++)
				for(int st1=mid1; st1<len[1]; st1++)
					for(int ch=0; ch<26; ch++)
						if(compress[0][mid0][st0][ch] && compress[1][mid1][st1][ch])
							f[st0+1][st1+1]=min(f[st0+1][st1+1], f[mid0][mid1]+1);
	int ans=f[len[0]][len[1]];
	if(ans==INT_INF) ans=-1;
	printf("%d\n", ans);
	return 0;
}