BZOJ 1048 [HAOI2007]分割矩阵 记忆化搜索

数据范围果断小，果断暴力。。。

 

 

View Code 
 1 #include <iostream>
 2 #include <cstring>
 3 #include <cstdio>
 4 #include <cstdlib>
 5 #include <algorithm>
 6 #include <cmath>
 7 
 8 #define N 14
 9 #define INF 1e9
10 
11 using namespace std;
12 
13 int a,b,n;
14 double sm[N][N],map[N][N],dp[N][N][N][N][N],sum,sumb;
15 
16 inline void read()
17 {
18     scanf("%d%d%d",&a,&b,&n);
19     for(int i=1;i<=a;i++)
20         for(int j=1;j<=b;j++)
21             scanf("%lf",&map[i][j]),sum+=map[i][j];
22     sumb=sum/n;
23     for(int i=1;i<=a;i++)
24         for(int j=1;j<=b;j++)
25             sm[i][j]=sm[i-1][j]+sm[i][j-1]-sm[i-1][j-1]+map[i][j];
26 }
27 
28 inline double dfs(int x1,int y1,int x2,int y2,int cs)
29 {
30     double &res=dp[x1][y1][x2][y2][cs];
31     if(res<INF) return res;
32     if(cs==0)
33     {
34         res=(sm[x2][y2]+sm[x1-1][y1-1]-sm[x1-1][y2]-sm[x2][y1-1])-sumb;
35         res=res*res;
36         return res;
37     }
38     for(int i=x1;i<x2;i++)
39         for(int j=0;j<=cs-1;j++)
40             res=min(res,dfs(x1,y1,i,y2,j)+dfs(i+1,y1,x2,y2,cs-j-1));
41     for(int i=y1;i<y2;i++)
42         for(int j=0;j<=cs-1;j++)
43             res=min(res,dfs(x1,y1,x2,i,j)+dfs(x1,i+1,x2,y2,cs-j-1));
44     return res;
45 }
46 
47 inline void go()
48 {
49     for(int i=0;i<=10;i++)
50         for(int j=0;j<=10;j++)
51             for(int k=0;k<=10;k++)
52                 for(int p=0;p<=10;p++)
53                     for(int q=0;q<=10;q++)
54                         dp[i][j][k][p][q]=INF;
55     dfs(1,1,a,b,n-1);
56     printf("%.2lf\n",sqrt(dp[1][1][a][b][n-1]/n));
57 }
58 
59 int main()
60 {
61     read(),go();
62     return 0;
63 }

 

 

posted @ 2013-03-05 21:34  proverbs  阅读(553)  评论(0)  编辑  收藏  举报