棋盘分割
时间限制:
1000 ms | 内存限制:
65535 KB
难度:
6
-
描述
-
将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的部分继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共有n块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行)
原棋盘上每一格有一个分值,一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘,并使各矩形棋盘总分的均方差最小。
均方差 ,其中平均值 ,x i为第i块矩形棋盘的总分。
请编程对给出的棋盘及n,求出O'的最小值。
-
输入
-
第1行为一个整数n(1 < n < 15)。
第2行至第9行每行为8个小于100的非负整数,表示棋盘上相应格子的分值。每行相邻两数之间用一个空格分隔。
输出
-
仅一个数,为O'(四舍五入精确到小数点后三位)。
样例输入
-
3 1 1 1 1 1 1 1 3 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 1 1 0 3
样例输出
-
1.633
来源
- IOI 2000 上传者
一开始理解错题意了,以为切完后所有块都要操作,后来才发现一分为二时,其中一块便不能再动了,呢就是由大化小的问题了,动态规划随便搞搞,感觉记忆化搜索会好写一点。。。。。
每次四种选择,分别搜索:
(1)横着切,选上边一块
(2)横着切,选下边一块
(3)竖着切,选左边一块
(4)竖着切,选右边一块
剩下的就是dp了,定义dp[kk][i][j][k][h]:左上角(i,j)到右下角(k,h)组成的矩形,切到第kk次的最优解。
#include<map> #include<stack> #include<queue> #include<vector> #include<math.h> #include<stdio.h> #include<iostream> #include<string.h> #include<stdlib.h> #include<algorithm> #include<functional> using namespace std; typedef long long ll; #define inf 1000000000 #define MOD 1000000007 #define maxn 5100005 #define lowbit(x) (x&-x) #define eps 1e-9 int dp[20][10][10][10][10],a[10][10][10][10],b[10][10]; int sm(int x,int y,int xx,int yy) { int ans=0,i,j; for(i=x;i<=xx;i++) for(j=y;j<=yy;j++) ans+=b[i][j]; return ans; } int dfs(int k,int x,int y,int xx,int yy) { if(dp[k][x][y][xx][yy]!=-1) return dp[k][x][y][xx][yy]; int i,j,t1,t2,t; dp[k][x][y][xx][yy]=inf; if(xx>x) { for(i=x;i<xx;i++) { t1=dfs(k-1,i+1,y,xx,yy); t2=dfs(k-1,x,y,i,yy); t=min(t1+a[x][y][i][yy],t2+a[i+1][y][xx][yy]); dp[k][x][y][xx][yy]=min(dp[k][x][y][xx][yy],t); } } if(yy>y) { for(i=y;i<yy;i++) { t1=dfs(k-1,x,i+1,xx,yy); t2=dfs(k-1,x,y,xx,i); t=min(t1+a[x][y][xx][i],t2+a[x][i+1][xx][yy]); dp[k][x][y][xx][yy]=min(dp[k][x][y][xx][yy],t); } } return dp[k][x][y][xx][yy]; } int main(void) { int x,y,xx,yy,n,i,j; scanf("%d",&n); for(i=1;i<=8;i++) for(j=1;j<=8;j++) scanf("%d",&b[i][j]); memset(dp,-1,sizeof(dp)); for(x=1;x<=8;x++) for(xx=x;xx<=8;xx++) for(y=1;y<=8;y++) for(yy=y;yy<=8;yy++) { int tmp=sm(x,y,xx,yy); dp[1][x][y][xx][yy]=a[x][y][xx][yy]=tmp*tmp; } dfs(n,1,1,8,8); double ans,sum; sum=1.0*sm(1,1,8,8); sum=(sum/n)*(sum/n); ans=sqrt(1.0*dp[n][1][1][8][8]/n-sum); printf("%.3f\n",ans); return 0; }
-
第1行为一个整数n(1 < n < 15)。