CD16 求最大子矩阵的大小

题目描述
给定一个整型矩阵 map，其中的值只有 0 和 1 两种，求其中全是 1 的所有矩形区域中，最大的矩形区域里 1 的数量。
输入描述:
第一行输入两个整数 n 和 m，代表 nm 的矩阵
接下来输入一个 nm 的矩阵
输出描述:
输出其中全是 1 的所有矩形区域中，最大的矩形区域里 1 的数量。
示例1
输入
复制
1 4
1 1 1 0
输出
复制
3
说明
最大的矩形区域有3个1，所以返回3

 #include<stdio.h>
 #include<iostream>
 #include<algorithm>
 #include<stack>
 using namespace std;
 
//思路：stack_1为一个单调栈，记 arr[i] 为当前需要放入栈中的元素,当前栈顶元素为 height[j]
// (1) 如height[j]>arr[i],直接放入。
// (2) height[j]<=arr[i],弹出当前元素，k为新的栈顶元素。思考位置j为柱子能扩展到那里
//       右边i-1,左边k+1  所以，j位置的柱子能扩出来的最大矩阵为[ (i-1)-(k+1)+1 ]*height[j]
//(注：当height[j]==arr[i],时右边不一定时i-1，但是此时height[j]为柱子的最大矩形和 height[i]为柱子的最大矩形完全一样)        
 
 

 int n,m;
 int arr[2007][2007]; 
 int height[2007];
 
 
 int maxfun( int height[] )
 {
 	stack<int> stack_1;//存height[]数组的下标 
 	int maxArea=-1,curArea=0;
 	for(int i=0;i<m;i++)
	 {
	 	while( !stack_1.empty() && height[ stack_1.top() ]>= height[i]  )
		 {
		 	int j=stack_1.top();
		 	stack_1.pop();
		 	int k= stack_1.empty()? -1 : stack_1.top();
		 	curArea=(i-1-k)*height[j];//以heirht[j]为柱子的最大面积 
		 	maxArea=max(maxArea , curArea );
		 }
		 stack_1.push( i );
	 } 
	 
	 while( !stack_1.empty() )
	 {
	 	int j=stack_1.top();
	 	stack_1.pop();
	 	int k= stack_1.empty()? -1: stack_1.top();
	 	curArea=(m-k-1)*height[j];//这里m是height[]数组的长度。
	    maxArea=max(maxArea , curArea );	
	 }
	 
	 return maxArea;
 }
 
 
 int main()
 {
 	int res=-1;
 	scanf("%d%d",&n,&m);
 	for(int i=0;i<n;i++)
	 {
	 	for(int j=0;j<m;j++)
	 	   scanf("%d",&arr[i][j]);
	 }
	
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<m;j++)
		{
			height[j]= arr[i][j]==0 ? 0 : height[j]+1; 
		}
//		cout<<"height[]数组为：";
//		for(int z=0;z<m;z++)
//		   cout<<height[z]<<"  "; 
		res=max(res,maxfun( height ) );
	}
	
	cout<<res<<endl;
	
	 
	
 	
 	
 	
 	return 0;
 }