【题解】HAOI2007分割矩阵

　　水题盛宴啦啦啦……做起来真的极其舒服，比某些毒瘤题好太多了……

　　数据范围极小 --> 状压 / 搜索 / 高维度dp；观察要求的均方差，开始考虑是不是能够换一下式子。我们用\(a_{x}\)来表示第 \(x\) 个矩阵的总值，则式子为：

　　\(ans = sqrt \frac{{\left ( \sum_{1}^{n} a_{x} - \bar{x} \right )^2}}{n}\)

　　转化一下，化成：

　　\(ans = sqrt \frac{{\left ( -n\bar{x}^2 + \sum_{1}^{n}a_{x}^2 \right )}}{n}\)

　　然后问题就变成了：使划分出来的矩阵的平方和最小。直接上记忆化搜索，BINGO~

#include <bits/stdc++.h>
using namespace std;
#define maxn 11
#define INF 99999999
#define db double
int n, m, K, sum[maxn][maxn];
int f[maxn][maxn][maxn][maxn][maxn];
int tot;
db aver;

int read()
{
    int x = ;
    char c;
    c = getchar();
    while(c < '' || c > '') c = getchar();
    while(c >= '' && c <= '') x = x *  + c - '', c = getchar();
    return x;
}

int dfs(int a, int b, int c, int d, int K)
{
    int *F = f[a][b][c][d];
    if(~F[K]) return F[K];
    else F[K] = INF;
    if(c * d - (a - ) * (b - ) < K) return F[K] = INF;
    if(K == )
    {
        F[K] = sum[c][d] - sum[a - ][d] - sum[c][b - ] + sum[a - ][b - ];
        return F[K] = F[K] * F[K];
    }
    for(int i = a; i < c; i ++)
        for(int j = ; j < K; j ++)
        {
            int x1 = dfs(a, b, i, d, j);
            int x2 = dfs(i + , b, c, d, K - j);
            F[K] = min(F[K], x1 + x2);
        }
    for(int i = b; i < d; i ++)
        for(int j = ; j < K; j ++)
        {
            int x1 = dfs(a, b, c, i, j);
            int x2 = dfs(a, i + , c, d, K - j);
            F[K] = min(F[K], x1 + x2);
        }
    return F[K];
}

int main()
{
    n = read(), m = read(), K = read();
    memset(f, -, sizeof(f));
    for(int i = ; i <= n; i ++)
        for(int j = ; j <= m; j ++)
        {
            int x = read(); tot += x;
            sum[i][j] = sum[i - ][j] + sum[i][j - ] - sum[i - ][j - ] + x;
        }
    db aver = (db) tot / K;
    dfs(, , n, m, K);
    printf("%.2f\n", sqrt(((db) f[][][n][m][K] - aver * aver * (db) (K)) / (db) K));
    return ;
}