创意吃鱼法

题目描述

回到家中的猫猫把三桶鱼全部转移到了她那长方形大池子中，然后开始思考：到底要以何种方法吃鱼呢（猫猫就是这么可爱，吃鱼也要想好吃法 ^_*）。她发现，把大池子视为01矩阵（0表示对应位置无鱼，1表示对应位置有鱼）有助于决定吃鱼策略。

在代表池子的01矩阵中，有很多的正方形子矩阵，如果某个正方形子矩阵的某条对角线上都有鱼，且此正方形子矩阵的其他地方无鱼，猫猫就可以从这个正方形子矩阵“对角线的一端”下口，只一吸，就能把对角线上的那一队鲜鱼吸入口中。

猫猫是个贪婪的家伙，所以她想一口吃掉尽量多的鱼。请你帮猫猫计算一下，她一口下去，最多可以吃掉多少条鱼？

输入输出格式

输入格式：
有多组输入数据，每组数据：

第一行有两个整数n和m（n,m≥1），描述池塘规模。接下来的n行，每行有m个数字（非“0”即“1”）。每两个数字之间用空格隔开。

对于30%的数据，有n,m≤100

对于60%的数据，有n,m≤1000

对于100%的数据，有n,m≤2500

输出格式：
只有一个整数——猫猫一口下去可以吃掉的鱼的数量，占一行，行末有回车。

输入输出样例

输入样例#1：
4 6
0 1 0 1 0 0
0 0 1 0 1 0
1 1 0 0 0 1
0 1 1 0 1 0
输出样例#1：
3
说明

右上角的

1 0 0 0 1 0 0 0 1
动态规划，f[i][j]表示以(i,j)为顶点的最大子矩阵中能吃到的鱼的条数。

f[i][j]可以从f[i-1][j-1]或f[i-1][j+1]转移而来，取其中最大值。

每次转移需要横向和纵向检查0的个数，取两者最小值k即得到了最大子矩阵中鱼的个数。

C代码：

动态规划，f[i][j]表示以(i,j)为顶点的最大子矩阵中能吃到的鱼的条数。

f[i][j]可以从f[i-1][j-1]或f[i-1][j+1]转移而来，取其中最大值。

每次转移需要横向和纵向检查0的个数，取两者最小值k即得到了最大子矩阵中鱼的个数。

C代码：

#include<stdio.h>
int n,m,i,j,k,k1,k2,max,temp,f[2501][2501];
int main(void)
{
    while(scanf("%d%d",&n,&m)==2)
    {
        max=0;//答案 
        for(i=1;i<=n;i++)
        {
            for(j=1;j<=m;j++)scanf("%d",&f[i][j]);//先输入一行，因为下面要横向检查0的个数。 
            for(j=1;j<=m;j++)
            {
                if(f[i][j]==1) 
                {
                    temp=f[i-1][j-1];//检查左上方 
                    if(temp>=1 && temp+1>f[i][j])
                    {
                        for(k1=1;k1<=temp;k1++)//横向往左检查0的个数k1，最多为temp个 
                            if(f[i][j-k1]>=1) break;
                        for(k2=1;k2<=temp;k2++)//纵向往上检查k2…… 
                            if(f[i-k2][j]>=1) break; 
                        f[i][j]=k1<k2?k1:k2; //取k1、k2较小的那个，就是满足条件的最大子矩阵。 
                    }
                    temp=f[i-1][j+1];//检查右上方 
                    if(temp>=1 && temp+1>f[i][j])
                    {
                        for(k1=1;k1<=temp;k1++)
                            if(f[i][j+k1]>=1) break;//唯一不同是这里的j+k1而不是j-k1,因为是往右检查 
                        for(k2=1;k2<=temp;k2++)
                            if(f[i-k2][j]>=1) break;
                        f[i][j]=k1<k2?k1:k2;
                    }
                }
                if(f[i][j]>max) max=f[i][j];
            }
        }
        printf("%d\n",max);        
    }
    return 0;
}


#include<iostream>
#include<cstring>
using namespace std;
int f[3000][3000];
int main()
{  int m,n,maxx=0,tmp;
   cin>>n>>m;
   memset(f,127,sizeof(f));
   for (int i=1;i<=n;++i)
   for (int j=1;j<=m;++j)
       cin>>f[i][j]; 
   for (int i=1;i<=n;++i)
   {  
       for (int j=1;j<=m;++j)
      {  if (f[i][j])
        {  tmp=f[i-1][j-1];
           if (tmp)
           {  int k,c=0;
              for (k=1;k<=tmp;++k)
                if (f[i-k][j]||f[i][j-k])
                  {   c=1;
                   break;
                  }
             if (c)
             f[i][j]=k;
             else
            f[i][j]=tmp+1;

           }
           tmp=f[i-1][j+1];
            if (tmp)
             {  int k,c=0;
              for (k=1;k<=tmp;++k)
                 if (f[i-k][j]||f[i][j+k])
                 {  c=1;
                   break;
                 }
                if (c)
                 f[i][j]=max(f[i][j],k);
                else f[i][j]=max(f[i][j],tmp+1); 
             }
            if (f[i][j]>maxx)
            maxx=f[i][j];
         }
        }
   }
  cout<<maxx;
 }