今天我们紧接着上次最后的问题继续讲解。
问题描述:
有一只蚂蚁沿着撗坐标走,可以向左也可以向右,每秒只能走一步。现在假设蚂蚁开始在X点,那么经过N秒钟后到达Y,问总共有多少种方式。
对于这个问题最直接的加法就是DFS搜索,但是每次搜索都会有很多多于的搜索,例如下面我手绘的那幅图的上面部分,4步从3走到5,我们可以走3-2-3-4-5,或者3-4-3-4-5,中间我们发现3-4-5这几步实际上是重复的,这也就是问什么效率较低的原因,而动态规划就是利用额外的空间保存中间的部分结果,从而达到高效。
对于动态规划我们所要做的也就是从子问题出发,利用有限的结果直接计算出下一步的结果。图的下面部分就是我的计算步骤。
我们的递归方程为c[i][j] = c[i-1][j-1] + c[i-1][j+1],其i ,j表示经过i步到达j这个位置(即j为上图的横坐标,i即为纵座标(没会出来)),c[i][j]就是在i,j这个位置的值等于左上角,和右上角的值的和。例如上上面从3到5经过4秒钟有4中方法,而4即为3 + 1.经过这番解释估计你已经明白了,下面我给出两种方法的代码。
#include<stdio.h>
#define c 6
#define r 6
const int STEPS =5;
const int POSITION = 3;
const int DESTINATION = 2;
int cnt = 0;
const int N = 5;
int MATRIX[r][c] = {0};
void caculate1(int STEP, int POSITION1)
{
if((STEP == 0) && (POSITION1 == DESTINATION))
{
cnt++;
return;
}
if((STEP < 0) ||(POSITION1 >N) || POSITION1 < 1)
return;
caculate1(STEP -1, POSITION1 + 1);
caculate1(STEP -1, POSITION1 - 1);
}
void caculate2()
{
int row;
int colum;
for(row = 1; row < r; row++)
for(colum = 1; colum < c; colum++)
{
if((colum -1) >=1)
{
MATRIX[row][colum] += MATRIX[row-1][colum -1];
}
if((colum + 1) < c)
{
MATRIX[row][colum] += MATRIX[row-1][colum + 1];
}
}
}
int main()
{
MATRIX[0][POSITION] = 1;
caculate1(STEPS, POSITION);
printf("the number of ways is:%d\n ", cnt);
printf("------------------------------------\n");
caculate2();
int i,j;
for(i = 0; i < r; i++)
for(j = 0; j < c; j++)
{
printf("%d ", MATRIX[i][j]);
if(j == c -1)
printf("\n" );
}
//printf("%d\n", cnt);
}
下图为测试结果:从3到2经过5步有9中方法
下期将带来进年来很火的一道问题:最短字符编辑距离。此问题作为斯坦福《自然语言处理》这门课的一章,在《算法导论》动态规划课后题中也出现过。
,