最大子矩阵（多维dp）

#include<bits/stdc++.h>
using namespace std;
int a[110][110];
int dp[110][110][20];
int s[110][3];
int main()
{
    int n,m,t;
    cin>>n>>m>>t;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            cin>>a[i][j];
            s[i][j]=s[i-1][j]+a[i][j];
        }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            for(int k=1;k<=t;k++){
                dp[i][j][k]=max(dp[i-1][j][k],dp[i][j-1][k]);
                for(int p=0;p<i;p++)
                    dp[i][j][k]=max(dp[i][j][k],dp[p][j][k-1]+s[i][1]-s[p][1]);
                for(int p=0;p<j;p++)
                    dp[i][j][k]=max(dp[i][j][k],dp[i][p][k-1]+s[j][2]-s[p][2]);
                if(i==j)
                    for(int p=0;p<i;p++)
                        dp[i][j][k]=max(dp[i][j][k],dp[p][p][k-1]+s[i][1]+s[j][2]-s[p][1]-s[p][2]);
            }
    int ans=0;
    for(int i=1;i<=n;i++){
    	for(int j=1;j<=n;j++){
    		ans=max(ans,dp[i][j][t]);
		}
	}
	printf("%d\n",ans);
    return 0;
}