BZOJ 1057 [ZJOI2007]棋盘制作

最新推荐文章于 2019-02-18 12:09:00 发布

原创最新推荐文章于 2019-02-18 12:09:00 发布 · 349 阅读

1 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#BZOJ-1057

动态规划-矩阵专栏收录该内容

2 篇文章

订阅专栏

本文介绍了一种使用悬线法动态规划解决最大正方形及矩形问题的方法。通过预处理水平方向最长连续1的长度l[i][j]和垂直方向最长连续1的长度r[i][j]，再结合高度h[i][j]，实现了对最大正方形和矩形的有效求解。

悬线法DP。

要求01相邻，方便起见，对于格子(i,j)，若(i+j)为奇数，就对它取反（异或1），题目的要求就可以变成相邻11或00而不是01。

第一问最大正方形比较好想，记f[i][j]为以(i,j)为右下角的最大正方形，可由f[i-1][j-1]通过一定的姿势转移过来。

第二问用到悬线法（orz好腻害），资料传送：《浅谈用极大化思想解决最大子矩形问题》

#include<cstdio>
#include<algorithm>
#define min3(x,y,z) min(min(x,y),z)
#define N 2005
using namespace std;
int n, m, map[N][N], ans1=0, ans2=0, f[N][N], l[N][N], r[N][N], h[N][N];
void solve1()//正方形 
{
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            if(map[i][j])l[i][j]=l[i][j-1]+1;
            else l[i][j]=0;
    for(int j = 1; j <= m; j++)
        for(int i = 1; i <= n; i++)
            if(map[i][j])r[i][j]=r[i-1][j]+1;
            else r[i][j]=0;

    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
        {
            f[i][j]=min3(f[i-1][j-1]+1,l[i][j],r[i][j]);
            ans1=max(ans1,f[i][j]);
        }
}
void solve2()//矩形
{
    for(int j = 1; j <= m; j++)
        if(map[1][j])h[1][j]=1;
        else h[1][j]=0;
    for(int i = 1; i <= n; i++)
        for(int j = m; j >= 1; j--)
            if(map[i][j])r[i][j]=r[i][j+1]+1;
            else r[i][j]=0;
    for(int i = 2; i <= n; i++)
        for(int j = 1; j <= m; j++)
        {
            if(!map[i][j]){h[i][j]=0;continue;}
            h[i][j]=h[i-1][j]+1;
            if(h[i][j]>1)
                l[i][j]=min(l[i][j],l[i-1][j]),
                r[i][j]=min(r[i][j],r[i-1][j]);
        }
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            ans2=max(ans2,h[i][j]*(l[i][j]+r[i][j]-1));
} 
int main()
{
    scanf("%d%d",&n,&m);
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
        {
            scanf("%d",&map[i][j]);
            if((i+j)&1)map[i][j]^=1;
        }
    solve1();
    solve2();
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++)
            map[i][j]^=1;
    solve1();
    solve2();
    printf("%d\n%d\n",ans1*ans1,ans2);
    return 0;
}