最大子矩阵问题实例解析_C 语言

问题:
求一个M*N的矩阵的最大子矩阵和。
比如在如下这个矩阵中:

 0 -2 -7 0
 9 2 -6 2
-4 1 -4 1
-1 8 0 -2

拥有最大和的子矩阵为:

 9 2
-4 1
-1 8

其和为15。
思路:
首先,这个子矩阵可以是任意大小的,而且起始点也可以在任何地方,所以,要把最大子矩阵找出来,我们要考虑多种情况。
假定原始矩阵的行数为M,那么对于子矩阵,它的行数可以是1到M的任何一个数,而且,对于一个K行(K < M)的子矩阵,它的第一行可以是原始矩阵的第1行到 M - K + 1 的任意一行。
例子:
对于上面的矩阵,如果子矩阵的行数是2,那么它可以是下面几个矩阵的子矩阵:

 0 -2 -7 0
 9 2 -6 2

或者

 9 2 -6 2
-4 1 -4 1

或者

-4 1 -4 1
-1 8 0 -2

在每一种情况里(我们这里有三种),我们还要找出一个最大的子矩阵,当然,这只是一种情况的最大子矩阵(局部最大),不一定是global最大。但是,如果我们知道每一种情况的最大,要找出global最大,那就小菜一碟儿了。
在讲在一个特殊情况下求最大子矩阵之前,先讲一个事实:
假设这个最大子矩阵的维数是一维,要找出最大子矩阵, 原理与求“最大子段和问题” 是一样的。最大子段和问题的递推公式是 b[j]=max{b[j-1]+a[j], a[j]},b[j] 指的是从0开始到j的最大子段和。

Java实现示例:
假设原始矩阵为:[9,  2, -6,  2], 那么b[] = {9, 11, 5, 7}, 那么最大字段和为11, 如果找最大子矩阵的话,那么这个子矩阵是 [9, 2]
求最大子段和的代码如下:

public int maxSubsequence(int[] array) {
 if (array.length == 0) {
 return 0;
 }
 int max = Integer.MIN_VALUE;
 int[] maxSub = new int[array.length];
 maxSub[0] = array[0];

 for (int i = 1; i < array.length; i++) {
 maxSub[i] = (maxSub[i-1] > 0) ? (maxSub[i-1] + array[i]) : array[i];
 if (max < maxSub[i]) {
 max = maxSub[i];
 }
 }
 return max;
 }

 但是,原始矩阵可以是二维的。假设原始矩阵是一个3 * n 的矩阵,那么它的子矩阵可以是 1 * k, 2 * k, 3 * k,(1 <= k <= n)。 如果是1*K,这里有3种情况:子矩阵在第一行,子矩阵在第二行,子矩阵在第三行。如果是 2 * k,这里有两种情况,子矩阵在第一、二行,子矩阵在第二、三行。如果是3 * k,只有一种情况。
为了能够找出最大的子矩阵,我们需要考虑所有的情况。假设这个子矩阵是 2 *k, 也就是说它只有两行,要找出最大子矩阵,我们要从左到右不断的遍历才能找出在这种情况下的最大子矩阵。如果我们把这两行上下相加,情况就和求“最大子段和问题” 又是一样的了。
为了找出在原始矩阵里的最大子矩阵,我们要遍历所有的子矩阵的可能情况,也就是说,我们要考虑这个子矩阵有可能只有1行,2行,。。。到n行。而在每一种情况下,我们都要把它所对应的矩阵部分上下相加才求最大子矩阵(局部)。
比如,假设子矩阵是一个3*k的矩阵,而且,它的一行是原始矩阵的第二行,那么,我们就要在

 9 2 -6 2
-4 1 -4 1
-1 8 0 -2

里找最大的子矩阵。
如果把它上下相加,我们就变成了 4, 11, -10,1, 从这个数列里可以看出,在这种情况下,最大子矩阵是一个3*2的矩阵,最大和是15.
为了能够在原始矩阵里很快得到从 i 行到 j 行 的上下值之和,我们这里用到了一个辅助矩阵,它是原始矩阵从上到下加下来的。
假设原始矩阵是matrix, 它每一层上下相加后得到的矩阵是total,那么我们可以通过如下代码实现:

int[][] total = matrix;
for (int i = 1; i < matrix[0].length; i++) {
 for (int j = 0; j < matrix.length; j++) {
 total[i][j] += total[i-1][j];
 }
}

