dp 最大全1子正方形

例题1：usaco Big Barn 巨大的牛棚

Big Barn 巨大的牛棚

Description

农夫约翰想要在他的正方形农场上建造一座正方形大牛棚。他讨厌在他的农场中砍树，想找一个能够让他在空旷无树的地方修建牛棚的地方。我们假定，他的农场划分成 N x N 的方格。输入数据中包括有树的方格的列表。你的任务是计算并输出，在他的农场中，不需要砍树却能够修建的最大正方形牛棚。牛棚的边必须和水平轴或者垂直轴平行。

考虑下面的方格，它表示农夫约翰的农场，‘.'表示没有树的方格，‘#'表示有树的方格

最大的牛棚是 5 x 5 的，可以建造在方格右下角的两个位置其中一个。

Input

Line 1: 两个整数： N （1 <= N <= 1000），农场的大小，和 T （1 <= T <= 10，000）有树的方格的数量

Lines 2..T+1: 两个整数（1 <= 整数 <= N), 有树格子的横纵坐标

Output

输出文件只由一行组成，约翰的牛棚的最大边长。

Sample Input

8 3

2 2

2 6

6 3

Sample Output

5

如图，此题不难看出结果是5，但是dp方程是什么呢？

这是个2维的图形，很显然一维的状态是不够的，设f[i][j]表示以i行j列为右下角(左上角)的最大正方形边长。

这个最大正方形与前边3个有关。

f[i][j]=min(f[i][j-1],f[i-1][j],f[i-1][j-1])+1;

代码如下:

#include<stdio.h>
int f[1001][1001];
int main()
{
	int n,shu,i,o,p,j,ans=0;
	scanf("%d%d",&n,&shu);
	for(i=1;i<=n;i++)
		for(j=1;j<=n;j++)
			f[i][j]=1;
	for(i=1;i<=shu;i++)
	{
		scanf("%d%d",&o,&p);
		f[o][p]=0;		
	}
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
		{
			if(f[i][j])
			{
				if(f[i-1][j]>f[i][j-1])
					f[i][j]=f[i][j-1]+1;
				else
					f[i][j]=f[i-1][j]+1;
			}
			if(f[i-1][j-1]+1<f[i][j])
				f[i][j]=f[i-1][j-1]+1;
		}
	}
	for(i=1;i<=n;i++)
	{
		for(j=1;j<=n;j++)
		{
			if(f[i][j]>ans)
				ans=f[i][j];	
		}
	}
	printf("%d",ans);
}

例2：Home on the Range 家的范围

Home on the Range 家的范围

Description

农民约翰在一片边长是N (2 <= N <= 250)英里的正方形牧场上放牧他的奶牛.(因为一些原因,他的奶牛只在正方形的牧场上吃草.)
遗憾的是,他的奶牛已经毁坏一些土地.( 一些1 平方英里的正方形)
农民约翰需要统计那些可以放牧奶牛的正方形牧场(至少是2x2 的,在这些较大的正方形中没有小于1x1 的部分被分割毁坏).
你的工作要在被供应的数据组里面统计所有不同的正方形放牧区域(>2x2)的个数.
当然,放牧区域可能是重叠.

Input

第 1 行: N,牧区的边长.
第 2 到n+1 行: N 个没有空格分开的字符.
0 表示 "那一个区段被毁坏了";1 表示 " 准备好被吃".

Output

输出那些存在的正方形的大小和个数,一种一行.

Sample Input

6

101111

001111

111111

001111

1 01101

111001

Sample Output

2 10

3 4

4 1

此题与上题有一样了，不解释。

#include<stdio.h>
int a[251][251],f[251][251],ans[251];
int main()
{
	int n,i,j,k,x;
	scanf("%d",&n);
	for(i=1;i<=n;i++)
		for(j=1;j<=n;j++)
			scanf("%1d",&a[i][j]);
    for(i=n-1;i>=1;i--)
    {
        for(j=n-1;j>=1;j--)
        {
            if(a[i][j])
            {
            	if(a[i+1][j]>a[i][j+1])
            		x=a[i][j+1]+1;
           		else
           			x=a[i+1][j]+1;
          		if(a[i+1][j+1]+1<x)
				  	x=a[i+1][j+1]+1;
				a[i][j]=x;	
            }
            if(a[i][j]>1)
                for(k=2;k<=a[i][j];k++)
                    ans[k]++;
        }
    }
    for(i=1;i<=251;i++)
    {
    	if(ans[i])
    		printf("%d %d\n",i,ans[i]);
    }
}

