二、递归(二)
2017-08-10 本文已影响12人
安东可
一、棋盘分割问题

【每次切完一个棋盘,就是需要舍弃掉一面,只能在一面操作】











#include "stdafx.h"
#include<iostream>
#include<algorithm>
#include<iomanip>
using namespace std;
int s[9][9]; //每个格子的分数
int sum[9][9]; //从(1,1)到(i.j)的分数和
int res[15][9][9][9][9];//记录表
//计算从(x1,y1)到(x2,y2)的矩形分数和
int calcSum(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];
}
//分成n分的最小平方和
int fun(int n, int x1, int y1, int x2, int y2){
int t, a, b, c, e, MIN = 1000000;
if (res[n][x1][y1][x2][y2] != -1){
return res[n][x1][y1][x2][y2];
}
if (n == 1){
t = calcSum(x1, y1, x2, y2);
res[n][x1][y1][x2][y2] = t* t;
}
//分成左右两块,选择较小的一块
for (a = x1; a<x2; a++) {
c = calcSum(a + 1, y1, x2, y2);
e = calcSum(x1, y1, a, y2);
t = min(fun(n - 1, x1, y1, a, y2) + c*c, fun(n - 1, a + 1, y1, x2, y2) + e*e);
if (MIN>t) MIN = t;
}
//分成上下两块
for (b = y1; b<y2; b++) {
c = calcSum(x1, b + 1, x2, y2);
e = calcSum(x1, y1, x2, b);
t = min(fun(n - 1, x1, y1, x2, b) + c*c, fun(n - 1, x1, b + 1, x2, y2) + e*e);
if (MIN>t) MIN = t;
}
res[n][x1][y1][x2][y2] = MIN;
return MIN;
}
int main() {
memset(sum, 0, sizeof(sum));
memset(res, -1, sizeof(res)); //初始化记录表
int n;
cin >> n;
for (int i = 1; i<9; i++)
for (int j = 1, rowsum = 0; j<9; j++) {
cin >> s[i][j];
rowsum += s[i][j];
sum[i][j] += sum[i - 1][j] + rowsum;
}
double result = n*fun(n, 1, 1, 8, 8) - sum[8][8] * sum[8][8];
cout << setiosflags(ios::fixed) << setprecision(3) << sqrt(result / (n*n)) << endl;
return 0;
}