如果我们要求第 i 行到第 j 行之间上下值的和,我们可以通过total[j][k] - total[i-1][k] 得到, k 的范围从1 到 matrix[0].length - 1。
有了这些知识点,我们只需要在所有的情况下,把它们所对应的局部最大子矩阵进行比较,就可以得到全局最大的子矩阵。代码如下:

public int subMaxMatrix(int[][] matrix) {

 int[][] total = matrix;
 for (int i = 1; i < matrix[0].length; i++) {
 for (int j = 0; j < matrix.length; j++) {
 total[i][j] += total[i-1][j];
 }
 }

 int maximum = Integer.MIN_VALUE;
 for (int i = 0; i < matrix.length; i++) {
 for (int j = i; j < matrix.length; j++) {
 //result 保存的是从 i 行 到第 j 行 所对应的矩阵上下值的和
        int[] result = new int[matrix[0].length];
 for (int f = 0; f < matrix[0].length; f++) {
  if (i == 0) {
  result[f] = total[j][f];
  } else {
  result[f] = total[j][f] - total[i - 1][f];
  }
 }
 int maximal = maxSubsequence(result);

 if (maximal > maximum) {
  maximum = maximal;
 }
 }
 }

 return maximum;
 }

C语言相关的实现

题目

    题目描述: 
    已知矩阵的大小定义为矩阵中所有元素的和。给定一个矩阵,你的任务是找到最大的非空(大小至少是1 * 1)子矩阵。 
    比如,如下4 * 4的矩阵  
      

 0 -2 -7 0
 9 2 -6 2
 -4 1 -4 1
 -1 8 0 -2

      
    的最大子矩阵是  
      

 9 2
 -4 1
 -1 8

      
    这个子矩阵的大小是15。 
    输入: 
    输入是一个N * N的矩阵。输入的第一行给出N (0 < N <= 100)。 
    再后面的若干行中,依次(首先从左到右给出第一行的N个整数,再从左到右给出第二行的N个整数……)给出矩阵中的N2个整数,整数之间由空白字符分隔(空格或者空行)。 
    已知矩阵中整数的范围都在[-127, 127]。 
    输出: 
    测试数据可能有多组,对于每组测试数据,输出最大子矩阵的大小。 
    样例输入:  
    4 
    0 -2 -7 0 
    9 2 -6 2 
    -4 1 -4  1 
    -1 8  0 -2  
    样例输出: 
    15 

AC代码

 #include <stdio.h>
 #include <stdlib.h> 

 int main(void)
 {
  int i, j, h, k, n, max, sum, cur, matrix[101][101]; 

  while (scanf("%d", &n) != EOF) {
   // 初始化接收矩阵
   for (i = 0; i < n; i ++) {
    for (j = 0; j < n; j ++)
     scanf("%d", *(matrix + i) + j);
   } 

   // 动态规划(类似于一维数组连续最大子序列和)
   max = matrix[0][0]; 

   for (i = 0; i < n; i ++) {
    // i,j确定上下界
    for (j = i; j < n; j ++) {
     // 初始化
     for (k = i, sum = 0; k <= j; k ++)
      sum += matrix[k][0];
     if (sum > max)
      max = sum; 

     for (h = 1; h < n; h ++) {
      for (k = i, cur = 0; k <= j; k ++)
       cur += matrix[k][h]; 

      if (sum >= 0)
       sum += cur;
      else
       sum = cur; 

      if (sum > max) max = sum;
     }
    }
   } 

   printf("%d\n", max);
  } 

  return 0;
 } 

以上是小编为您精心准备的的内容,在的博客、问答、公众号、人物、课程等栏目也有的相关内容,欢迎继续使用右上角搜索按钮进行搜索java
, c
最大子矩阵
c语言高级实例解析、c语言实例解析精粹pdf、c语言高级实例解析pdf、c语言实例解析精粹、易语言 解析html实例,以便于您获取更多的相关知识。

时间: 2024-10-28 11:52:06

最大子矩阵问题实例解析_C 语言的相关文章

C++普通函数指针与成员函数指针实例解析_C 语言

C++的函数指针(function pointer)是通过指向函数的指针间接调用函数.相信很多人对指向一般函数的函数指针使用的比较多,而对指向类成员函数的函数指针则比较陌生.本文即对C++普通函数指针与成员函数指针进行实例解析. 一.普通函数指针 通常我们所说的函数指针指的是指向一般普通函数的指针.和其他指针一样,函数指针指向某种特定类型,所有被同一指针运用的函数必须具有相同的形参类型和返回类型. int (*pf)(int, int); // 声明函数指针 这里,pf指向的函数类型是int (

