Description
将一个a*b的数字矩阵进行如下分割:将原矩阵沿某一条直线分割成两个矩阵,再将生成的两个矩阵继续如此
分割(当然也可以只分割其中的一个),这样分割了(n-1)次后,原矩阵被分割成了n个矩阵。(每次分割都只能
沿着数字间的缝隙进行)原矩阵中每一位置上有一个分值,一个矩阵的总分为其所含各位置上分值之和。现在需要
把矩阵按上述规则分割成n个矩阵,并使各矩阵总分的均方差最小。请编程对给出的矩阵及n,求出均方差的最小值
。
Input
第一行为3个整数,表示a,b,n(1<a,b<=10,1<n<=10)的值。
第二行至第n+1行每行为b个小于100的非负整数,表示矩阵中相应位置上的分值。每行相邻两数之间用一个空
格分开。
Output
仅一个数,为均方差的最小值(四舍五入精确到小数点后2位)
Sample Input
5 4 4
2 3 4 6
5 7 5 1
10 4 0 5
2 0 2 3
4 1 1 1
2 3 4 6
5 7 5 1
10 4 0 5
2 0 2 3
4 1 1 1
Sample Output
0.50
HINT
直接记忆化搜索乱搞
#include <stdio.h>
#include <math.h>
#include <algorithm>
using namespace std;
const int maxn=17;
int q[maxn][maxn];
double sum[maxn][maxn];
double sum1[maxn][maxn];
int r,c;
double mian(int a,int b,int c,int d)
{
return sum[a][b]-sum[a][d-1]-sum[c-1][b]+sum[c-1][d-1];
}
double dp[11][11][11][11][11];
double ave;
double dfs(int a,int b,int c,int d,int e)
{
double& r=dp[a][b][c][d][e];
if(r!=-1) return r;
double y=mian(c,d,a,b);
if(e==1)
{
return r=(y-ave)*(y-ave);
}
r=999999999.0;
for(int i=a;i<c;i++)
for(int k=1;k<e;k++)
r=min(r,dfs(a,b,i,d,k)+dfs(i+1,b,c,d,e-k));
for(int i=b;i<d;i++)
for(int k=1;k<e;k++)
r=min(r,dfs(a,b,c,i,k)+dfs(a,i+1,c,d,e-k));
return r;
}
int main(){
int n;
scanf("%d%d%d",&r,&c,&n);
for(int i=1;i<=r;i++)
for(int k=1;k<=c;k++)
{
scanf("%d",&q[i][k]);
}
for(int i=1;i<=r;i++)
for(int k=1;k<=c;k++)
sum1[i][k]=sum1[i][k-1]+q[i][k];
for(int i=1;i<=r;i++)
for(int k=1;k<=c;k++)
for(int j=i;j<=r;j++)
for(int z=k;z<=c;z++)
for(int p=1;p<=n;p++)
dp[i][k][j][z][p]=-1;
for(int i=1;i<=r;i++)
{
for(int k=1;k<=c;k++)
{
sum[i][k]=sum[i-1][k]+sum1[i][k];
}
}
ave=(double)sum[r][c]/n;
printf("%.2lf",sqrt(dfs(1,1,r,c,n)/n));
return 0;
}