[bzoj1057][ZJOI2007]棋盘制作

传送门

图上的棋子有两种类型：
1.黑格行列奇偶性相同，白格不同
2.白格行列奇偶性相同，黑格不同

判断同不同的话，直接异或就可以了。
然后问题就转化为求最大01子矩阵了
感觉这个真的挺巧妙的，要是不转化的话，大概会写吐？

以下是一些废话。。（没错我越来越话痨了
统计矩形本来写了悬线法，然后有个地方忘记了怎么写，就写错了，结果竟然只有vijos上TLE（没错我就是有每个oj交一遍的癖好）。洛谷的数据实在太水，bzoj大概是由于算总时间的吧。然后重新学了悬线法觉得好烦啊就写了单调队列。然而我以前也没写过单调队列的方法，就听别人讲过一次，但是忘得差不多了，而且当时好像也一知半解，只知其然而不知其所以然。于是我又成功写错了。。但是后来自己仔细推了一下，在晚自习最后十分钟切了。果然一个算法还是要自己想一遍写一遍才能真正理解啊！（然而我就是太懒）

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=2001;
int n,m,a[N][N],f[N][N],up[N][N],jx[N],h[N],L,ans1,ans2;

void read(int &x){
    char ch=getchar();x=0;
    for(;ch<'0'||ch>'9';ch=getchar());
    for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0';
}

int main(){
    read(n);read(m);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
     	read(a[i][j]);
     	a[i][j]=(i%2)^(j%2)^a[i][j];
     }
//	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) cout<<a[i][j]<<' ';
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
     	if (a[i][j]) continue;
     	int mn=min(f[i-1][j],f[i][j-1]);
     	if (!a[i-mn][j-mn]) f[i][j]=mn+1;
     	 else f[i][j]=mn;
     	ans1=max(ans1,f[i][j]);
     }
    memset(f,0,sizeof f);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
     	if (!a[i][j]) continue;
     	int mn=min(f[i-1][j],f[i][j-1]);
     	if (a[i-mn][j-mn]) f[i][j]=mn+1;
     	 else f[i][j]=mn;
     	ans1=max(ans1,f[i][j]);
     }
     
    for(int o=0;o<=1;o++){
        memset(up,0,sizeof up);
        for(int i=1;i<=n;i++){
            L=0;
            for(int j=1;j<=m;j++){
                if (a[i][j]==o) up[i][j]=up[i-1][j]+1;
                jx[++L]=j;h[L]=up[i][j];
                for(;L>1&&h[L-1]>=h[L];L--){
                    ans2=max(ans2,h[L-1]*(j-jx[L-1]));
                    h[L-1]=h[L];
                }
            }
            ans2=max(ans2,h[L]*(m-jx[L]+1));
            for(int j=L-1;j>0;j--) ans2=max(ans2,h[j]*(m-jx[j]+1));
        }
    }
    cout<<ans1*ans1<<endl<<ans2;
    return 0;
}

还有一个辣鸡并且是假的悬线法，但是在洛谷和bzoj可以过

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int N=2001;
int n,m,a[N][N],f[N][N],up[N][N],L[N][N],r[N][N],ans1,ans2;

void read(int &x){
    char ch=getchar();x=0;
    for(;ch<'0'||ch>'9';ch=getchar());
    for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<3)+(x<<1)+ch-'0';
}

int main(){
    read(n);read(m);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
     	read(a[i][j]);
     	a[i][j]=(i%2)^(j%2)^a[i][j];
     }
//	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) cout<<a[i][j]<<' ';
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
     	if (a[i][j]) continue;
     	int mn=min(f[i-1][j],f[i][j-1]);
     	if (!a[i-mn][j-mn]) f[i][j]=mn+1;
     	 else f[i][j]=mn;
     	ans1=max(ans1,f[i][j]);
     }
    memset(f,0,sizeof f);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++){
     	if (!a[i][j]) continue;
     	int mn=min(f[i-1][j],f[i][j-1]);
     	if (a[i-mn][j-mn]) f[i][j]=mn+1;
     	 else f[i][j]=mn;
     	ans1=max(ans1,f[i][j]);
     }
     
    for(int o=0;o<=1;o++){
        memset(up,0,sizeof up);
        for(int i=1;i<=n;i++)
         for(int j=1;j<=m;j++)
          if (a[i][j]==o) up[i][j]=up[i-1][j]+1;
        memset(L,0,sizeof L);
        for(int i=1;i<=n;i++)
         for(int j=1;j<=m;j++)
          if (up[i][j-1]>=up[i][j]&&a[i][j]==o){
          	int k;
          	for(k=j-1;up[i][k-1]>=up[i][j];k-=L[i][k-1]+1);
          	L[i][j]=j-k;
          }
        memset(r,0,sizeof r);
        for(int i=1;i<=n;i++)
         for(int j=1;j<=m;j++)
          if (up[i][j+1]>=up[i][j]&&a[i][j]==o){
          	int k;
          	for(k=j+1;up[i][k+1]>=up[i][j];k+=r[i][k+1]+1);
          	r[i][j]=k-j;
          }
        for(int i=1;i<=n;i++)
         for(int j=1;j<=m;j++)
          ans2=max(ans2,(L[i][j]+1+r[i][j])*up[i][j]);
    }
//	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) cout<<r[i][j]<<' ';
/*	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) cout<<L[i][j]<<' ';
    cout<<endl;
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) cout<<up[i][j]<<' ';*/
    cout<<ans1*ans1<<endl<<ans2;
    return 0;
}