（C语言）最大子矩阵（乘积）

示例：

原理：

这个算法实现起来不算太复杂，具体解法如下：

首先将矩阵存入到一个二维数组里面。由于所求的是二阶子阵，所以只需要记录二阶子阵的左上角数据的坐标（下标）即可，利用双层循环，追踪左上角数据的坐标，并且将乘积赋值给 now，判断now与max的大小关系，如果now>max，则将该数据的坐标 (x,y) 存到一个数组sum[2]里面。（随时追踪并更新乘积和坐标）
最后将max，最大子矩阵按具体要求输出（因为sum[2]记录了最大子矩阵的左上角数据的坐标，很容易输出最大子矩阵）

话不多说，上代码！！！

代码如下：

#include<stdio.h>
int main(){
    int m,n,sum[2];
    scanf("%d%d",&m,&n);
    long int a[m][n];//存储矩阵信息
    long long int now=0,max=0;
    //由于每个数据最大为10000，需要使用long long int
    for(int i=0;i<m;i++){
        for(int j=0;j<n;j++){
            scanf("%ld",&a[i][j]);
        }
    }
    for(int i=0;i<m-1;i++){
        for(int j=0;j<n-1;j++){
            now=a[i][j]*a[i][j+1]*a[i+1][j]*a[i+1][j+1];
            if(now>max){
                max=now;
                sum[0]=i;
                sum[1]=j;
            }
            else{
                continue;
            }
        }
    }
    printf("%lld\n",max);
    int x=sum[0];
    int y=sum[1];
    printf("%ld %ld\n",a[x][y],a[x][y+1]);
    printf("%ld %ld",a[x+1][y],a[x+1][y+1]);
    return 0;
}

亿点细节：

在定义二维数组a[m][n]时,需要用 long int，如果用 int 定义，会发生溢出，具体原因如下:

long int 定义的变量，在被赋值时，不能使用过多的int 变量进行连乘，可能会产生溢出。
例如本算法中的这行代码：

now=a[i][j]*a[i][j+1]*a[i+1][j]*a[i+1][j+1];

尽管每个因子都没有超出int的范围，但多个因子连乘的结果可能会超出int 能够表示的范围，造成溢出。（其中int所能表示的最大值为1874919424 [10位]）
如果还不太明白，可以自行搜索相关知识，小编不再赘述。

上述代码是小编自己写的，如有雷同，还请见谅。