那个均方差,可以通过展开、合并Σ,发现最终只有Xi^2会对答案造成影响,其他都是定值,所以求出最小的和的平方就行。
其实这才是这题最难的部分,以下都是码农部分。
f[x1][y1][x2][y2][k] 从(x1,y1)到(x2,y2)这个矩形,割k次的最小平方和。
边界 k==1时,就是矩形本身面积的平方和。
其余情况,分两种,横着割和竖着割,一块继续割,另一块的面积平方加入答案。
(强烈建议纸笔画图标出坐标,不用纠结下标+1/-1问题)
维数比较多,所以用了记忆化,好写一点。
洛谷AC代码
//Stay foolish,stay hungry,stay young,stay simple
#include<iostream>
using namespace std;
int n;
int a[9][9],s[9][9];
int f[9][9][9][9][16];
inline int S(int x1,int y1,int x2,int y2){
x1--;y1--;//it includes (x1,y1) ,so..
return s[x2][y2]-s[x2][y1]-s[x1][y2]+s[x1][y1];
}
int slove(int x1,int y1,int x2,int y2,int k){
int &ans=f[x1][y1][x2][y2][k];
if(k==1) return S(x1,y1,x2,y2)*S(x1,y1,x2,y2);//actually there is no need to calc it twice
if(ans) return ans;
ans=1<<30;//big enough!
if(x1<x2)
for(int i=x1;i<x2;i++){
int t1=slove(x1,y1,i,y2,k-1);
int s1=S(i+1,y1,x2,y2);
int t2=slove(i+1,y1,x2,y2,k-1);
int s2=S(x1,y1,i,y2);
ans=min(ans,min(t1+s1*s1,t2+s2*s2));
}
if(y1<y2)//maybe useful?
for(int i=y1;i<y2;i++){
int t1=slove(x1,y1,x2,i,k-1);
int s1=S(x1,i+1,x2,y2);
int t2=slove(x1,i+1,x2,y2,k-1);
int s2=S(x1,y1,x2,i);
ans=min(ans,min(t1+s1*s1,t2+s2*s2));
}
return ans;
}
int main(){
cin>>n;
for(int i=1;i<=8;i++)
for(int j=1;j<=8;j++)
cin>>a[i][j];
for(int i=1;i<=8;i++)
for(int j=1;j<=8;j++)
s[i][j]=s[i-1][j]+s[i][j-1]-s[i-1][j-1]+a[i][j];
cout<<slove(1,1,8,8,n);
return 0;
}