根据这2个公式可以 得到 O^2 = sum(x1^2 + x2^2 + x3^2 +...xi^2)/n - 平均值^2
所以我们就是要使得总分的平方和尽量小……
对于一次切割,可以横着切,可以竖着切,可以使剩下了的包含左上角,也可以使剩下了的包含右上角。
具体方程什么的看代码吧。
G++提交的时候记得把.3lf改成.3f
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
double sum[9][9];
double a[9][9];
int n;
double ave;
double d[16][9][9][9][9];
double s(int x1,int y1,int x2,int y2)
{
double tmp=sum[x2][y2]-sum[x1-1][y2]-sum[x2][y1-1]+sum[x1-1][y1-1];
return tmp*tmp;
}
double dfs(int k,int x1,int y1,int x2,int y2)
{
if(k==1)
{
d[k][x1][y1][x2][y2]=s(x1,y1,x2,y2);
return d[k][x1][y1][x2][y2];
}
if(d[k][x1][y1][x2][y2]) return d[k][x1][y1][x2][y2];
double t1=100000000,t2=100000000;
for(int a=x1;a<x2;a++)
{
t1=min(t1,min(dfs(k-1,x1,y1,a,y2)+s(a+1,y1,x2,y2),dfs(k-1,a+1,y1,x2,y2)+s(x1,y1,a,y2)));
}
for(int b=y1;b<y2;b++)
{
t2=min(t2,min(dfs(k-1,x1,y1,x2,b)+s(x1,b+1,x2,y2),dfs(k-1,x1,b+1,x2,y2)+s(x1,y1,x2,b)));
}
d[k][x1][y1][x2][y2]=min(t1,t2);
return d[k][x1][y1][x2][y2];
}
int main()
{
while(scanf("%d",&n)!=EOF)
{
memset(sum,0,sizeof(sum));
memset(d,0,sizeof(d));
for(int i=1;i<=8;i++)
for(int j=1;j<=8;j++)
{
scanf("%lf",&a[i][j]);
sum[i][j]=a[i][j]+sum[i-1][j]+sum[i][j-1]-sum[i-1][j-1];
}
ave=sum[8][8]/n;
dfs(n,1,1,8,8);
printf("%.3f\n",sqrt(d[n][1][1][8][8]/n-ave*ave));
}
return 0;
}