棋盘分割

将一个８*８的棋盘进行如下分割：将原棋盘割下一块矩形棋盘并使剩下部分也是矩形，再将剩下的部分继续如此分割，这样割了(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

比赛的时候，没想出来怎么做这道题，后来的时候发现用记忆化搜索可以做这道题，但是重复度蛮高的。就是把现在的棋盘分为两份然后求方差一次，然后再分，再求方差，一直重复这个操作，直到n-1次，看看怎么样的方差和最小。这样求每次先把割完的四个角的坐标记录下来，然后求出面积，然后去计算这个面积平均分到每个小块上面是多少，让这个平均的面积最小就完成了我们的目的。

  dp[i1][j1][i2][j2][i]=min(dp[i1][j1][i2][j2][i], min(dp[i1][j1][i2][k][i-1]+dp[i1][k+1][i2][j2][0], dp[i1][j1][i2][k][0]+dp[i1][k+1][i2][j2][i-1]));

  dp[i1][j1][i2][j2][i]=min(dp[i1][j1][i2][j2][i], min(dp[i1][j1][k][j2][i-1]+dp[k+1][j1][i2][j2][0], dp[i1][j1][k][j2][0]+dp[k+1][j1][i2][j2][i-1]));

上面这两个就是动态转移方程。

#include<stdio.h>
#include<math.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define inf 0x3f3f3f3f
int n,a[9][9]={0},book[9][9]={0},dp[9][9][9][9][15]={0};
int cn(int y1,int x1,int y2,int x2)
{
    int u=book[y2][x2]-book[y2][x1-1]-book[y1-1][x2]+book[y1-1][x1-1];
    return u*u;
}
int main()
{
    scanf("%d", &n);
    for(int i=1; i<=8; i++)
        for(int j=1; j<=8; j++)
            scanf("%d",&a[i][j]);
    for(int i=1; i<=8; i++)
    {
        for(int j=1; j<=8; j++)
            book[i][j]=book[i][j-1]+a[i][j];
        for(int j=1; j<=8; j++)
            book[i][j]+=book[i-1][j];
    }
    for(int i1=1; i1<=8; i1++)
        for(int j1=1; j1<=8; j1++)
            for(int i2=i1; i2<=8; i2++)
                for(int j2=j1; j2<=8; j2++)
                    dp[i1][j1][i2][j2][0]=cn(i1,j1,i2,j2);
    for(int i=1; i<n; i++)
        for(int i1=1; i1<=8; i1++)
            for(int j1=1; j1<=8; j1++)
                for(int i2=i1; i2<=8; i2++)
                    for(int j2=j1; j2<=8; j2++)
                    {
                        dp[i1][j1][i2][j2][i]=inf;
                        for(int k=j1; k<j2; k++)
                            dp[i1][j1][i2][j2][i]=min(dp[i1][j1][i2][j2][i], min(dp[i1][j1][i2][k][i-1]+dp[i1][k+1][i2][j2][0], dp[i1][j1][i2][k][0]+dp[i1][k+1][i2][j2][i-1]));
                        for(int k=i1; k<i2; k++)
                            dp[i1][j1][i2][j2][i]=min(dp[i1][j1][i2][j2][i], min(dp[i1][j1][k][j2][i-1]+dp[k+1][j1][i2][j2][0], dp[i1][j1][k][j2][0]+dp[k+1][j1][i2][j2][i-1]));
                    }
    printf("%.3f\n", sqrt(double(dp[1][1][8][8][n-1])/n-double(book[8][8]*book[8][8])/n/n));
    return 0;
}