洛谷 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<iostream>
#include<stdio.h>
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;
}