bzoj1177 [Apio2009] Oil 分类讨论

题目大意：
给定n*m的矩阵，求三个互不相交的边长为k的子正方形的最大权值和。

题目分析：
因为这写矩形互不相交，所以我们可以把他们所在的区域划分成块。
YY一下，情况也就6种：（出题人真是良（sang）心（bing）啊）
1：先横着切一下，再把上面竖着切一下
2：先横着切一下，再把下面竖着切一下
3：先横着切一下，再把下面横着切一下
4：先竖着切一下，再把左面横着切一下
5：先竖着切一下，再把右面横着切一下
6：先竖着切一下，再把右面竖着切一下
（我写的代码和这个顺序不太一样QWQ）
这样的话我们只要求出每一大块内最大的边长为k的正方形的权值就可以了。
我们发现分出的块都包含大矩形的四个角之一，我们就分别预处理出（i，j）到左上，左下，右上，右下范围内最大的正方形，然后三块求和与ans取max就行了。
对于3和6这两种情况，我们可以枚举中间的那个正方形的位置，然后加上两边的。
这样再预处理出每一行最大的正方形和每一列最大的正方形即可。

（我觉的这个代码我写的十分整齐，但是这并不能遮盖它奇丑无比这个事实=。=）
代码如下：

#include <cstdio>
#define N 1520
using namespace std;
inline int Max(int x,int y) { return x>y?x:y; }
int n,m,k,ans;
int zs[N][N],ys[N][N],zx[N][N],yx[N][N];
int h[N],l[N];
int a[N][N];
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",&a[i][j]);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            a[i][j]+=a[i-1][j];
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            a[i][j]+=a[i][j-1];
    for(int i=n;i>=k;i--)
        for(int j=m;j>=k;j--)
            a[i][j]+=a[i-k][j-k]-a[i][j-k]-a[i-k][j];
    for(int i=k;i<=n;i++)
        for(int j=k;j<=m;j++)
            zs[i][j]=Max(Max(zs[i][j-1],zs[i-1][j]),a[i][j]);
    for(int i=k;i<=n;i++)
        for(int j=m-k+1;j>=1;j--)
            ys[i][j]=Max(Max(ys[i][j+1],ys[i-1][j]),a[i][j+k-1]);
    for(int i=n-k+1;i>=1;i--)
        for(int j=k;j<=m;j++)
            zx[i][j]=Max(Max(zx[i][j-1],zx[i+1][j]),a[i+k-1][j]);
    for(int i=n-k+1;i>=1;i--)
        for(int j=m-k+1;j>=1;j--)
            yx[i][j]=Max(Max(yx[i+1][j],yx[i][j+1]),a[i+k-1][j+k-1]);
    for(int i=k;i<=n;i++)
        for(int j=k;j<=m;j++)
            h[i]=Max(h[i],a[i][j]);
    for(int i=k;i<=n;i++)
        for(int j=k;j<=m;j++)
            l[j]=Max(l[j],a[i][j]);
    for(int i=k;i<=n-k-k;i++)
        ans=Max(ans,zs[i][m]+h[i+k]+zx[i+k+1][m]);
    for(int j=k;j<=m-k-k;j++)
        ans=Max(ans,zx[1][j]+l[j+k]+yx[1][j+k+1]);
    for(int i=k;i<=n-k;i++)
        for(int j=k;j<=m-k;j++)
            ans=Max(ans,zs[i][j]+ys[i][j+1]+yx[i+1][1]);
    for(int i=k;i<=n-k;i++)
        for(int j=k;j<=m-k;j++)
            ans=Max(ans,ys[i][1]+zx[i+1][j]+yx[i+1][j+1]);
    for(int j=k;j<=m-k;j++)
        for(int i=k;i<=n-k;i++)
            ans=Max(ans,zs[i][j]+zx[i+1][j]+yx[1][j+1]);
    for(int j=k;j<=m-k;j++)
        for(int i=k;i<=n-k;i++)
            ans=Max(ans,zx[1][j]+ys[i][j+1]+yx[i+1][j+1]);
    printf("%d\n",ans);
    return 0;
}