本题和矩阵相乘问题非常类似。
问题:
定义于字母表∑{a,b,c)上的乘法表如表所示:
依此乘法表,对任一定义于∑上的字符串,适当加括号表达式后得到一个表达式。
例如,对于字符串x=bbbba,它的一个加括号表达式为(b(bb))(ba)。依乘法表,该表达式的值为a。
试设计一个动态规划算法,对任一定义于∑上的字符串x=x1x2…xn,计算有多少种不同的加括号方式,使由x导出的加括号表达式的值为a。
输入:
输入一个以a,b,c组成的任意一个字符串 str。
输出:
计算出的加括号方式数。
分析:
乘法表问题直观理解就是通过加括号使得最终运算结果为a,该问题与矩阵连乘问题类似,矩阵连乘是每一次加括号选择运算量最小的,写成数学表达式有:
而乘法表问题则是计算所有加括号运算结果为a的情况数,并不要求输出加括号方式。那么可以从乘法的最小单元两个符号相乘的所有可能开始,接着在两个符号相乘的基础上计算三个符号相乘的所有可能。直到计算N长度的符号1-N的所有可能。 可以定义一个三维数组 result[n][n][3],n为输入字符串的长度, result[i][j][0]为从字符串中第i个元素到第j个元素的子串表达式值为a的加括号方式数, result[i][j][1]为从字符串中第i个元素到第j个元素的子串表达式值为b的加括号方式数,result[i][j][2]为从字符串中第i个元素到第j个元素的子串表达式值为c的加括号方式数。
设常量a,b,c 分别为 1, 2 ,3 。n 为字符串的长度。
设字符串的第 i 到 第 j 位乘积为 a 的加括号法有 result[i][j][a] 种,
字符串的第 i 到 第 j 位乘积为 b 的加括号法有 result[i][j][b] 种,
字符串的第 i 到 第 j 位乘积为 c 的加括号法有 result[i][j][c] 种。
则原问题的解是: result[1][n][a] 。
设 k 为 i 到 j 中的某一个字符,则对于 k 从 i 到 j :
result[i][j][a] += result[i][k][a] * result[k + 1][j][c] + result[i][k][b] * result[k + 1][j][c] + result[i][k][c] * result[k + 1][j][a];
result[i][j][b] += result[i][k][a] * result[k + 1][j][a] + result[i][k][a] * result[k + 1][j][b] + result[i][k][b] * result[k + 1][j][b];
result[i][j][c] += result[i][k][b] * result[k + 1][j][a] + result[i][k][c] * result[k + 1][j][b] + result[i][k][c] * result[k + 1][j][c];
#include<iostream>
#include<cstdio>
#include<string>
#include<algorithm>
using namespace std;
const int maxn=50;
int result[maxn][maxn][3]={0};
int main()
{
ios::sync_with_stdio(false);
freopen("in.txt","r",stdin);
int len;//字符串的长度
string str;
cin>>str;
len=str.size();
for (int i = 1; i <= len; i++)
{
//初始化
if(str[i-1]=='a') result[i][i][0]=1;
else result[i][i][0]=0;
if(str[i-1]=='b') result[i][i][1]=1;
else result[i][i][1]=0;
if(str[i-1]=='c') result[i][i][2]=1;
else result[i][i][2]=0;
}
for (int r = 2; r <= len; r++)
{//接着从长度为 2 的子字符串计算,直至计算好整串 str
for (int i = 1; i <= len-r+1; i++)
{
int j = i + r - 1;//计算str[i:j]
for (int k = i; k <= j; k++)
{//根据题目中的表,计算加括号法
result[i][j][0] += result[i][k][0] * result[k + 1][j][2] + result[i][k][1] * result[k + 1][j][2] + result[i][k][2] * result[k + 1][j][0];
result[i][j][1] += result[i][k][0] * result[k + 1][j][0] + result[i][k][0] * result[k + 1][j][1] + result[i][k][1] * result[k + 1][j][1];
result[i][j][2] += result[i][k][1] * result[k + 1][j][0] + result[i][k][2] * result[k + 1][j][1] + result[i][k][2] * result[k + 1][j][2];
}
}
}
cout << result[1][len][0] << endl;
return 0;
}