【DP+记忆化搜索】PKU-1191-棋盘分割

题目

sum[x0][y0][x1][y1]表示的是左上角为（x0，y0），右上角为（x1，y1）的矩形的每个点分数总和的平方；dp[n][x0][y0][x1][y1]表示的是左上角为（x0，y0），右上角为（x1，y1）的矩形再经n刀切到的方差的平方加上平均值的平方的和再乘以n……

状态转移方程为dp(n, x0, y0, x1, y1)=Min( Min( dp[n-1, x0, y0, a ,y1]+sum[a+1, y0, x1, y1], dp[n-1, a+1, y0, x1, y1]+sum[x0, y0, a, y1]), Min(dp[n-1, x0, y0, x1, b]+sum[x0, b+1, x1, y1], dp[n-1, x0, b+1, x1, y1]+sum[x0, y0, x1, b]))。其中x0<=a<x1, y0<=b<y1。

#include<iostream>
#include<cmath>
using namespace std;
template<class T> T Min(T x,T y){return x<y?x:y;}
int sum[10][10][10][10],dp[15][10][10][10][10];
int dfs(int n,int x0,int y0,int x1,int y1)
{
	int i,d,&ret=dp[n][x0][y0][x1][y1];
	if(ret!=-1)return ret;
	if(n==0)return ret=sum[x0][y0][x1][y1];
	d=999999999;
	for(i=x0;i<x1;i++)
	{
		d=Min(d,dfs(n-1,x0,y0,i,y1)+sum[i+1][y0][x1][y1]);                           
		d=Min(d,dfs(n-1,i+1,y0,x1,y1)+sum[x0][y0][i][y1]);
	}
	for(i=y0;i<y1;i++)
	{
		d=Min(d,dfs(n-1,x0,y0,x1,i)+sum[x0][i+1][x1][y1]);
		d=Min(d,dfs(n-1,x0,i+1,x1,y1)+sum[x0][y0][x1][i]);
	}
	return ret=d;
}
int main()
{
	//freopen("a.txt","r",stdin);
	int i,j,k,l,m,n,sum0,map[10][10];
	double d;
	while(scanf("%d",&n)!=EOF)
	{
		for(i=1;i<=8;i++)
			for(j=1;j<=8;j++)
				scanf("%d",&map[i][j]);
		memset(dp,-1,sizeof(dp));
		memset(sum,0,sizeof(sum));
		sum0=0;
		for(i=1;i<=8;i++)
			for(j=1;j<=8;j++)
			{
				sum0+=map[i][j];
				for(k=i;k<=8;k++)
					for(l=j;l<=8;l++)
					{
						for(m=i;m<k;m++)sum[i][j][k][l]+=map[m][l];
						for(m=j;m<l;m++)sum[i][j][k][l]+=map[k][m];
						sum[i][j][k][l]+=sum[i][j][k-1][l-1]+map[k][l];
					}
			}
		for(i=1;i<=8;i++)
			for(j=1;j<=8;j++)
				for(k=i;k<=8;k++)
					for(l=j;l<=8;l++)
						sum[i][j][k][l]*=sum[i][j][k][l];
		d=sum0*1.0/n;
		d*=d;
		printf("%.3f\n",sqrt(dfs(n-1,1,1,8,8)*1.0/n-d));
	}
	return 0;
}