书法家（c++）

这是一道DP题（不要问我为什么），我们应该如何考虑呢？

根据房学长的原话讲：

此题属于根据数据范围猜做法系列

首先感受一下题意，像个DP

然后观察数据范围，感觉因该是个三次方的DP

然后发现n比m小不少，应该是O(n^2 * m)的。。。

（此处省略一万句吐槽）

言归正传，观察一下这个题的图：

然后我们可以发现一个非常奇妙的性质，

我们可以把这个图分成十一部分，按列枚举，可以发现，第十一部分可以由第十部分转移而来，也可以由自己转移而来（因为第十一部分的长度不限），第十部分同理，其他同理，注意第七部分只能由第六部分转移而来（第七部分只能有一列），第五部分只能由第四部分转移而来（同理），这样我们就可以尝试写DP了：

用fi表示第i块，枚举列进行转移，大概有以下转移：

empty->1,1->1,1->2,2->2,2->3,3->3,3->4,4->4,4->5,5->6,6->6,6->7,7->8,8->8,8->9,9

->9,9->10,10->10,10->11,11->11,11

最后还是房学长的话：大力DP就好啦！

附上代码：

#include <iostream>    
#include <cstdio>    
#include <cstring>    
#define clear(x) memset(x,-63,sizeof(x))    
using namespace std;    
const int maxn=155;    
const int maxm=505;    
int n,m,f4,f8,ans,INF;    
int a[maxm][maxn],s[maxn];    
int f1[maxn][maxn],f2[maxn][maxn],f3[maxn][maxn], f5[maxn][maxn],f6[maxn][maxn],f7[maxn][maxn],f9[maxn][maxn],f10[maxn][maxn],f11[maxn][maxn],s1[maxn][maxn],s2[maxn][maxn];    
inline int read(void){    
    int x=0,f=1; char ch=getchar();    
    for(;ch<'0'||ch>'9';ch=getchar()) if(ch=='-') f=-1;    
    for(;ch>='0'&&ch<='9';x=(x<<3)+(x<<1)+(ch^48),ch=getchar());    
    return x*f;     
}    
int main(int argc, char const *argv[]){    
    //初始化    
    n=read(),m=read();    
    for(int i=1;i<=n;i++)    
        for(int j=1;j<=m;j++)    
            a[j][n-i+1]=read();    
    clear(f1); clear(f2); clear(f3); clear(f5);    
    clear(f6); clear(f7); clear(f9); clear(f10);    
    clear(f11); clear(s1); clear(s2);    
    INF=-f1[0][0]; f4=f8=ans=-INF;    
    //DP    
    for(int j=1;j<=m;j++){    
        //前缀和    
        for(int i=1;i<=n;i++) s[i]=s[i-1]+a[j][i];    
        //第十一部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+2;r<=n;r++)    
                ans=max(ans,f11[l][r]=max(f11[l][r],f10[l][r])+a[j][l]+a[j][r]);    
        //第十部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+2;r<=n;r++)    
                f10[l][r]=max(f10[l][r],f9[l][r])+s[r]-s[l-1];    
        //第九部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+2;r<=n;r++)    
                f9[l][r]=max(f9[l][r],f8)+a[j][l]+a[j][r];    
        //第八部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+2;r<=n;r++)    
                f8=max(f8,f7[l][r]);    
        //第七部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+2;r<=n;r++)    
                f7[l][r]=f6[l][r]+s[r]-s[l-1];    
        //第六部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+2;r<=n;r++)    
                f6[l][r]=max(f6[l][r],f5[l][r])+a[j][l]+a[j][r];    
        //第五部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+2;r<=n;r++)    
                f5[l][r]=f4+s[r]-s[l-1];    
        //第四部分    
        for(int l=1;l<=n;l++)    
            for(int r=l+1;r<=n;r++)    
                f4=max(f4,f3[l][r]);    
        //这里有点不一样，自行体会    
        for(int l=1;l<=n;l++){    
            int tmp=-INF;    
            for(int r=l+1;r<=n;r++)    
                tmp=max(tmp,f2[l][r-1]),f3[l][r]=max(f3[l][r],tmp)+s[r]-s[l-1];    
        }    
        //同理,但要反过来处理    
        for(int r=1;r<=n;r++){    
            int tmp=s2[r+1][r];    
            for(int l=r;l;l--)    
                tmp=max(tmp,s2[l][r]),f2[l][r]=max(s1[l-1][r],tmp)+s[r]-s[l-1];    
        }    
        //第一部分    
        for(int l=1;l<=n;l++)    
            for(int r=l;r<=n;r++)    
                f1[l][r]=max(0,f1[l][r])+s[r]-s[l-1];    
        //处理前缀和    
        for(int l=1;l<=n;l++)    
            for(int r=n;r;r--)    
                s2[l][r]=max(f2[l][r],s2[l][r+1]);    
        for(int r=1;r<=n;r++)    
            for(int l=1;l<=n;l++)    
                s1[l][r]=max(f1[l][r],s1[l-1][r]);    
    }    
    cout<<ans<<endl;    
    return 0;    
}

 原题传送带