BLMOOC1191棋盘分割

总时间限制:

1000ms

内存限制:

65536kB

描述

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

原棋盘上每一格有一个分值，一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘，并使各矩形棋盘总分的均方差最小。
均方差，其中平均值，xi为第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

#include<bits/stdc++.h>
using namespace std;
int s[9][9];//每个格子的分数
int sum [9][9];//(1,1)到(i,j)的矩形分数之和
int res[15][9][9][9][9];//fun函数的记录表

int calsum(int x1,int y1,int x2,int y2){//计算(x1,y1)到(x2,y2)的矩形的分数之和
    return sum[x2][y2]-sum[x2][y1-1]-sum[x1-1][y2]+sum[x1-1][y1-1];
}

int fun(int n,int x1,int y1,int x2,int y2){
    int temp,i,j,cut,remain,MINN=10000000;
    if(res[n][x1][y1][x2][y2]!=-1)
        return res[n][x1][y1][x2][y2];
    if(n==1){
        temp=calsum(x1,y1,x2,y2);
        res[n][x1][y1][x2][y2]=temp*temp;
        return temp*temp;
    }
    for(i=x1;i<x2;i++){//沿x方向切
        cut=calsum(i+1,y1,x2,y2);
        remain=calsum(x1,y1,i,y2);
        temp=min(fun(n-1,x1,y1,i,y2)+cut*cut,
                 fun(n-1,i+1,y1,x2,y2)+remain*remain);
        if(MINN>temp)
            MINN=temp;
    }
    for(j=y1;j<y2;j++){//沿y方向切
        cut=calsum(x1,j+1,x2,y2);
        remain=calsum(x1,y1,x2,j);
        temp=min(fun(n-1,x1,y1,x2,j)+cut*cut,
                 fun(n-1,x1,j+1,x2,y2)+remain*remain);
        if(MINN>temp)
            MINN=temp;
    }
    res[n][x1][y1][x2][y2]=MINN;
    return MINN;
}

int main(){
    memset(sum,0,sizeof(sum));
    memset(res,-1,sizeof(res));
    int n;
    scanf("%d",&n);
    for(int i=1;i<9;i++)
    for(int j=1,rowsum=0;j<9;j++){
        scanf("%d",&s[i][j]);
        rowsum+=s[i][j];
        sum[i][j]+=sum[i-1][j]+rowsum;
    }
    double result=n*fun(n,1,1,8,8)-sum[8][8]*sum[8][8];
    cout<<setiosflags(ios::fixed)<<setprecision(3)<<sqrt(result/(n*n))<<endl;
    return 0;
}

最后输出也可以：

 printf("%.3f",sqrt(result/(n*n)));