棋盘分割
Time Limit: 1000MS | Memory Limit: 10000K | |
Total Submissions: 14692 | Accepted: 5248 |
Description
将一个8*8的棋盘进行如下分割:将原棋盘割下一块矩形棋盘并使剩下部分也是矩形,再将剩下的部分继续如此分割,这样割了(n-1)次后,连同最后剩下的矩形棋盘共有n块矩形棋盘。(每次切割都只能沿着棋盘格子的边进行)
原棋盘上每一格有一个分值,一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘,并使各矩形棋盘总分的均方差最小。
均方差 ,其中平均值 ,x i为第i块矩形棋盘的总分。
请编程对给出的棋盘及n,求出O'的最小值。
原棋盘上每一格有一个分值,一块矩形棋盘的总分为其所含各格分值之和。现在需要把棋盘按上述规则分割成n块矩形棋盘,并使各矩形棋盘总分的均方差最小。
均方差 ,其中平均值 ,x i为第i块矩形棋盘的总分。
请编程对给出的棋盘及n,求出O'的最小值。
Input
第1行为一个整数n(1 < n < 15)。
第2行至第9行每行为8个小于100的非负整数,表示棋盘上相应格子的分值。每行相邻两数之间用一个空格分隔。
第2行至第9行每行为8个小于100的非负整数,表示棋盘上相应格子的分值。每行相邻两数之间用一个空格分隔。
Output
仅一个数,为O'(四舍五入精确到小数点后三位)。
Sample Input
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
Sample Output
1.633
参考:http://wenku.baidu.com/link?url=Omo0iEun_QSOQRA6vsSMii65boarphsLPNt3mBpJgq-8nXhH6F01FcHormzS3qVXFREP7e0ad91zaimzKbEnhxq-m-szyOEo160NpuLCFVm
#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
using namespace std;
#define INF 0x3f3f3f3f
int sum[10][10];
int s[10][10];
int dp[20][10][10][10][10];
int n;
int calsum(int x1, int y1, int x2, int y2)
{
return sum[x2][y2] - sum[x2][y1 - 1] -sum[x1 - 1][y2] + sum[x1-1][y1-1];
}
int dfs( int k, int x1, int y1, int x2, int y2 )
{
if(dp[k][x1][y1][x2][y2] != -1)
return dp[k][x1][y1][x2][y2];
if(k == 1)
{
int t = calsum(x1, y1, x2, y2);
dp[k][x1][y1][x2][y2] = t * t; return dp[k][x1][y1][x2][y2];
}
int Min = INF;
for( int i = x1; i < x2; i++ )
{
int t = calsum(i + 1, y1, x2, y2);
int e = calsum( x1, y1, i, y2);
int up = dfs(k-1, x1, y1, i, y2) + t * t ;
int down = dfs(k-1, i + 1, y1, x2, y2) + e * e;
int aa = min(up, down);
if(aa < Min) Min = aa;
}
for( int i = y1; i < y2; i++)
{
int t = calsum(x1, i + 1, x2, y2);
int e = calsum(x1, y1, x2, i);
int left = dfs(k - 1, x1, y1, x2, i) + t*t;
int right = dfs(k - 1, x1, i + 1, x2, y2) + e * e;
int aa = min(left, right);
if(Min > aa) Min = aa;
}
dp[k][x1][y1][x2][y2] = Min;
return dp[k][x1][y1][x2][y2];
}
int main()
{
scanf("%d", &n);
memset(sum, 0, sizeof(sum));
memset(dp, -1, sizeof(dp));
for( int i = 1; i <= 8; i++)
{
int rowsum = 0;
for( int j = 1; j <= 8; j++)
{
scanf("%d", &s[i][j]);
rowsum += s[i][j];
sum[i][j] += sum[i-1][j] + rowsum;
}
}
double aver = 1.0*sum[8][8]/n;
double ss = 1.0*dfs(n, 1, 1, 8, 8)/n;
double last = sqrt( ss - aver * aver);
printf("%.3lf\n",last);
}