列表

详情


NC19966. [HAOI2007]分割矩阵

描述

将一个a*b的数字矩阵进行如下分割:将原矩阵沿某一条直线分割成两个矩阵,再将生成的两个矩阵继续如此 分割(当然也可以只分割其中的一个),这样分割了(n-1)次后,原矩阵被分割成了n个矩阵。(每次分割都只能 沿着数字间的缝隙进行)原矩阵中每一位置上有一个分值,一个矩阵的总分为其所含各位置上分值之和。现在需要把矩阵按上述规则分割成n个矩阵,并使各矩阵总分的均方差最小。请编程对给出的矩阵及n,求出均方差的最小值 。

输入描述

第一行为3个整数,表示a,b,n(1 < a,b ≤ 10,1 < n ≤ 10)的值。
第二行至第n+1行每行为b个小于100的非负整数,表示矩阵中相应位置上的分值。
每行相邻两数之间用一个空格分开。

输出描述

仅一个数,为均方差的最小值(四舍五入精确到小数点后2位)

示例1

输入:

5 4 4
2 3 4 6
5 7 5 1
10 4 0 5
2 0 2 3
4 1 1 1

输出:

0.50

原站题解

上次编辑到这里,代码来自缓存 点击恢复默认模板

C++(clang++ 11.0.1) 解法, 执行用时: 12ms, 内存消耗: 1908K, 提交时间: 2022-11-21 23:01:47

#include <bits/stdc++.h>
#define sqr(x) ((x)*(x))
using namespace std;
const int N =  17;
int n, m, k;
int a[N][N], s[N][N];
double f[N][N][N][N][N], ave;

double sum(int a, int b, int c, int d){
    return double(s[c][d] - s[a-1][d] - s[c][b-1] + s[a-1][b-1]);
}

double dp(int a, int b, int c, int d, int num){
    if(f[a][b][c][d][num]) return f[a][b][c][d][num];
    else if(num == 1) return sqr(sum(a, b, c, d) - ave);
    f[a][b][c][d][num] = 1e9;
    for(int i = b; i < d; i++)
        for(int j = 1; j < num; j++){
            f[a][b][c][d][num] = min(f[a][b][c][d][num], 
                   dp(a, b, c, i, j) + dp(a, i + 1, c, d, num - j));
        }
    for(int i = a; i < c; i++)
        for(int j = 1; j < num; j++){
            f[a][b][c][d][num]=min(f[a][b][c][d][num],
                    dp(a, b, i, d, j)+dp(i+1, b, c, d, num - j));
        }
    return f[a][b][c][d][num];
}

int main()
{
    scanf("%d%d%d", &n, &m, &k);
    for(int i = 1; i <= n; i++)
        for(int j = 1; j <= m; j++){
            scanf("%d", &a[i][j]);
            s[i][j] += s[i - 1][j] + s[i][j - 1] - s[i - 1][j - 1] + a[i][j];
        }
    ave = (double) s[n][m] / (double)k;
    printf("%.2lf", sqrt(dp(1, 1, n, m, k) / (double)k));
    
}

上一题