单位矩阵

Description
X是一个灰常讨厌数学的学生，所以想请你帮他解决一个简单的关于矩阵的问题。给定一个只有0和1构成的大小为M*N的矩阵，在其中找到最大的子单位矩阵。单位矩阵是指主对角线上元素全部为1，其余元素全部为0的方阵。主对角线是指左上角到右下角方向的对角线。数据规模(1<=M，N<=1000)。

Input
多组数据输入，每组第一行输入M和N，空格分开，接下来M行输入矩阵数据，每行N个元素，输入以文件尾(EOF)结束。

Output
对于每组数据，输出该矩阵中最大子单位矩阵的行数，每次输出独占一行。

Sample Input

5 6
010000
010001
010000
001000
100111
2 5
01000
00110

Sample Output

3
2

这道题目真是有点尴尬呀，直接暴力（n^3 m^3）或者( n^2 m^2),尝试了一下ETl了，真是太难受了，感觉没有什么好的方法，于是参考了一个大犇的写法，真是太巧妙了，dp，这个东西本蒻蒻还没有能力能够很好的运用，这里就流下的么技术的眼泪v_v
这里的话我们可以用dp一层一层往下找，当然就得先初始化0，找的过程就是，求左面的长度和高度的最小值与左斜上方的最小值，这个用到的就是木桶效应，有了这个状态方程问题不就迎刃而解了嘛.

#include<bits/stdc++.h>
using namespace std;
int dp[1005][1005];
int up[1005];
int main(){
    int n,m,left,ma;
    char ch;
    while(~scanf("%d%d",&n,&m)){
        ma=0;
        memset(dp,0,sizeof(dp));
        memset(up,0,sizeof(up));
        for(int i=1;i<=n;i++){
            left=0;
            for(int j=1;j<=m;j++){
                cin>>ch;
                if(ch=='1'){
                    dp[i][j]=min(dp[i-1][j-1],min(left,up[j]))+1;
                    left=0;
                    up[j]=0;
                    ma=max(ma,dp[i][j]);
                }
                else{
                    left++;
                    up[j]++;
                }
            }
        }
        printf("%d\n",ma);
    }
    return 0;
}