C++实现判断字符串是否回文实例解析_C 语言

本文实例解析了C++判断字符串是否回文的实现过程,通过数据结构中的相关例子,回文判断中采用过滤空格字符.有效字符依次入栈等方法实现该功能. 具体实例代码如下: #include <iostream> using namespace std; #define Max_String_Len 100 #include "SqStack.h" //判断字符串是否回文 bool ispalindrome(char *in_string) { SqStack <char>

C++日志记录类实例解析_C 语言

本文所述实例是从一个Red Hat开源项目里面扒出来的,非常实用!读者还可以根据自身需求加以修改!完整源码如下: Log.h文件部分: #ifndef __LOG_H__ #define __LOG_H__ #include <stdio.h> #include <tchar.h> #include <crtdbg.h> #include <windows.h> #include <time.h> #include <sys/timeb.h

C++函数模板与类模板实例解析_C 语言

本文针对C++函数模板与类模板进行了较为详尽的实例解析,有助于帮助读者加深对C++函数模板与类模板的理解.具体内容如下: 泛型编程(Generic Programming)是一种编程范式,通过将类型参数化来实现在同一份代码上操作多种数据类型,泛型是一般化并可重复使用的意思.泛型编程最初诞生于C++中,目的是为了实现C++的STL(标准模板库). 模板(template)是泛型编程的基础,一个模板就是一个创建类或函数的蓝图或公式.例如,当使用一个vector这样的泛型类型或者find这样的泛型函数

C++中vector的用法实例解析_C 语言

本文实例展示了C++中的vector用法,分享给大家供大家参考.具体如下: 一.概述 vector是C++标准模板库中的部分内容,它是一个多功能的,能够操作多种数据结构和算法的模板类和函数库.vector是一个容器,它能够存放各种类型的对象,简单地说,vector是一个能够存放任意类型的动态数组,可以动态改变大小. 例如: // c语言风格 int myHouse[100] ; // 采用vector vector<int> vecMyHouse(100); 当如上定义后,vecMyHouse

C++命名空间实例解析_C 语言

命名空间是C++非常重要的概念,本文就以实例形式对其进行深入分析,具体内容如下: 通常来说,在C++中,命名空间(namespace)的目的是为了防止名字冲突.每个命名空间是一个作用域,在所有命名空间之外,还存在一个全局命名空间(global namespace),全局命名空间以隐式的方式声明,它并没有名字.在命名空间机制中,原来的全局变量,就是位于全局命名空间中(可以用::member的形式表示). 一.定义命名空间 1.每个命名空间都是一个作用域 和其他作用域类似,在命名空间中的每个名字必须

C语言安全之数组长度与指针实例解析_C 语言

1.C语言编码需要保证变长数组的长度参数位于合法范围之内 例如以下代码: void func(size_t s) { int vla[s]; /*...*/ } /*...*/ func(size); /*...*/ 解决方案如下: enum {MAX_ARRAY = 1024}; void func(size_t s) { if(s < MAX_ARRAY && s != 0) { int vla[s]; /*...*/ } else { //错误处理 } } /*...*/ fu

C++设计模式编程中Facade外观模式的使用实例解析_C 语言

外观模式提供了一个统一的接口,用来访问子系统的一群接口.外观定义了一个高层接口,让子系统更容易使用.外观模式让接口变得简单,简化了子系统的接口.外观模式十分简单,简而言之,就是简化你的类的接口,将一系列的复杂的过程封装到内部,对外只提供最简单的接口.结构图: 适用场景:当你要为一个复杂子系统提供一个简单接口时.子系统往往因为不断演化而变得越来越复杂.大多数模式使用时都会产生更多更小的类.这使得子系统更具可重用性,也更容易对子系统进行定制,但这也给那些不需要定制子系统的用户带来一些使用上的困难.F

C++形参与实参的区别实例解析_C 语言

本文以实例阐述了C++中形参与实参的区别,有助于读者加深对于C++形参与实参的认识. 形参出现在函数定义中,在整个函数体内都可以使用, 离开该函数则不能使用.实参出现在主调函数中,进入被调函数后,实参变量也不能使用. 形参和实参的功能是作数据传送.发生函数调用时, 主调函数把实参的值传送给被调函数的形参从而实现主调函数向被调函数的数据传送. 1.形参变量只有在被调用时才分配内存单元,在调用结束时, 即刻释放所分配的内存单元.因此,形参只有在函数内部有效. 函数调用结束返回主调函数后则不能再使用该