bzoj 3111 [Zjoi2013]蚂蚁寻路 二维dp

题意：告诉你蚂蚁的行走路线是 右转右转左转左转 -----  (最后三次转弯要求全是右转）题目告诉你K 蚂蚁左转的次数>>1
最后蚂蚁必须回到起点，问你蚂蚁的行走路线包含的值最大
解：
1、分析蚂蚁行走的路线发现蚂蚁走的路线 是一个类型 手掌一样的形状（n指）
2、转化为求k个底相同高低错落（大小大小---）的矩形的最大和
3、

 for(int i=1; i<=n; ++i)
        for(int j=1; j<=m; ++j) {
            for(int k=1; k<=K; ++k) {
                for(int h=i; h<=n; ++h) {
                    dp[j][k][h]=max(dp[j-1][k][h],g[j-1][k-1][h+1-2*(k&1)])+ss[j][h]-ss[j][i-1];
                }

（1）i 是底 J为列 K为矩形个数 h为最后这一列的高度
（2）k个矩形要么从k个转移or k-1转移  （k分奇偶    奇从h小的转移  偶从h大的转移）
（3）g是个优化 （我们每次转移的想从K-1转移 一定是  前缀里最大的 or 后缀最大的）（复杂度降一维）

（4）初始化：mem(dp,-inf) ，对于g 全设置成-inf是不合理的， 因为max()那一项无论怎样都会是-inf 是不合理的
所以我们让矩形个数 k为0时 g为0，所以当k==1时dp可以从（0 +sum[h]-sum[i-1]）转移过来了。
程序中注释部分是我最开始的写法 。

#include<bits/stdc++.h>
#define mod 51061
#define en '\n'
#define mem(a,b)  memset(a,b,sizeof(a))
#define max(a,b) (a>b? a:b)
#define ll unsigned int
const int maxn = 1e2 +5;
const int inf=0x3f3f3f3f;
using namespace std;
int read() {
    int x=0,f=1;
    char ch=getchar();
    while(ch<'0'||ch>'9') {
        if(ch=='-')
            f=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9') {
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}

int n,m,K;
int dp[maxn][maxn][maxn],g[maxn][maxn][maxn],s[maxn][maxn],ss[maxn][maxn];
signed main() {
#ifdef local
    freopen("input2.txt","r",stdin);
#endif // local
#define int register int
    n=read(),m=read();
    K=read();
    K=K<<1|1;
    for(int i=n; i>=1; --i) {
        for(int j=1; j<=m; j++) {
            s[i][j]=read();
        }
    }
    for(int j=1; j<=m; j++) {
        for(int i=1; i<=n; i++) {
            ss[j][i]=ss[j][i-1]+s[i][j];
        }
    }
    ///8888888888888888888888888888888888888888888888888
    int ans=-inf;
    mem(dp,-inf);
    for(int j=0;j<=m;++j)
        for(int k=1;k<=K;++k)
            for(int h=0;h<=n+1;++h)g[j][k][h]=-inf;

    //mem(g,-inf);
//    for(int h=1; h<=n; ++h) {
//        for(int j=0; j<=m; ++j)
//            g[j][0][h]=0;
//    }
    for(int i=1; i<=n; ++i)
        for(int j=1; j<=m; ++j) {
            for(int k=1; k<=K; ++k) {
                for(int h=i; h<=n; ++h) {
                    dp[j][k][h]=max(dp[j-1][k][h],g[j-1][k-1][h+1-2*(k&1)])+ss[j][h]-ss[j][i-1];
                }
                g[j][k][i-1]=-inf;
                if(k&1) {
                    for(int h=n; h>=i; --h)
                        g[j][k][h]=max(g[j][k][h+1],dp[j][k][h]);
                } else {
                    for(int h=i; h<=n; ++h)
                        g[j][k][h]=max(g[j][k][h-1],dp[j][k][h]);
                }
                if(k==K) {
                    for(int h=i; h<=n; ++h) {
                        ans=max(ans,dp[j][k][h]);
                    }
                }
            }
        }
    cout<<ans<<en;
    return 0;
}