agc 033 D-Complexity

https://atcoder.jp/contests/agc033/tasks/agc033_d

場上沒想出來...做法：可以發現答案上界是log(H*W)的，所以考慮暴力枚答案，那麼的當前狀態可以表示爲選取(r,c)和(r',c)兩個點能擴展到最右的c'在哪裏以及選取(r,c)和(r,c')兩個點能擴展到最下面的r'在哪裏，所以一層裏狀態數是n^3的。對於轉移考慮由於它的複雜度計算是取max，所以每次更新可以O(1)的類似倍增的更新，不過橫着切和豎着切還要相互之間更新最遠位置，這一步複雜度貌似我不太會證明，不過也有n^4/4的上界5s不虛，而且應該是不可能卡滿的......

代碼：

#include<bits/stdc++.h>
using namespace std;
const int N=190;

int n,m;
char mp[N][N];
int f[N][N][N],g[N][N][N],sum[N][N];
int nf[N][N][N],ng[N][N][N];

bool all_same(int h1,int h2,int l1,int l2)
{
	int x=(h2-h1+1)*(l2-l1+1),y=sum[h2][l2]-sum[h2][l1-1]-sum[h1-1][l2]+sum[h1-1][l1-1];
	return (y==0)||(x==y);
}

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++)
		scanf("%s",mp[i]+1);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			sum[i][j]=sum[i][j-1]+sum[i-1][j]-sum[i-1][j-1]+(mp[i][j]=='.');
	for(int h1=1;h1<=n;h1++)
	{
		for(int h2=h1;h2<=n;h2++)
		{
			for(int l1=1;l1<=m;l1++)
			{
				int l=l1-1,r=m,mid;
				while(l<r)
				{
					mid=(l+r+1)/2;
					if(all_same(h1,h2,l1,mid))l=mid;
					else r=mid-1;
				}
				f[h1][h2][l1]=l;
			}
		}
		for(int l1=1;l1<=m;l1++)
		{
			for(int l2=l1;l2<=m;l2++)
			{
				int l=h1-1,r=n,mid;
				while(l<r)
				{
					mid=(l+r+1)/2;
					if(all_same(h1,mid,l1,l2))l=mid;
					else r=mid-1;
				}
				g[l1][l2][h1]=l;
			}
		}
	}
	for(int ans=0;;++ans)
	{
		if(f[1][n][1]==m)return printf("%d\n",ans),0;
		for(int h1=1;h1<=n;h1++)
		{
			for(int h2=h1;h2<=n;h2++)
			{
				for(int l1=1,nx;l1<=m;l1++)
				{
					nx=f[h1][h2][l1];
					if(nx<m)nx=f[h1][h2][nx+1];
					nf[h1][h2][l1]=nx;
				}
			}
			for(int l1=1;l1<=m;l1++)
			{
				for(int l2=l1,nx;l2<=m;l2++)
				{
					nx=g[l1][l2][h1];
					if(nx<n)nx=g[l1][l2][nx+1];
					ng[l1][l2][h1]=nx;
				}
			}
		}
		for(int h1=1;h1<=n;h1++)
		{
			for(int h2=h1;h2<=n;h2++)
			{
				for(int l1=1,r;l1<=m;l1++)
				{
					r=nf[h1][h2][l1];
					while(r<m)
					{
						if(ng[l1][r+1][h1]>=h2)++r;
						else break;
					}
					f[h1][h2][l1]=r;
				}
			}
			for(int l1=1;l1<=m;l1++)
			{
				for(int l2=l1,r;l2<=m;l2++)
				{
					r=ng[l1][l2][h1];
					while(r<n)
					{
						if(nf[h1][r+1][l1]>=l2)++r;
						else break;
					}
					g[l1][l2][h1]=r;
				}
			}
		}
	}
}