最大子图形问题
CODEVS1159最大全0子矩阵
题目描述 Description
在一个0,1方阵中找出其中最大的全0子矩阵，所谓最大是指O的个数最多。
思路：这个题最朴素的n^6的算法，超时美美的。。。然后想优化，从一个点向上方、左方、右方扩展，首先更新这个点向上能有多少个0h0，然后找左右h比h0大的作为左右边界，然后计算这个矩形的面积，最后输出最大值。。。这种构造的美丽算法，真心。。。
比较： 最大全0子正方形：f[i][j](以i，j为右下角的最大正方形的边长)=min（f[i-1][j],f[i][j-1],f[i-1][j-1]）,这里用了正方形的特性，所以和矩形的求法不同。
这属于dp中的重要分支，最大子图形问题，详细的讲解可以参考下面的网址。。。真心丧病。。。
http://www.docin.com/p-46970779.html
code： #include<iostream> #include<cstdio> using namespace std; int li[2001]={0},ri[2001]={0},hi[2001]={0},a[2001][2001]={0}; int main() { int n,i,j,ans=0; scanf("%d",&n); for (i=1;i<=n;++i) for (j=1;j<=n;++j) scanf("%d",&a[i][j]); for (i=1;i<=n;++i) { for (j=1;j<=n;++j) { if (a[i][j]==0) ++hi[j]; else hi[j]=0; li[j]=ri[j]=j; } for (j=2;j<=n;++j) if (a[i][j]==0) while (hi[li[j]-1]>=hi[j]) li[j]=li[li[j]-1]; for (j=n-1;j>=1;--j) if (a[i][j]==0) while (hi[ri[j]+1]>=hi[j]) ri[j]=ri[ri[j]+1]; for (j=1;j<=n;++j) if (a[i][j]==0) ans=max(ans,(ri[j]-li[j]+1)*hi[j]); } cout<<ans<<endl; }
CODEVS1259最大正方形子矩阵
在一个01矩阵中，包含有很多的正方形子矩阵，现在要求出这个01矩阵中，最大的正方形子矩阵，使得这个正方形子矩阵中的某一条对角线上的值全是1，其余的全是0。
思路：做了好几个有关的最大子阵的问题，发现还是有些困难，做这个题想了好久，发现其实很简单，利用全0子矩阵的思路和正方形的思路就可以比较简单的写出dp方程。预处理一个点上方hi，左方li和右方ri0的个数（不包含这个点本身）。
f[i][j]=min(f[i-1][j-1]+1,min(li[i][j]+1,hi[i][j]+1))
第一个比较容易错的地方就来，每次对于能更新的f[i][j]的位置要求在map中为1，否则就不能构成要求的正方形；
其次就是题目中要求对角线为1，一个正方形有两条对角线，都应该考虑到，根据f数组的更新可以同理写出；
在预处理的时候细心，注意这个点和周围点的关系就可以了。
#include<iostream> #include<cstdio> using namespace std; int map[1010][1010]={0},li[1010][1010]={0},ri[1010][1010]={0},hi[1010][1010]={0},f[1010][1010]={0},g[1010][1010]={0}; int main() { int n,m,i,j,ans=0; cin>>n>>m; for (i=1;i<=n;++i) { for (j=1;j<=m;++j) { scanf("%d",&map[i][j]); hi[i][j]=hi[i-1][j]; li[i][j]=li[i][j-1]; if (map[i-1][j]==0) ++hi[i][j]; else hi[i][j]=0; if (map[i][j-1]==0&&j>1) ++li[i][j]; else li[i][j]=0; } for (j=m-1;j>=1;--j) { ri[i][j]=ri[i][j+1]; if (map[i][j+1]==0) ++ri[i][j]; else ri[i][j]=0; } } for (i=1;i<=n;++i) { for (j=1;j<=m;++j) if (map[i][j]==1) { f[i][j]=min(f[i-1][j-1]+1,min(li[i][j]+1,hi[i][j]+1)); if (f[i][j]>ans) ans=f[i][j]; } for (j=m;j>=1;--j) if (map[i][j]==1) { g[i][j]=min(g[i-1][j+1]+1,min(ri[i][j]+1,hi[i][j]+1)); if (g[i][j]>ans) ans=g[i][j]; } } cout<<ans<<endl; }
从网上看到了一个很全的子图形问题的Word，分享一下：http://www.docin.com/p-351795539.html
tyvj1563最大正方形
思路：这个题目很特殊，要求最大正方形中相邻两点的颜色不同，所以可以读图的时候进行处理，隔一个变一次（使这个位置上0、1互换），这样就转化成了矩阵中的最大全0或全1子正方形了，非常简单的dp解决
f[i][j]=min(f[i-1][j-1],min(f[i-1][j],f[i][j-1]))+1;(全0，先判断map[i][j]=0) 全1同理；
最大子图形问题的众多变式都应能转化为基本的求正方形、矩形等问题，得以比较简单的解决。看来noip后要找时间好好研究了。。。
bzoj1057 棋盘制作
题目大意：求一个01交错的最大正方形和矩形面积。
思路：把（i+j）%2==0的位置0/1互换一下，就变成了求最大0/1子矩阵问题，正方形可以根据长宽中取较小作为边长就可以了。
#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #define maxnode 2005 using namespace std; int map[maxnode][maxnode]={0},hi[maxnode]={0},li[maxnode]={0},ri[maxnode]={0}; int fang(int x){return x*x;} int main() { int i,j,ans1=0,ans2=0,k,n,m; scanf("%d%d",&n,&m); for (i=1;i<=n;++i) for (j=1;j<=m;++j) { scanf("%d",&map[i][j]); if ((i+j)%2) map[i][j]^=1; } for (k=0;k<=1;++k) { memset(hi,0,sizeof(hi)); for (i=1;i<=n;++i) { for (j=1;j<=m;++j) { if (map[i][j]==k) ++hi[j]; else hi[j]=0; li[j]=ri[j]=j; if (j>1&&map[i][j]==k) while(hi[li[j]-1]>=hi[j]) li[j]=li[li[j]-1]; } for (j=m-1;j>=1;--j) if (map[i][j]==k) while(hi[ri[j]+1]>=hi[j]) ri[j]=ri[ri[j]+1]; for (j=1;j<=m;++j) if (map[i][j]==k) { ans1=max(ans1,fang(min(hi[j],ri[j]-li[j]+1))); ans2=max(ans2,hi[j]*(ri[j]-li[j]+1)); } } } printf("%d\n%d\n",ans1,ans2); }
