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

Sample Output

HINT

Source

我真是弱的可以了0.0 一個記憶化搜索搞了這麽長時間0.0我可是在yy那裏立了軍令狀的啊QAQ

 1 #include "
 
bits/stdc++.h"
 2 using namespace std;
 3 typedef long long LL;
 4 const int MAX=15;
 5 int n,m,t;
 6 double f[MAX][MAX][MAX][MAX][MAX];
 7 int s[MAX][MAX];
 8 double ave;
 9 inline int read(){
10     int x=1,an=0;char c=getchar();
11     while (c<‘0‘ || c>‘9‘) {if (c==‘-‘) x=-1;c=getchar();}
12     while
 
 (c>=‘0‘ && c<=‘9‘) {an=an*10+c-‘0‘;c=getchar();}
13     return an*x;
14 }
15 inline double mn(double x,double y){return x<y?x:y;}
16 double dfs(int a,int b,int c,int d,int x){
17     if (f[a][b][c][d][x]!=-1) return f[a][b][c][d][x];
18     if (x==1){
19         double zt;
20         zt=(s[c][d]-s[a-1][d]-s[c][b-1]+s[a-1][b-1])*1.0;
21         return f[a][b][c][d][x]=(zt-ave)*(zt-ave);
22     }
23     int i,j;double tmp=100000000.0;
24     for (i=a;i<c;i++){
25         for (j=1;j<x;j++){
26             tmp=mn(tmp,dfs(a,b,i,d,j)+dfs(i+1,b,c,d,x-j));
27         }
28     }
29     for (i=b;i<d;i++){
30         for (j=1;j<x;j++){
31             tmp=mn(tmp,dfs(a,b,c,i,j)+dfs(a,i+1,c,d,x-j));
32         }
33     }
34     return f[a][b][c][d][x]=tmp;
35 }
36 int main(){
37     freopen ("matrix.in","r",stdin);
38     freopen ("matrix.out","w",stdout);
39     int i,j;
40     n=read();m=read();t=read();
41     memset(s,0,sizeof(s));
42     for (int i1=0;i1<=n+1;i1++) for (int i2=0;i2<=m+1;i2++) for (int i3=0;i3<=n+1;i3++) for (int i4=0;i4<=m+1;i4++) for (int i5=1;i5<=t;i5++) f[i1][i2][i3][i4][i5]=-1;
43     for (i=1;i<=n;i++){
44         for (j=1;j<=m;j++){
45             s[i][j]=read();
46             s[i][j]=s[i][j]+s[i-1][j]+s[i][j-1]-s[i-1][j-1];
47         }
48     }
49     ave=s[n][m]*1.0/t;
50     printf("%.2lf",sqrt(dfs(1,1,n,m,t)/(t*1.0)));
51     return 0;
52 }

BZOJ-1048: [HAOI2007]分割矩陣 (記憶化搜索)