[bzoj1177]Oil

考虑将三个矩形按某种方式划分为再三个大矩形中找最大值，容易发现只有6种划分方式，分为两类：
1.4种，考虑第一条横/竖和第二条在第一条的两侧，这一类情况只需要预处理出左上/左下/右上/右下的最大子矩阵即可
2.2种，两条横线/竖线划分，以横线为例，只需要处理处每一行的答案，然后，然后简单维护一下即可

#include<bits/stdc++.h>
using namespace std;
#define N 2005
int n,m,k,x,y,ans,f[N],a[N][N],ll[N][N],lr[N][N],rl[N][N],rr[N][N];
int S(int x,int y){
    return a[x][y]-a[x-k][y]-a[x][y-k]+a[x-k][y-k];
}
int main(){
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            scanf("%d",&x);
            a[i][j]=a[i][j-1]+a[i-1][j]-a[i-1][j-1]+x;
        }
    for(int i=k;i<=n;i++)
        for(int j=k;j<=m;j++)
            ll[i][j]=max(max(ll[i][j-1],ll[i-1][j]),S(i,j));
    for(int i=k;i<=n;i++)
        for(int j=m-k+1;j;j--)
            lr[i][j]=max(max(lr[i][j+1],lr[i-1][j]),S(i,j+k-1));
    for(int i=n-k+1;i;i--)
        for(int j=k;j<=m;j++)
            rl[i][j]=max(max(rl[i][j-1],rl[i+1][j]),S(i+k-1,j));
    for(int i=n-k+1;i;i--)
        for(int j=m-k+1;j;j--)
            rr[i][j]=max(max(rr[i][j+1],rr[i+1][j]),S(i+k-1,j+k-1));
    for(int i=k;i<=n-k;i++)
        for(int j=k;j<=m-k;j++){
            ans=max(ans,ll[n][j]+lr[i][j+1]+rr[i+1][j+1]);
            ans=max(ans,ll[i][j]+lr[i][j+1]+rr[i+1][1]);
            ans=max(ans,ll[i][m]+rl[i+1][j]+rr[i+1][j+1]);
            ans=max(ans,ll[i][j]+rl[i+1][j]+rr[1][j+1]);
        }
    for(int i=k;i<=n;i++)
        for(int j=k;j<=m;j++)f[i]=max(f[i],S(i,j));
    x=0;
    for(int i=k;i<=n;i++){
        x=max(x,f[i]);
        int y=0;
        for(int j=i+k;j<=n-k;j++){
            y=max(y,f[j]);
            ans=max(ans,x+y+rr[j+1][1]);
        }
    }
    x=0;
    memset(f,0,sizeof(f));
    for(int i=k;i<=n;i++)
        for(int j=k;j<=m;j++)f[j]=max(f[j],S(i,j));
    for(int i=k;i<=m;i++){
        x=max(x,f[i]);
        int y=0;
        for(int j=i+k;j<=m-k;j++){
            y=max(y,f[j]);
            ans=max(ans,x+y+rr[1][j+1]);
        }
    }
    printf("%d",ans);
}