P1398 [NOI2013]书法家

传送门

就是个普及组 $dp$ 合集，把 $NOI$ 从左到右拆成 $9$ 个部分，每个部分都可以分别 $dp$

除了 $N$ 的中间部分比较恶心以外其他都还好，自己推一下然后就知道转移，就 $N$ 的中间优化转移比较不好写

随便吧，反正 $9$ 个 $dp$ 都挺简单的，量变导致质变，我在想那一年的选手是不是都被恶心到了...反正我是被恶心死了

$luogu$ 上这一题空间限制太小了，原题是 $512MB$ ，所以这份代码在 $luogu$ 过不去（我都滚动数组了啊...）

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=,f=; char ch=getchar();
    while(ch<''||ch>'') { if(ch=='-') f=-; ch=getchar(); }
    while(ch>=''&&ch<='') { x=(x<<)+(x<<)+(ch^); ch=getchar(); }
    return x*f;
}
const int N=,M=;
int n,m,sum[N][M],ans;
int f[][M][N][N],mx[M][N][N],mxx,mxxx[N][N];
inline int calc(int x,int l,int r) { return sum[r][x]-sum[r][x-]-sum[l-][x]+sum[l-][x-]; }
int main()
{
    n=read(),m=read();
    for(int i=;i<=n;i++)
        for(int j=;j<=m;j++) sum[i][j]=sum[i-][j]+sum[i][j-]-sum[i-][j-]+read();
    memset(f,~0x3f,sizeof(f)); memset(mx,~0x3f,sizeof(mx));
    for(int j=;j<=n;j++)
        for(int k=j;k<=n;k++) f[][][j][k]=;
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++)
            for(int k=j;k<=n;k++)
                f[][i][j][k]=max(,f[][i-][j][k])+calc(i,j,k);
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++)
            for(int k=n-;k>=j;k--)
                mx[i][j][k]=max(mx[i][j][k+],f[][i-][j][k+]);
    for(int i=;i<=m;i++)
    {
        memset(mxxx,~0x3f,sizeof(mxxx));
        for(int k=n;k>=;k--)
            for(int j=;j<=k;j++) mxxx[j][k]=max(mxxx[j-][k],f[][i-][j][k]);
        for(int j=;j<=n;j++)
        {
            mxx=mxxx[j-][j-];
            for(int k=j;k<=n;k++)
            {
                mxx=max(mxx,mxxx[j][k]);
                f[][i][j][k]=max(mx[i][j][k],mxx)+calc(i,j,k);
            }
        }
    }
    memset(mx,~0x3f,sizeof(mx)); memset(f[],~0x3f,sizeof(f[]));
    for(int i=;i<=m;i++)
        for(int k=;k<=n;k++)
            for(int j=k-;j>=;j--)
                mx[i][j][k]=max(mx[i][j+][k],f[][i-][j+][k]);
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++)
            for(int k=j;k<=n;k++)
                f[][i][j][k]=max(mx[i][j][k],f[][i-][j][k])+calc(i,j,k);
    mxx=f[][][][]; memset(f[],~0x3f,sizeof(f[]));
    for(int i=;i<=m;i++)
    {
        for(int j=;j<=n;j++)
            for(int k=j+;k<=n;k++)
                f[][i][j][k]=mxx+calc(i,j,k);
        for(int j=;j<=n;j++)
            for(int k=j;k<=n;k++)
                mxx=max(mxx,f[][i-][j][k]);
    }
    memset(f[],~0x3f,sizeof(f[]));
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++)
            for(int k=j+;k<=n;k++)
                f[][i][j][k]=max(f[][i-][j][k],f[][i-][j][k])+calc(i,j,k)-calc(i,j+,k-);
    memset(f[],~0x3f,sizeof(f[]));
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++)
            for(int k=j+;k<=n;k++)
                f[][i][j][k]=f[][i-][j][k]+calc(i,j,k);
    mxx=f[][][][]; memset(f[],~0x3f,sizeof(f[]));
    for(int i=;i<=m;i++)
    {
        for(int j=;j<=n;j++)
            for(int k=j+;k<=n;k++)
                f[][i][j][k]=max(mxx,f[][i-][j][k])+calc(i,j,k)-calc(i,j+,k-);
        for(int j=;j<=n;j++)
            for(int k=j;k<=n;k++)
                mxx=max(mxx,f[][i-][j][k]);
    }
    memset(f[],~0x3f,sizeof(f[]));
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++)
            for(int k=j+;k<=n;k++)
                f[][i][j][k]=max(f[][i-][j][k],f[][i-][j][k])+calc(i,j,k);
    ans=f[][][][]; memset(f[],~0x3f,sizeof(f[]));
    for(int i=;i<=m;i++)
        for(int j=;j<=n;j++)
            for(int k=j+;k<=n;k++)
            {
                f[][i][j][k]=max(f[][i-][j][k],f[][i-][j][k])+calc(i,j,k)-calc(i,j+,k-);
                ans=max(ans,f[][i][j][k]);
            }
    printf("%d\n",ans);
    return ;
}