棋盘分割

题目：

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

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

Input

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

Output

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

Sample Input

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

Sample Output

1.633

 

思路：

8*8的棋盘，数据量较小，可以直接暴力。

要使各矩形棋盘总分的均方差最小，只需求分割后各块平方和的最小值。

先计算出每个区域只分成一块时的平方和，然后暴力枚举把每一块分成1，2..c份储存在dp[i][j][h][r][c]数组中。

最后dp[1][1][8][8][n]即为分成n份后各块平方和的最小值。

最后算出各矩形棋盘总分的均方差最小的值。

 

代码：

#include<math.h>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#define N 9
using namespace std;
int n,dp[N][N][N][N][2*N],k[N][N],s[N][N];
int main()
{
    while(~scanf("%d",&n))
    {
        memset(dp,0x3f3f3f3f,sizeof(dp));
        for(int i=1; i<=8; i++)
            for(int j=1; j<=8; j++)
                scanf("%d",&k[i][j]);
        for(int i=1; i<=8; i++)
            for(int j=1; j<=8; j++)
                s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+k[i][j];
        for(int i=1; i<=8; i++)
            for(int j=1; j<=8; j++)
                for(int h=1; h<=8; h++)
                    for(int r=1; r<=8; r++)
                    {
                        int t=s[h][r]-s[h][j-1]-s[i-1][r]+s[i-1][j-1];
                        dp[i][j][h][r][1]=t*t;
                    }
        for(int c=2; c<=n; c++)
            for(int i=1; i<=8; i++)
                for(int j=1; j<=8; j++)
                    for(int h=i; h<=8; h++)
                        for(int r=j; r<=8; r++)
                        {
                            for(int t=i+1; t<=h; t++)
                                dp[i][j][h][r][c]=min(dp[i][j][h][r][c],(min(dp[i][j][t-1][r][c-1]+dp[t][j][h][r][1],dp[i][j][t-1][r][1]+dp[t][j][h][r][c-1])));
                            for(int t=j+1; t<=r; t++)
                                dp[i][j][h][r][c]=min(dp[i][j][h][r][c],(min(dp[i][j][h][t-1][c-1]+dp[i][t][h][r][1],dp[i][j][h][t-1][1]+dp[i][t][h][r][c-1])));
                        }
        double a=s[8][8]*1.0/n;
        printf("%.3f\n",sqrt(dp[1][1][8][8][n]*1.0/n-a*a));
    }
    return 0;
}