矩阵连乘源码:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cstring>
#include <climits>
#include <cmath>
#include <cctype>

typedef long long ll;
using namespace std;

int m[110][110];
int a[110];
int n;

int course(int i,int j)
{
if(m[i][j] != -1)//已经算出来这个值,直接返回
{
return m[i][j];
}
if(i == j)
{
return 0;
}
if(i == j-1)
{
m[i][j] = a[i] * a[i+1] * a[i+2];
return m[i][j];
}
int u = course(i,i) + course(i+1,j) + a[i] * a[i+1] * a[j+1];
int k,t;
for(k=i+1; k<j; k++)
{
t = course(i,k) + course(k+1,j) + a[i] * a[k+1] * a[j+1];
if(t < u)
u = t;
}
m[i][j] = u;
return m[i][j];
}

int main()
{
int i;
while(scanf("%d",&n) != EOF)
{
memset(m,-1,sizeof(m));//设置备忘录。刚开始初始值为-1
for(i=1; i<=n+1; i++)
{
scanf("%d",&a[i]);
}
int sum = course(1,n);
printf("%d\n",sum);
}
return 0;
}

题目地址:http://acm.nyist.net/JudgeOnline/problem.php?pid=536

思路:在模板上的简单改动

AC代码:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cstring>
#include <climits>
#include <cmath>
#include <cctype>

typedef long long ll;
using namespace std;

int m[110][110];
int a[110];
int n;

int course(int i,int j)
{
if(m[i][j] != -1)//已经算出来这个值,直接返回
{
return m[i][j];
}
if(i == j)
{
return 0;
}
if(i == j-1)
{
m[i][j] = a[i] * a[i+1] * a[i+2];
return m[i][j];
}
int u = course(i,i) + course(i+1,j) + a[i] * a[i+1] * a[j+1];
int k,t;
for(k=i+1; k<j; k++)
{
t = course(i,k) + course(k+1,j) + a[i] * a[k+1] * a[j+1];
if(t < u)
u = t;
}
m[i][j] = u;
return m[i][j];
}

int main()
{
int i;
while(scanf("%d",&n) != EOF)
{
int b;
memset(m,-1,sizeof(m));//设置备忘录。刚开始初始值为-1
for(i=1; i<=n-1; i++)
{
scanf("%d%d",&a[i],&b);
}
scanf("%d%d",&a[n],&a[n+1]);
int sum = course(1,n);
printf("%d\n",sum);
}
return 0;
}


石子合并

题目地址:http://acm.nyist.net/JudgeOnline/problem.php?pid=737

思路:

区间DP问题(矩阵连乘,石子合并,括号匹配)_i++

从最中间的一条线向右上角斜着算,和矩阵连乘类似,可以用四边形不等式优化

四边形不等式详解地址:http://blog.163.com/dqx_wl/blog/static/2396821452015111133052112/

手动模拟一下程序的大概过程

合并第0堆和第1堆

dp[0][1]=7+6=13

合并第1堆和第2堆

dp[1][2]=6+5=11

合并第0堆到第2堆有2种方法

1.先合并第0堆和第1堆再合并第2堆

13+13+5=31

2.先合并第1堆和第2堆再合并第0堆

11+11+7=29

取最小值,29,合并第0堆到第2堆的最小代价为29,即dp[0][2]=29

照这个思路算下去,先求出区间长度为2的部分,再根据长度为2的部分,求出长度为3的部分,最后求出整个区间的最小值

AC代码:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cstring>
#include <climits>
#include <cmath>
#include <cctype>

typedef long long ll;
using namespace std;

int a[210];
int dp[210][210];
int s[210][210];

int main()
{
int n,i,j,k,r;
while(scanf("%d",&n) != EOF)
{
for(i=1; i<=n; i++)
{
scanf("%d",&a[i]);
}
for(i=1; i<=n; i++)
{
s[i][i] = a[i];
for(j=i+1; j<=n; j++)
{
s[i][j] = s[i][j-1] + a[j];//得出合并区间i到j的结果
}
}
for(r=2; r<=n; r++)//求区间长度为r的最小值
{
for(i=1; i<=n-r+1; i++)//枚举区间长度的开头
{
j = i + r - 1;//区间长度固定,根据开头求出结尾
dp[i][j] = INT_MAX;
for(k=i; k<j; k++)//枚举出合并这段区间以前的最小花费
{
if(dp[i][j] > dp[i][k] + dp[k+1][j])
{
dp[i][j] = dp[i][k] + dp[k+1][j];
}
}
dp[i][j] += s[i][j];//加上这一段合并的花费
}
}
printf("%d\n",dp[1][n]);
}
return 0;
}

题目地址:​​点击打开链接​

思路:注意数组是从a+1开始存的

AC代码:

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <stack>
#include <map>
#include <cstring>
#include <climits>
#include <cmath>
#include <cctype>

typedef long long ll;
using namespace std;

char a[110];
int dp[110][110];

int main()
{
int i,j,k,l;
while(scanf("%s",a+1))
{
memset(dp,0,sizeof(dp));
if(!strcmp(a+1,"end"))
break;
int len = strlen(a+1);
for(l=2; l<=len; l++)
{
for(i=1; i<=len-l+1; i++)
{
j = i + l - 1;
if((a[i] == '(' && a[j] == ')') || ( a[i] == '[' && a[j] == ']'))
{
dp[i][j] = dp[i+1][j-1] + 2;
}
for(k=i; k<j; k++)
{
dp[i][j] = max(dp[i][j],dp[i][k] + dp[k+1][j]);
}
}
}
printf("%d\n",dp[1][len]);
}
return 0;
}