洛谷 P1169 [ZJOI2007]棋盘制作(bzoj P1057 [ZJOI2007]棋盘制作)

传送门

学习一发悬线法，但是好像出了点意外的事，在洛谷上A了却在BZOJ上炸了然后和标解拍了没问题。。有毒！！(好吧发现了一个弱智错误有一个m打成了n

悬线法：

设line[i][j]表示[i][j]位置向上的最大悬线(就是到障碍点的最大长度)，fx[i][j]表示[i][j]位置悬线向左最大能移动距离，fy[i][j]表示[i][j]位置悬线向右最大能移动距离。怎么求fx和fy呢？？如果悬线长度为1那么就是这个位置向左或向右达到障碍的距离，若悬线不为1那么就是fx[i][j]=min(fx[i-1][j],fx1[i][j])，fy[i][j]=min(fy[i-1][j],fy1[i][j])，fx1，fx2表示[i][j]位置向左向右遇到第一个障碍的距离。然后显然答案就是max(line[i][j]*(fx[i][j]+fy[i][j]-1))。

UPDATE（7.19）：我们来理解一下fx1,fy1,fx,fy其实都可以理解为单调栈

代码：

#include
#include
using namespace std;
const int oo=2147483647;
int g[2005][2005],f[2005][2005],line[2005][2005],fx1[2005][2005],fy1[2005][2005],fx[2005][2005],fy[2005][2005];
int main()
{
	int n,m,ans1=-oo,ans2=-oo;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&g[i][j]);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			line[i][j]=1;
			fx1[i][j]=1;
			fy1[i][j]=1;
			fx[i][j]=1;
			fy[i][j]=1;
			f[i][j]=1;
		}
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			if(g[i][j-1]!=g[i][j])
				fx1[i][j]=fx1[i][j-1]+1;
			if(g[i-1][j]!=g[i][j])
				line[i][j]=line[i-1][j]+1;
		}
	for(int i=1;i<=n;i++)
		for(int j=m;j>=1;j--)
			if(g[i][j+1]!=g[i][j])
				fy1[i][j]=fy1[i][j+1]+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if(line[i][j]!=1)
			{
				fx[i][j]=min(fx[i-1][j],fx1[i][j]);
				fy[i][j]=min(fy[i-1][j],fy1[i][j]);
			}
			else
			{
				fx[i][j]=fx1[i][j];
				fy[i][j]=fy1[i][j];
			}
	for(int i=2;i<=n;i++)
		for(int j=2;j<=m;j++)
			if(g[i][j]==g[i-1][j-1]&&g[i][j]!=g[i-1][j]&&g[i][j]!=g[i][j-1])
				f[i][j]=min(f[i-1][j-1],min(f[i-1][j],f[i][j-1]))+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
		{
			ans1=max(ans1,f[i][j]);
			ans2=max(ans2,line[i][j]*(fx[i][j]+fy[i][j]-1));
		}
	printf("%d\n%d\n",ans1*ans1,ans2);
	return 0;
}