从矩阵左上角至右下角(无权值)问题

来源:互联网 发布:mysql更换引擎 编辑:程序博客网 时间:2024/05/21 19:31

问题来源

来自于Euler Project的第十五个题目,题目不复杂,但是可以有很多种思路和解法。此外,我宿舍一姐们面试刚好碰到此题。题目链接

问题描述

一个 M*N 的矩阵,从矩阵的左上角到矩阵的右下角,只能向右走或者向下走,有多少种解法?

例如一个2*2的矩阵,有下面六种解法。


问题解法

假设矩阵有M行N列。

1、数学解法

使用排列组合有两种思路。因为只能向右走或者向下走。

(1)在M+N次行走后,才能到达终点,也就是右下角,而在这M+N次行走中,有M次是向下,N次向右。所以是一个选择问题。


对于Euler这题而言,就是C(40, 20),在google中输入40 choose 20,答案就出来啦。


(2)水平行走记作0,竖直行走记作1。每一种行走足迹可以作为一个0,1串,其中n个0,m个1。可以看做0000000000000(n个0)1111111111111(m个1)的重排列。

也就是


如果实在想不出公式,可以写出前面一些特殊的然后找规律,这也是没有办法的办法了,但是不一定有效。


2、递归解法。

递归,动态规划,说白了,就是把一个问题分解成子问题,然后找出最小子问题。动态规划还需要找出问题的求解顺序。

因此,我们发现了这个问题的递归式:

matrix[i][j]=matrix[i-1][j] + matrix[i][j-1]

以及初始条件:matrix[i][0]=matrix[0][j]=1

[cpp] view plaincopy
  1. #include <stdio.h>  
  2. #include <stdlib.h>  
  3. #include <time.h>  
  4. __int64 digui(int row, int col){  
  5.         if(row == 1 || col == 1) return 1;  
  6.         else return digui(row-1, col) + digui(row, col-1);  
  7. }  
  8. int main(){  
  9.     clock_t begin = clock();  
  10.     printf("%lld\n", digui(18,18));  
  11.     clock_t end = clock();  
  12.     double cost = (double)(end - begin) / CLOCKS_PER_SEC;  
  13.     printf("%lf seconds\n", cost);  
  14.     system("pause");  
  15. }  

但是,递归因为不记录中间结果,会重复计算很多中间值,超级费时!

比如这个程序在我的PC上跑,16*16的时候,需要3.727000秒。而17*17规模,需要13.601000秒,当问题规模到18*18时候计算时间已经到了52.803000秒了,基本到21*21是很慢很慢的了。

3、动态规划解法

和递归的原理差不多,需要多注意一个问题就是计算顺序的问题。这个题就按顺序一行一行扫描就行了。

[cpp] view plaincopy
  1. #include<stdio.h>  
  2. #define XMAX 21  
  3. #define YMAX 21  
  4. int main() {  
  5.   __int64 matrix[XMAX][YMAX];  
  6.   for(int i=0;i<XMAX;i++) matrix[i][0]=1;  
  7.   for(int i=0;i<YMAX;i++) matrix[0][i]=1;  
  8.   for(int i=1;i<XMAX;i++)   
  9.     for(int j=1;j<YMAX;j++)  
  10.       matrix[i][j]=matrix[i-1][j]+matrix[i][j-1];  
  11.   printf("%lld",matrix[20][20]);  
  12. }  

原文地址:点击打开链接