问题 x: 棋盘分割----区间dp

第二天叫醒我的不是闹钟，是梦想！

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

原棋盘上每一格有一个分值，一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘，并使各矩形棋盘总分的均方差最小。
均方差，其中平均值，xi为第i块矩形棋盘的总分。
请编程对给出的棋盘及n，求出O’的最小值。
输入
第1行为一个整数n(1 < n < 15)。
第2行至第9行每行为8个小于100的非负整数，表示棋盘上相应格子的分值。每行相邻两数之间用一个空格分隔。
输出
仅一个数，为O’（四舍五入精确到小数点后三位）。
样例输入 Copy
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
样例输出 Copy
1.633

//二维差分+区间dp。分为横切和竖切

#include<bits/stdc++.h>
using namespace std;
const double INF=1e9;
const int N=15,M=9;
int n;
double x;
int m=8;
double f[M][M][M][M][N];
int s[M][M];
double get(int x1,int y1,int x2,int y2)
{
    double sum=s[x2][y2]-s[x2][y1-1]-s[x1-1][y2]+s[x1-1][y1-1]-x;
    return (double)(sum*sum)/n;
}
double dp(int x1,int y1,int x2,int y2,int k)
{
    double &v=f[x1][y1][x2][y2][k];
    if(v>=0) return v;
    if(k==1) return v=get(x1,y1,x2,y2);
    v=INF;
    for(int i=x1;i<x2;i++)
    {
        v=min(v,dp(x1,y1,i,y2,k-1)+get(i+1,y1,x2,y2));
        v=min(v,dp(i+1,y1,x2,y2,k-1)+get(x1,y1,i,y2));
    }
    for(int i=y1;i<y2;i++)
    {
        v=min(v,dp(x1,y1,x2,i,k-1)+get(x1,i+1,x2,y2));
        v=min(v,dp(x1,i+1,x2,y2,k-1)+get(x1,y1,x2,i));
    }
    return v;
}
int main()
{
    cin>>n;
    for(int i=1;i<=m;i++)
        for(int j=1;j<=m;j++)
        {
            cin>>s[i][j];
            s[i][j]+=s[i-1][j]+s[i][j-1]-s[i-1][j-1];
        }
        memset(f,-1,sizeof f);
     x=(double ) s[m][m]/n;
    printf("%.3lf\n",sqrt(dp(1,1,8,8,n)));
}