POJ - 1191 棋盘分割

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

题解：每个矩形用最左上和最右下两个点表示，比如最初的矩形表示为（1，1，8，8），dp[x1][y1][x2][y2][n]，n表示切n块。注意切块只能横切或竖切并且切到底，这样才能是矩形（头疼）

看的这个https://blog.csdn.net/v5zsq/article/details/46580505懂大概。。。。

代码：

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#define inf 0x3f3f3f3f
using namespace std;
int dp[10][10][10][10][20];
int mp[10][10],s[10][10];
int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        int i,j,k,l;
        for(i=1; i<=8; i++)
            for(j=1; j<=8; j++)
                scanf("%d",&mp[i][j]);
        memset(dp,inf,sizeof(dp));
        memset(s,0,sizeof(s));
        for(i=1; i<=8; i++)
            for(j=1; j<=8; j++)
                s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+mp[i][j];//表示（1，1）到（i，j）这个矩形的值
        for(i=1; i<=8; i++)
            for(j=1; j<=8; j++)
                for(k=1; k<=8; k++)
                    for(l=1; l<=8; l++)
                        dp[i][j][k][l][1]=(s[k][l]-s[k][j-1]-s[i-1][l]+s[i-1][j-1])*(s[k][l]-s[k][j-1]-s[i-1][l]+s[i-1][j-1]);
        for(int num=2; num<=n; num++)//切几块
            for(i=1; i<=8; i++)
                for(j=1; j<=8; j++)
                    for(k=i; k<=8; k++)
                        for(l=j; l<=8; l++)
                        {
                            for(int o=i+1; o<=k; o++)//横切
                                dp[i][j][k][l][num]=min(dp[i][j][k][l][num],(min(dp[i][j][o-1][l][num-1]+dp[o][j][k][l][1],dp[i][j][o-1][l][1]+dp[o][j][k][l][num-1])));
                            for(int o=j+1; o<=l; o++)//竖切
                                dp[i][j][k][l][num]=min(dp[i][j][k][l][num],(min(dp[i][j][k][o-1][num-1]+dp[i][o][k][l][1],dp[i][j][k][o-1][1]+dp[i][o][k][l][num-1])));
                        }
        double a=dp[1][1][8][8][n],b=s[8][8]*1.0/n;
        printf("%.3f\n",sqrt(a*1.0/n-b*b));
    }
    return 0;
}