九度[1102]-最小面积子矩阵

九度[1102]-最小面积子矩阵

题目描述：
一个N*M的矩阵，找出这个矩阵中所有元素的和不小于K的面积最小的子矩阵（矩阵中元素个数为矩阵面积）

输入
每个案例第一行三个正整数N,M<=100，表示矩阵大小，和一个整数K
接下来N行，每行M个数，表示矩阵每个元素的值

输出
输出最小面积的值。如果出现任意矩阵的和都小于K，直接输出-1。

样例输入
4 4 10
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16

样例输出
1

解题思路：
暴力解法，利用缓存并剪枝来减少运行时间。
大神们还有更厉害的做法（最短连续子序列），详见链接

AC代码：

#include 
#include 
#include 
#define MAX 1000000000
const int maxn = 110;
int N, M, K;
int dat[maxn][maxn], sumRow[maxn][maxn];

int sum(int x1, int y1, int x2, int y2){
    if(x1 < x2) {
        int t = x1;
        x1 = x2;
        x2 = t;
    }
    if(y1 < y2) {
        int t = y1;
        y1 = y2;
        y2 = t;
    }
    int ans = 0;
    for(int i = x2; i <= x1; i++){
        ans += sumRow[i][y1] - sumRow[i][y2] + dat[i][y2];
    }
    return ans;
}

int bruce(){
    int min = MAX;
    for(int i = 0; i < N; i++){
        for(int j = 0; j < M; j++){
            for(int ii = 0; ii < N; ii++){
                for(int jj = 0; jj < M; jj++){
                    if((abs(i-ii)+1)*(abs(j-jj)+1) < min){
                        if(sum(i, j, ii, jj) >= K){
                            //printf("%d %d %d %d\n", i, j, ii, jj);
                            min = (abs(i-ii)+1)*(abs(j-jj)+1);
                        }
                    }
                }
            }
        }
    }
    if(min < MAX) return min;
    else return -1;
}

int main(){
    //freopen("C:\\Users\\Administrator\\Desktop\\test.txt", "r", stdin);
    while(scanf("%d%d%d", &N, &M, &K) != EOF){
        for(int i = 0; i < N; i++){
            for(int j = 0; j < M; j++){
                scanf("%d", &dat[i][j]);
                if(j > 0) sumRow[i][j] = sumRow[i][j-1] + dat[i][j];
                else sumRow[i][j] = dat[i][j];
            }
        }
        printf("%d\n", bruce());
    }
    //fclose(stdin);
    return 0;
} 

九度[1102]-算法：最短连续子序列