九度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++;}elsebreak;}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++;}elsec--;}}printf("%d\n",maxk);}
}