NYOJ 87-棋盘分割（记忆化搜索）

棋盘分割

时间限制： 1000 ms  |  内存限制： 65535 KB

难度： 6

描述

将一个８*８的棋盘进行如下分割：将原棋盘割下一块矩形棋盘并使剩下部分也是矩形，再将剩下的部分继续如此分割，这样割了(n-1)次后，连同最后剩下的矩形棋盘共有n块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行) 

原棋盘上每一格有一个分值，一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘，并使各矩形棋盘总分的均方差最小。 
均方差 ，其中平均值 ，x _i为第i块矩形棋盘的总分。 
请编程对给出的棋盘及n，求出O'的最小值。 

输入

第1行为一个整数n(1 < n < 15)。 
第2行至第9行每行为8个小于100的非负整数，表示棋盘上相应格子的分值。每行相邻两数之间用一个空格分隔。

输出

仅一个数，为O'（四舍五入精确到小数点后三位）。

样例输入

3
1 1 1 1 1 1 1 3
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 0
1 1 1 1 1 1 0 3

样例输出

1.633

来源

IOI 2000

上传者

hzyqazasdf

一开始理解错题意了，以为切完后所有块都要操作，后来才发现一分为二时，其中一块便不能再动了，呢就是由大化小的问题了，动态规划随便搞搞，感觉记忆化搜索会好写一点。。。。。

每次四种选择，分别搜索：

（1）横着切，选上边一块

（2）横着切，选下边一块

（3）竖着切，选左边一块

（4）竖着切，选右边一块

剩下的就是dp了，定义dp[kk][i][j][k][h]:左上角（i，j）到右下角（k,h）组成的矩形，切到第kk次的最优解。

#include<map>          
#include<stack>          
#include<queue>          
#include<vector>          
#include<math.h>          
#include<stdio.h>          
#include<iostream>          
#include<string.h>          
#include<stdlib.h>  
#include<algorithm> 
#include<functional>  
using namespace std;          
typedef long long  ll;          
#define inf  1000000000         
#define MOD 1000000007           
#define  maxn  5100005
#define  lowbit(x) (x&-x)          
#define  eps 1e-9
int dp[20][10][10][10][10],a[10][10][10][10],b[10][10];
int sm(int x,int y,int xx,int yy)
{
	int ans=0,i,j;
	for(i=x;i<=xx;i++)
		for(j=y;j<=yy;j++)
			ans+=b[i][j];
	return ans;
}
int dfs(int k,int x,int y,int xx,int yy)
{
	if(dp[k][x][y][xx][yy]!=-1)
		return dp[k][x][y][xx][yy];
	int i,j,t1,t2,t;
	dp[k][x][y][xx][yy]=inf;
	if(xx>x)
	{
		for(i=x;i<xx;i++)
		{
			t1=dfs(k-1,i+1,y,xx,yy);
			t2=dfs(k-1,x,y,i,yy);
			t=min(t1+a[x][y][i][yy],t2+a[i+1][y][xx][yy]);
			dp[k][x][y][xx][yy]=min(dp[k][x][y][xx][yy],t);
		}
	}
	if(yy>y)
	{
		for(i=y;i<yy;i++)
		{
			t1=dfs(k-1,x,i+1,xx,yy);
			t2=dfs(k-1,x,y,xx,i);
			t=min(t1+a[x][y][xx][i],t2+a[x][i+1][xx][yy]);
			dp[k][x][y][xx][yy]=min(dp[k][x][y][xx][yy],t);
		}
	}
	return dp[k][x][y][xx][yy];
}
int main(void)
{
	int x,y,xx,yy,n,i,j;
	scanf("%d",&n);
	for(i=1;i<=8;i++)
		for(j=1;j<=8;j++)
			scanf("%d",&b[i][j]);
	memset(dp,-1,sizeof(dp));
	for(x=1;x<=8;x++)
		for(xx=x;xx<=8;xx++)
			for(y=1;y<=8;y++)
				for(yy=y;yy<=8;yy++)
				{
					int tmp=sm(x,y,xx,yy);
					dp[1][x][y][xx][yy]=a[x][y][xx][yy]=tmp*tmp;
				}
	dfs(n,1,1,8,8);
	double ans,sum;
	sum=1.0*sm(1,1,8,8);
	sum=(sum/n)*(sum/n);
	ans=sqrt(1.0*dp[n][1][1][8][8]/n-sum);
	printf("%.3f\n",ans);
	return 0;
}