[POJ1050]To the Max（最大子矩阵，DP）

题目链接：http://poj.org/problem?id=1050

发现这个题没有写过题解，现在补上吧，思路挺经典的。

思路就是枚举所有的连续的连续的行，比如1 2 3 4 12 23 34 45 123 234 345...然后把这些行对应列相加缩成一行，之后就是求最大子序列和了。

     /*
 ━━━━━┒ギリギリ♂ eye！
 ┓┏┓┏┓┃キリキリ♂ mind！
 ┛┗┛┗┛┃＼○／
 ┓┏┓┏┓┃ /
 ┛┗┛┗┛┃ノ)
 ┓┏┓┏┓┃
 ┛┗┛┗┛┃
 ┓┏┓┏┓┃
 ┛┗┛┗┛┃
 ┓┏┓┏┓┃
 ┛┗┛┗┛┃
 ┓┏┓┏┓┃
 ┃┃┃┃┃┃
 ┻┻┻┻┻┻
 */
 #include <algorithm>
 #include <iostream>
 #include <iomanip>
 #include <cstring>
 #include <climits>
 #include <complex>
 #include <fstream>
 #include <cassert>
 #include <cstdio>
 #include <bitset>
 #include <vector>
 #include <deque>
 #include <queue>
 #include <stack>
 #include <ctime>
 #include <set>
 #include <map>
 #include <cmath>
 using namespace std;
 #define fr first
 #define sc second
 #define cl clear
 #define BUG puts("here!!!")
 #define W(a) while(a--)
 #define pb(a) push_back(a)
 #define Rint(a) scanf("%d", &a)
 #define Rll(a) scanf("%lld", &a)
 #define Rs(a) scanf("%s", a)
 #define Cin(a) cin >> a
 #define FRead() freopen("in", "r", stdin)
 #define FWrite() freopen("out", "w", stdout)
 #define Rep(i, len) for(int i = 0; i < (len); i++)
 #define For(i, a, len) for(int i = (a); i < (len); i++)
 #define Cls(a) memset((a), 0, sizeof(a))
 #define Clr(a, x) memset((a), (x), sizeof(a))
 #define Full(a) memset((a), 0x7f7f7f, sizeof(a))
 #define lrt rt << 1
 #define rrt rt << 1 | 1
 #define pi 3.14159265359
 #define RT return
 #define lowbit(x) x & (-x)
 #define onenum(x) __builtin_popcount(x)
 typedef long long LL;
 typedef long double LD;
 typedef unsigned long long ULL;
 typedef pair<int, int> pii;
 typedef pair<string, int> psi;
 typedef pair<LL, LL> pll;
 typedef map<string, int> msi;
 typedef vector<int> vi;
 typedef vector<LL> vl;
 typedef vector<vl> vvl;
 typedef vector<bool> vb;

 const int maxn = ;
 int G[maxn][maxn];
 int n;
 int dp[maxn];

 int lss() {
     int tmp = -0x7f7f7f, ret = -0x7f7f7f;
     For(i, , n+) {
         if(tmp > ) tmp += dp[i];
         else tmp = dp[i];
         ret = max(tmp, ret);
     }
     return ret;
 }

 int main() {
     // FRead();
     while(~Rint(n)) {
         For(i, , n+) For(j, , n+) Rint(G[i][j]);
         int ret = ;
         For(i, , n+) {
             Cls(dp);
             For(j, i, n+) {
                 For(k, , n+) dp[k] += G[j][k];
                 ret = max(ret, lss());
             }
         }
         printf("%d\n", ret);
     }
     RT ;
 }