九度OJ 1497：面积最大的全1子矩阵（DP）

#include<stdio.h>

int dp[1005][1005];//dp[i][j]代表位置(i,j)左上角矩阵之和
int a[1005][1005];

int main()
{
    int n,m;
    for(int i=0;i<=1000;i++)
        dp[0][i]=dp[i][0]=0;
    while(scanf("%d%d",&n,&m)>0)
    {
        for(int i=1;i<=n;i++)
            for(int j=1;j<=m;j++)
            scanf("%d",&a[i][j]);
        int sum=0;
        for(int i=1;i<=n;i++)//处理
        {
            sum=0;
            for(int j=1;j<=m;j++)
            {
               sum+=a[i][j]; dp[i][j]=dp[i-1][j]+sum;
            }
        }
        
        int maxk=0;
        for(int i=n;i>0;i--)
        for(int j=m;j>0&&maxk<i*j;j--)//枚举位置
        if(a[i][j])
        {
            int r=1,c=1;
            while(j-c>=0)//先只有一行时向左扩展
            {
                if(dp[i][j]-dp[i][j-c]-dp[i-r][j]+dp[i-r][j-c]==r*c)
                {
                     if(maxk<r*c) maxk=r*c; c++;
                }
                else
                    break;
            }
            while(i-r>=0&&c>0)//有多行时向现扩展，有可能列要往回缩
            {
                if(dp[i][j]-dp[i][j-c]-dp[i-r][j]+dp[i-r][j-c]==r*c)
                {
                    if(maxk<r*c) maxk=r*c; r++;
                }
                else
                    c--;
            }
        }
        printf("%d\n",maxk);
    }
}