Tips
1、每一次分割都有四种方法,留上,留下,留左,留右
2、f(k,棋盘)={f(1,留下的棋盘)+f(k-1,待分割的棋盘)} k>=2
3、对于每次fun(n,x1,y1,x2,y2)来说,可能多次使用这个值,每次都重新计算太花时间,所以使用记录表,用res[n][x1][y1][x2][y2]来记录fun(n,x1,y1,x2,y2)
4、使用setiosflags(ios::fixed)<<setprecision(3)确定精度
#include<iostream>
#include<stdio.h>
#include<iomanip>
#include<math.h>
using namespace std;
int s[9][9];
int sum[9][9];
int res[15][9][9][9][9];
int calSum(int x1,int y1,int x2,int y2)
{
return sum[x2][y2]+sum[x1-1][y1-1]-sum[x1-1][y2]-sum[x2][y1-1];
}
int fun(int n,int x1,int y1,int x2,int y2)
{
int t,a,b,c,e,MIN=10000000;
if(res[n][x1][y1][x2][y2]!=-1)
return res[n][x1][y1][x2][y2];
if (n==1) //迭代终止条件
{
t=calSum(x1,y1,x2,y2);
res[n][x1][y1][x2][y2]=t*t;
return t*t;
}
for(a=x1;a<x2;a++)
{
c=calSum(a+1,y1,x2,y2); //右边的块
e=calSum(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=calSum(x1,b+1,x2,y2);//下边的块
e=calSum(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;
}