bzoj 1768: [Ceoi2009]logs

Description

有一个N*M的01矩阵，现在你可以的任意交换其中的列，要求找一个最大的仅由1组成的矩阵。 N<=15000,M<=1500 1 ≤ N ≤ 15000 1 ≤ M ≤ 1500 30% of the test cases will have N,M ≤ 1024

Input

N,M 以下N行每行M个字符。

Output

Sample Input

10 6
001010
111110
011110
111110
011110
111111
110111
110111
000101
010101

Sample Output

21

HINT

By permuting the columns such that columns 2, 
4 and 5 are adjacent you have a rectangle of area 

21 (rows 2-8 and columns 2, 4, 5). 

ceoi是中欧么...连题解都只有一篇还看不太懂

http://blog.sina.com.cn/s/blog_7cd3ac450100skzu.html

我们用f[i][j]表示该列最多向下延伸的长度。

若s[i][j]=0 则f[i][j]=0，否则f[i][j]=f[i-1][j]+1 【因为统一列的顺序不变】

然后用t[i]表示横向最多延伸的长度,q数组记录1的位置

然后因为内存只有64M所以要滚动数组

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
using namespace std;
int i,j,n,ans,l,r,m;
string s;
int f[2][1501],q[2][1501];
bool d[2][1501];
int t[2];
int main()
{
     scanf("%d%d",&n,&m);
     for(i=1;i<=n;i++)
     {
          l=i&1;
		  r=(i+1)&1;
		  t[l]=0;
		  cin>>s;
          for(j=1;j<=m;j++)
          {
               if(s[j-1]=='1')
               {
                    d[l][j]=true;
                    if(!d[r][j])
                    {
                         t[r]++;
                         q[r][t[r]]=j;
                         f[l][j]=1;
                    }
					else
					     f[l][j]=f[r][j]+1;
			   }
               else
               {
                    d[l][j]=false;
					f[l][j]=0;
			   }
          }
          for(j=1;j<=t[r];j++)
          {
               if(d[l][q[r][j]])
               {
                    t[l]++;
                    q[l][t[l]]=q[r][j];
                    if(f[l][q[r][j]]*t[l]>ans)
					     ans=f[l][q[r][j]]*t[l];
               }
          }
     }
     printf("%d\n",ans);
     return 0;
}