最大子矩阵问题（dp）

最大子矩阵问题（dp）

题目描述

在一个N*M的大矩阵中，求出数字和最大的子矩阵。

输入

第一行有两个数字N和M（1<=N,M<=100），接下来是一个N*M的矩阵。

输出

仅一个数即最大子矩阵的数字之和，数据保证答案在longint范围内。

样例输入 Copy

3 3
1 -5 4
-2 3 2
1 3 3 

样例输出 Copy

11 

提示

对于30%的数据1<=N,M<=10;
对于70%的数据1<=N,M<=50;

对于100%的数据1<=N,M<=100。

感觉就是最大字段和的翻版。二维形式，记录每一列的前缀和。动态的求最大值。

#pragma comment(linker , "/STACK:1024000000,1024000000")
#pragma GCC optimize(3 , "Ofast" , "inline")
 
#include <bits/stdc++.h>
 
using namespace std;
 
#define rep(i , a , b) for(register int i=(a);i<=(b);i++)
#define per(i , a , b) for(register int i=(a);i>=(b);i--)
#define ms(s) memset(s, 0, sizeof(s))
#define squ(x) (x)*(x)
 
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll , ll> pi;
typedef unordered_map<int , int> un_map;
 
template<class T>
inline void read (T &x) {
    x = 0;
    int sign = 1;
    char c = getchar ();
    while (c < '0' || c > '9') {
        if ( c == '-' ) sign = - 1;
        c = getchar ();
    }
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar ();
    }
    x = x * sign;
}
 
const int maxn = 1010;
const int inf = 0x3f3f3f3f;
const ll INF = ll (1e18);
const int mod = 1e9 + 7;
const double PI = acos (- 1);
 
 
int n,m;
ll mat[111][111];
ll dp[111][111];
 
int main()
{
    read(n);read(m);
    rep(i,1,n) {
        rep(j,1,m) {
            read(mat[i][j]);
            dp[i][j]=dp[i-1][j]+mat[i][j];
        }
    }
    ll ans = mat[1][1];
    rep(i,1,n) {
        rep(j,i,n) {
            ll tmp = 0;
            rep(k,1,m) {
                tmp+=dp[j][k]-dp[i-1][k];
                if(tmp<0) tmp=0;
                if(tmp>ans) ans=tmp;
            }
        }
    }
    cout<<ans<<endl;
    return 0;
}

求子矩阵为r*c的最大值。

题目描述

在幻想乡，八云紫是最像妖怪的妖怪。她的境符「四重结界」能够消除结界里面所有的弹幕。假设现在八云紫面前有一块N行M列的弹幕群，每一个单位面积内有分值有num[i][j]的弹幕。八云紫可以将一块R行C列的弹幕消除，并且得到这一块区域内所有弹幕的分值(累加)。现在八云紫想要取得尽可能多的分值，请你计算出她最多能够得到的分值。

输入

第1行：4个正整数N,M,R,C 第2..N+1行：每行M个正整数，第i+1行第j个数表示num[i][j]

输出

第1行：1个整数，表示八云紫能够取得的最大得分

样例输入 Copy

3 5 2 3
5 2 7 1 1
5 9 5 1 5
3 5 1 5 3 

样例输出 Copy

33 

提示

对于60%的数据：1 <= N,M <= 200
对于100%的数据：1 <= N,M <= 1,000 1 <= R <= N, 1 <= C <= M 1 <= num[i][j] <= 1000
保证结果不超过2,000,000,000

二维前缀和。

#pragma comment(linker, "/STACK:1024000000,1024000000")
#pragma GCC optimize(3,"Ofast","inline")
#include <bits/stdc++.h>
 
using namespace std;
 
#define rep(i , a , b) for(register int i=(a);i<=(b);i++)
#define per(i , a , b) for(register int i=(a);i>=(b);i--)
#define ms(s) memset(s, 0, sizeof(s))
#define squ(x) (x)*(x)
 
typedef long long ll;
typedef unsigned long long ull;
typedef pair<ll , ll> pi;
typedef unordered_map<int,int> un_map;
template<class T>
inline void read (T &x) {
    x = 0;
    int sign = 1;
    char c = getchar ();
    while (c < '0' || c > '9') {
        if ( c == '-' ) sign = - 1;
        c = getchar ();
    }
    while (c >= '0' && c <= '9') {
        x = x * 10 + c - '0';
        c = getchar ();
    }
    x = x * sign;
}
 
const int maxn = 1e6 + 10;
const int inf = 0x3f3f3f3f;
const ll INF = ll(1e18);
const int mod = 1e9+7;
const double PI = acos(-1);
//#define LOCAL
 
ll dp[1005][1005];
int n,m,r,c;
int main(int argc, char * argv[])
{
    #ifdef LOCAL
    freopen("D.in", "r", stdin);
    //freopen("/home/yesky/桌面/date.out", "w", stdout);
    #endif
 
 
    while (~scanf("%d%d%d%d",&n,&m,&r,&c))
    {
        ll ans = 0;
        rep(i,1,n) {
            rep(j,1,m) {
                read(dp[i][j]);
                dp[i][j]+=dp[i-1][j]+dp[i][j-1]-dp[i-1][j-1];
                if(i>=r&&j>=c) ans=max(ans,dp[i][j]-dp[i][j-c]-dp[i-r][j]+dp[i-r][j-c]);
            }
        }
        printf("%lld\n",ans);
    }
 
    return 0;
}