例3：最大正方形

Description

背景描述：忙完了学校的事，v神终于可以做他的“正事”：陪女朋友散步。一天，他和女朋友走着走着，不知不觉就来到了一个千里无烟的地方。v神正要往回走，如发现了一块牌子，牌子上有有一行小字和一张图，小字说道：“找到图上最大的交错正方形之后和我联系，这块地就是你的了。”在房价疯长的年代，v神当然不愿错过这个机会，于是开始找了起来……以v神的能力当然找不出来了，你能帮v神找出来吗？

题目描述：

图上有一个矩阵，由N*M个格子组成，这些格子由两种颜色构成，黑色和白色。请找到面积最大的且内部是黑白交错（即两个相连的正方形颜色不能相同）的正方形。

样例解释：

(1,1)到(2,2)这个正方形是满足条件的，它的边长是2

 

数据范围约定：

对于30%的数据，N <= 20

对于60%的数据，N <=300

对于100%的数据，N <= 1500

Input

第一行两个整数N和M，分别表示行数和列数。接下来有N行，每行M个数，0或1分别表示这个格子是黑色或白色。

Output

仅有一行，表示满足条件最大正方形的 边长。

Sample Input

3 3

0 1 0

1 0 0

1 1 1

Sample Output

2

这是一道奇偶染色问题，需要考虑两种，把0变成1，剩下留给读者思考。

#include<stdio.h>
int a[1501][1501],f[1501][1501],f2[1501][1501];
int main()
{
    int n,m,i,j,min,ans=0,min2,ans2=0;
    scanf("%d%d",&n,&m);
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
        {
            scanf("%d",&a[i][j]);
            f[i][j]=a[i][j];
        }
    }
    for(i=1;i<=n;i++)
    {
    	for(j=1;j<=m;j++)
    	{
            if(a[i][j]==0)
            {
                if((i+j)%2==0)
                    f[i][j]=1;
                else
                    f[i][j]=0;
            }
            else
            {
                if((i+j)%2==0)
                    f[i][j]=0;
                else
                    f[i][j]=1;
            }  		
	    }
    }
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
        {
            if(f[i][j]==1)
            {
                if(f2[i-1][j-1]<=f2[i-1][j]&&f2[i-1][j-1]<=f2[i][j-1])
                {
                    f2[i][j]=f2[i-1][j-1]+1;
                }
                else if(f2[i-1][j]<=f2[i-1][j-1]&&f2[i-1][j]<=f2[i][j-1])
                {
                    f2[i][j]=f2[i-1][j]+1;
                }
                else if(f2[i][j-1]<=f2[i-1][j-1]&&f2[i][j-1]<=f2[i-1][j])
                {
                    f2[i][j]=f2[i][j-1]+1;
                }
            }
        }
    }
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
        {
            if(f2[i][j]>ans)
                ans=f2[i][j];
        }
    }    
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
        {
            if(f[i][j]==0)
            {
                if(f2[i-1][j-1]<=f2[i-1][j]&&f2[i-1][j-1]<=f2[i][j-1])
                {
                    f2[i][j]=f2[i-1][j-1]+1;
                }
                else if(f2[i-1][j]<=f2[i-1][j-1]&&f2[i-1][j]<=f2[i][j-1])
                {
                    f2[i][j]=f2[i-1][j]+1;
                }
                else if(f2[i][j-1]<=f2[i-1][j-1]&&f2[i][j-1]<=f2[i-1][j])
                {
                    f2[i][j]=f2[i][j-1]+1;
                }
            }
            else
            {
                f2[i][j]=0;
            }
        }
    }
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=m;j++)
        {
            if(f2[i][j]>ans2)
                ans2=f2[i][j];
        }
    }
    if(ans>ans2)
    	printf("%d",ans);
   	else
   		printf("%d",ans2);
}