Matrix Swapping II（ＤＰ，最大矩形问题）

http://acm.hdu.edu.cn/diy/diy_previewproblem.php?cid=19572&pid=1020
Matrix Swapping II

Time Limit : 9000/3000ms (Java/Other)   Memory Limit : 32768/32768K (Java/Other)
Total Submission(s) : 1   Accepted Submission(s) : 1
Font: Times New Roman | Verdana | Georgia
Font Size: ← →
Problem Description
Given an N * M matrix with each entry equal to 0 or 1. We can find some rectangles in the matrix whose entries are all 1, and we define the maximum area of such rectangle as this matrix’s goodness.

We can swap any two columns any times, and we are to make the goodness of the matrix as large as possible.
Input
There are several test cases in the input. The first line of each test case contains two integers N and M (1 ≤ N,M ≤ 1000). Then N lines follow, each contains M numbers (0 or 1), indicating the N * M matrix
Output
Output one line for each test case, indicating the maximum possible goodness.
Sample Input
3 4
1011
1001
0001
3 4
1010
1001
0001
Sample Output
4
2

Note: Huge Input, scanf() is recommended.
Source
2009 Multi-University Training Contest 2 - Host by TJU
解析：求最大矩形，跟hdu1505，1506类似，
只不过这里的整列可以进行交换。
故可以借助一个数组保存当前情况，并且对其进行排序，然后按照hdu1506的方法进行计算
*/

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<iostream>
using namespace std;
const int maxn=1000+5;
double  h[maxn],temp[maxn];
int r[maxn],l[maxn];
int n,m;
char ch[maxn];
void init()//寻找大于等于h[i]的边界点
{
    int i;
    for(i=1;i<=m;i++)
     {r[i]=l[i]=i;
      h[i]=temp[i];
     }
     sort(h+1,h+m+1);
    h[0]=h[m+1]=-1;
    for(i=1;i<=m;i++)
     {
         if(h[l[i]-1]>=h[i])//寻找左边界点
          {
              while(h[l[i]-1]>=h[i])
                {l[i]=l[l[i]-1];

                }
            }
    }
 for(i=m;i>=1;i--)
    {
        if(h[r[i]+1]>=h[i])
         {
             while(h[r[i]+1]>=h[i])//寻找右边界点
              {r[i]=r[r[i]+1];
              }
         }
    }
}
int main()
{
    int i,j,k;
     int t;
     int ans;
     while(scanf("%d%d",&n,&m)!=EOF)
       {
        memset(h,0,sizeof(h));
        memset(temp,0,sizeof(temp));
        ans=0;
        for(i=1;i<=n;i++)
        {
             scanf("%s",ch);
            for(j=0;j<m;j++)//以行为底进行其高度
             {
                 if(ch[j]=='1')
                 temp[j+1]++;
                 else//如果不连续则将本行j+1列置为0
                 temp[j+1]=0;
             }
             init();
             for(k=1;k<=m;k++)
              {
                  t=h[k]*(r[k]-l[k]+1);
                  if(ans<t)
                  ans=t;
              }

        }
        printf("%d\n",ans);
    }
    return 0;
}