[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 ;

 }